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I. INTRODUCTION 


The Space Tug Thermal Control Follow-On program was conducted to 
further explore some of the thermal control concepts proposed for 
use in Space Tug in a breadboard test program. The objectives of 
the test program were to demonstrate the thermal control capabil- 
ities of a louver/battery configuration and a thermal conditioning 
panel/heat pipe radiator configuration. An additional objective 
was added to model the header pipe and radiator of the second 
test and correlate the analysis with the test results. 


These three objectives were achieved and are reported separately 
within the three appendices of this report. 


II. 


SUMMARY 


The louver /battery test was conducted to demonstrate the louver’s 
ability to control the Tug battery temperature and also compare 
actuator mechanisms. Two flight type louver assemblies were pro- 
cured and tested. Each of these assemblies were actuated by 
bimetallic springs. Results of this test were compared with a 
Martin Marietta built louver assembly which was actuated by a 
fluidic actuator. Each unit exhibited similiar effective em- 
ittances and performance when full open and closed. The bimetallic 
units function well throughout the testing while fluidic actuator 
exhibited a sluggish response. 

The heat pipe radiator/thermal conditioning panel test was conducted 
to demonstrate the performance of an all heat pipe thermal control 
system. The radiator and thermal conditioning panel were provided 
GFP for this test. The two units were connected together thermally 
via a procured fixed conductance heat pipe. Steady-state and 
transient tests were performed with all test articles performing 
flawlessly. The major conclusion drawn from this test was that 
physical interfaces between the test articles should be minimized. 
Bolted and clamped joints for such a design result in Ats that 
reduce the system capability to reject heat. The test demonstrated 
the radiators capability to reject the design heat loads at design 
environments and also demonstrated the capability to reject heat 
loads at lower than the design cold case environment. Appendix B 
discusses the test program and results. 

The development of an integrated thermal model of the radiator, 
header pipe, and thermal conditioning panel provided the means to 
evaluate the above test results. The analysis of heat pipe sys- 
tems in steady-state and transient conditions has recently been 
a state-of-the-art problem. While the results of the analysis 
indicate that system performance can be obtained with confidence, 
the analysis is sensitive individual pipe performance and inter- 
face conductances. Appendix G presents the analysis results and 
overall conclusions. 
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I. 


INTRODUCTION 


Application of a louver system to thermal control of typical space 
tug components is desirable for many reasons . They are partic- 
ularly useful because of high reliability, a semi-passive nature 
requiring no electrical power, and their ability to readily com- 
pensate for a change in the component power dissipation on the 
spacecraft environment. 

Interior louvers are even more reliable and predictable than ex- 
posed louvers because they are unhampered by the complex effects 
of solar impingement and direct exposure to the space environment. 
However, the heat rejection is reduced by the addition of the 
cover which acts as a radiation shield to the environment. The 
choice to use covered louvers on the tug was based on the thermal 
and mission requirements. The heat dissipation of typical tug 
components, with the exception of the fuel cell, are within the 
control range of covered louvers with moderate surface areas. 
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IX. SUMMARY 


Two louver systems were tested to demonstrate their thermal control 
capabilities as applied to spacecraft components. Both louver 
systems are of the interior type and differ primarily in the way 
the blades are actuated. System 1 uses bimetallic spiral wound 
springs while System 2 uses a fluidic actuator developed for the 
Viking thermal switch. 

The primary objective of both series of tests was to demonstrate 
the louver systems capability of maintaining thermal control of 
typical spacecraft components. Both steady-state and transient 
environmental conditions were simulated along with various compo- 
nent power dissipations. A post test model correlation was then 
carried out to verify the analytical techniques and to establish 
interior louver system operational characteristics. 
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III. 


TEST DESCRIPTION AND OBJECTIVES 


Two covered louver systems were tested in a thermal vacuum en- 
vironment in order to establish the operational characteristics 
of covered louver systems and to verify the analytical modeling 
techniques. The primary difference between the two louver sys- 
tems tested was the way in which the blades were actuated. Louver 
System 1 used the more conventional bimetallic spiral wound actu- 
ators while Louver System 2 used a fluidic actuator built by 
Martin Marietta for use on the Viking thermal switch. 

Each louver system was subjected to a series of tests designed 
tb demonstrate the capability of the louver system to maintain 
the thermal control of a simulated tug component. Both steady- 
state and transient conditions were simulated. Steady-state 
tests were performed by supplying a constant power input to the 
component simulator while maintaining the louver cover at a fixed 
temperature. Transient simulation consisted of varying the power 
supplied to the lower cover consistent with an absorbed environ- 
mental flux typical of a tug earth orbital mission. 


IV. 


BIMETALLIC ACTUATED LOUVER SYSTEM TEST ARTICLE 


Photographs of the bimetallic louver system taken during the test 
article fabrication are shown in Figures A^l through A-4. Table 
A-1 gives physical characteristics of the test article. The louver 
system is composed of two 20.32x40.64 cm (8x16 in.) flight-type 
thermal controller modules designed and built by Northrop Elec- 
tronics Division, Palos Verdes Peninsula, California. 

Each module consists of 16 highly polished aluminum blade halves 
bonded together to form eight louver blades. The eight louver 
blades are controlled in place by four bimetallic actuator ele- 
ments housed in the center channel of the lower system. The shaft 
end of each louver blade fits into a nylon spool in the center of 
each bimetallic element. The elements respond to a temperature 
change of the baseplate via a radiation coupling through slots 
in the bottom of the housing. The actuator is designed to pro- 
duce 90° of blade rotation over the temperature range of 288. 9°K 
to 302. 8°K (60°F to 85°F) . 

The modules were fastened to a 45.72x45.72 cm (18x18 in.) base- 
plate using low conductance nylon screws. The baseplate was 0.32 
cm (1/8 in.) aluminum. Two aluminum blocks weighing a total of 
6.41 kg (14.14 lb) provided the thermal mass for the component 
simulator. The blocks were bolted directly to the baseplate using 
a thermal heatsink compound to ensure a good contact. Heater tape 
was attached to the blocks to simulate an evenly distributed 
component power input. A fluid loop utilizing GH 2 was also at- 
tached to the blocks to provide for cooling the test article. 

This was necessary to establish initial temperatures at the start 
of each test. Guard heaters were put on the copper tubes to con- 
trol the heat flow from the test article during testing. 

The cover also measured 45.72x45.72 cm (18x18 in.) and was made 
from 0.16 cm (1/16 in.) thick aluminum. It was stood off from 
the baseplate using six teflon blocks measuring 1.27 cm (% in.) 
square by 6.35 cm (2H in.) long. Heater tape was used to provide 
up to 100 W of power to control the cover temperature. 

The baseplate and both sides of the cover were painted with 3 M 
Black Velvet paint (e ;^:;i0.90) . The entire assembly, with the 
exception of the cover, was wrapped with 20 layers of aluminized 
Mylar, each being separated by a layer of netting. 

Instrumentation used on the assembly consisted of 24 Type E chromel/ 
constantan thermocouples located as shown in Figure A-5. These 
thermocouples are tabulated in Table A-2. In addition to temper- 
ature measurement, voltage and current measurements were made for 
each of the heaters. These included the simulated component heater, 
the cover heater, and the fluid loop guard heaters. A photograph 
of the louver assembly mounted in the chamber is shorn in Figure A-6 . 
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Louver System Attached to Baseplate 






Figure A-Z Component Simulator Blooks and Fluid Loop 
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Figttpe A-4 Completed Bimetallic Actuated Test Article with Cover Installed 




Table A-1 

Summary of Bimetall-ia Actuated Louver Test 
Article Characteristics 


Item 

Characteristics 

Louver System 

Two Northrop flight-type thermal controller 
modules 

Dimensions 

20.3x40.6x4,9 cm (8,00x16.00x1.93 in.) 
per module 

Blades 

Polished aluminum, e 0.06, 8 per module 

18.3x4.9 cm (7.20x1.93 in.) 

Actuators 

Spirally wound bimetallic element; four 
per module 

Weight 

0.27 kg (0.6 lb) per module 

Baseplate 

Aluminum 

Dimensions 

45.7x45.7x0.32 cm (18.0x18.0x0.125 in.) 

Coating 

3M Black Velvet paint, e ^ 0.90 

Area 

Total 0.209 m2 (2.250 ft^) 

eff. Radiating 0.149 m^ (1.602 ft^) 

Weight 

1.84 kg (4.05 lb) 

Component Simulator 

Two aluminum blocks 

15 .2x30.5x2 .5 cm (6.0x12.0x1.0 in.) each 

Weight 

3.495 kg (7.705 lb) each 

Cover 

Aluminum 

Dimension 

45.7x45.7x0.16 cm (18x18x0.063 in.) 

Coating 

3M Black Velvet paint, e 0.90 (both ^ 

sides) 

Insulation 

20 Alternate layers of perforated alum- 
inized Mylar and tissue* glass 









. 

Table A-2 


Bimetallio Actuated Test Article Thermocouple Locations 

T/C Number 

Location 


1 

Baseplate 


2 

Baseplate 


3 

Baseplate 


4 

Baseplate 


5 

Frame Side (Inside) 


6 

Frame Side (Inside) 


7 

Frame End (Inside) 


8 

Frame End (Inside) 


9 

Actuator Center Channel 


10 

Actuator Center Channel 


11 

Component Simulator 


12 

Component Simulator 


13 

Baseplate underneath Actuator Element 


14 

Baseplate underneath Actuator Element 


15 

Cover 


16 

Cover 


17 

Cover 


18 

Cover 


19 

Guard Heater 


20 

Guard Heater 


21 

Guard Heater 


22 

Guard Heater 


23 

Outside Insulation 


24 

Outside Insulation 






Figure A-6 BimtalHo Actuated Test Article Mounted in Chamber 
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FLUID ACTUATED LOUVER SYSTEM TEST ARTICLE 


The fluid actuated louver system test article is shown in Figure 
A-7 which was a photograph taken during disassembly of the test 
article after completion of the test. The radiation character- 
istics of both louver systems are quite similar. That is, they 
have the same number of blades, approximately the same radiating 
area, and the same optical properties. The major distinction 
between the two systems is the blade actuation mechanism. Table 
A-3 gives the physical characteristics of the test article. The 
actual louver system was built by Martin Marietta in 1970 and was 
based on a louver design used on the Mariner spacecraft. 

A photograph of the fluid actuator used for the test is shown in 
Figure A-8. The actuator was developed by Martin Marietta for 
the Viking program to actuate the Viking lander thermal switch. 

A schematic of the actuator is shown in Figure A-9 . The rack 
and pinion gear mechanism which controls the louver blade rota- 
tion is powered by the actuator. The actuator utilizes the vapor 
pressure of Freon-12 as the pressure agent. For a given tempera- 
ture, the vapor pressure of Freon is constant and thus permits 
the actuator bellows to stroke without pressure loss due to volume 
increase. A single rod transfiits the actuator motion to the 
rack and pinion mechanism. The actuator is bolted through foot- 
pads that are integral to the actuator body to the component 
side of the louver baseplate. In addition, thermal grease was 
applied to the bolf,ed joints to provide good heat transfer from 
the baseplate. The actuator body is a stainless steel pressure 
chamber which is divided into two sections by metallic bellows. 

The bellows serve as flexible barriers to retain the Freon and 
were sized to provide the spring force to close the thermal con- 
tacts on the Viking thermal switch. 

The actuator was designed to transmit motion to the output shaft 
over a temperature range of 275°K (34°F) to 286. 1°K (55“F) . The 
desirable operating range for the louver system was 288. 7°K (60°F) 
to 302.5°K (85°F) which would require that the activation tem- 
perature range of the actuator be modified. However, after eval- 
uating the modification alternatives available, which included 
introducing a noncondensible pressurized gas into the bellows or 
adding an external pressurization mechanism to the housing, the 
decision was made to test the assembly at the design temperature 
range of the Viking actuator. Bench tests of the actuator by 
itself and attached to the louver system were performed to estab- 
lish the relationship between actuator temperature and output 
shaft stroke. The results of this test, shown in Figure A-10, 
indicate a very linear relationship between the actuator temper- 
ature and output shaft movement. 
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Table A--Z 

Summary of Fluid Actuated Louver System Test Article Characteristics 


Louver System 


Dimensions 

Blades 


Blade Coating 


Actuator 


Weight 

Baseplate 


Weight 

Coating 

Component Simula toi 
Dimensions 
Weight 
Cover 

Dimensions 

Coatings 

Insulation 


Characteristics 

Basic Mariner louver system design built by 
Martin Marietta and modified to allow blade 
positioning by a rack and pinion gear mech- 
anism. 

36.0x48.6x4.8 cm (14.7x19.125x1.90 in.) 

16 Aluminum-type H 45-1119 
20.96x4.67 cm (8.25x1.84 in.) 

Aluminized pol^rster tape . 3M Company 
designation Y8437, E 0 . 06 

Viking thermal switch actuator built by 
Martin Marietta. Weight 0.91 kg (2.0 lb), 
dimensions 38 cm long, 3.81 cm dia (15x1.5 
in . ) . 

1.23 kg (2.71 lb) 

Aluminum 

Total 0.195 m? (2.01 ft^) 

Eff. Radiating 0.165 m^ (1.78 ft^) 

1.64 kg (3.62 lb) 

3M Black Velvet paint, e{^0,09 

Aluminum block with integral fluid loop 

48.5x36.8x1.3 cm (19.1x14.5x0,5 in.) 

4.88 kg (10.75 lb) 


46.8x36.8x0.16 cm (19.125x14.5x0,063 in.) 

3M Black Velvet paint, e?s0,90 both sides 

20 alternate layers of perforated aluminized 
Mylar and tissue glass 







Figure A 
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The louver system instrumentation is shown in Figure A-11 and 
tabulated in Table A-4. In addition to the test article temper- 
ature and power measurements, a potentiometer was attached to an 
extra pinion gear to record the blade rotation produced by the 
actuator. The paint coating and insulation system were identical 
to those used on the bimetallic louver system test article. An 
aluminum plate with an integral fluid loop was fabricated to 
provide both the simulated component mass and enable cooling to 
establish initial temperatures. A photograph of the test article 
mounted in the chamber is shown in Figure A-12. 
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A T/Cs Located on Blades 

T/Cs 4, 6, 11, 12 on Guard Heaters 
T/Cs 21 through 24 on External Cover 





TC Number 

Location 

1 

Baseplate 

2 

Baseplate 

3 

Baseplate 

4 

Baseplate 

5 

Guard Heater 

6 

Guard Heater 

7 

Frame End (Inside) 

8 

Frame End (Inside) 

9 

Frame Side (Inside) 

10 

Frame Side (Inside) 

11 

Guard Heater 

12 

Guard Heater 

13 

Louver Blade 

14 

Louver Blade 

15 

Outside Insulation 

16 

Outside Insulation 

17 

Comp on en t S imula tor 

18 

Component Simulator 

19 

Baseplate underneath Actuator 

20 

Actuator 

21 

External Cover 

22 

External Cover 

23 

External Cover 

24 

External Cover 









Figure A-12 Fluid Actuated Test Article Mounted in Chamber 
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‘TEST FACILITIES AND DATA FECOEDERS 


Photographs of the test articles mounted in the 6x15 thermal 
vacuum chamber are shoxm in Figures A-6 and A-12. The test arti- 
cles are suspended by low conductance stainless steel wire at- 
tached to the teflon cover mounting blocks. Guard heaters and 
insulation were used to prevent heat leaks from the inlet and 
outlet tubes to the fluid loops. 

A description of the thermal vacuum facility is presented in 
Table A-5. A photograph of the data acquisition systems is shown 
in Figure A-13. A Bristol strip chart recorder with a tempera- 
ture range from 61®K (-350®F) to 672°K (750°F) was used to re- 
cord thermocouple data. This wide temperature range was neces*^ 
sary because the next standard recorder scale, 144°K (-200°F) to 
478°K (400°F) , would not be adequate to record temperatures on 
the outer insulation. The accuracy of the temperature readings 
was ± 1.4°K (±2.5®F). Heater power which included two guard 
heaters, a cover heater, and a baseplate heater for both test 
articles was individually controlled using Varac controllers. 
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Table A-5 6x15 ft Vaoum Chamber Facility 


Nominal Chamber Size 

6 ft in diameter by 15,5 ft long 

Maximum Test Specimen 
Size 

5 ft in diameter by 11 ft long 

Ultimate Pressure 

1x10 ^ torr 

Pumping System 

Roughing - 

Portable cart capable of pumping the 
chamber to 5x10 ^ torr. Cart contains 
a cold trap and a 650 CFH booster pump 
backed up with a 60 CFM compound pump. 


Hi-vacuum *- 

Ion/ titanium sublimation system with net 
pumping speed in excess of 24,000 liters/ 1 
sec 

Cold Wall 

LN 2 cooled - 

99% optically dense shroud capable of 
absorbing a heat load of 26,000 Btu/hr 

Solar Simulation 

Target volume - 

1 sun - 42 in, diameter by 60 in. long 

2.5 sun - hexagonal, 30 in. across flats 
by 30 in. long 

Intensity ~ 

1 sun - 1400 W/m^ nominal with a range 
of 1250 W/m^ to 1600 W/m^ 

2.5 sun - 3500 W/m^ nominal with a range 
of 2000 W/m^ to 4500 W/m^ 









H[%s 



1 ^ 






»jl A\^%iVA'Vvii %>. . 


VII, THERMAL MODEL DESCRIPTION 


A therinal network model was constructed for each louver system 
test article for use in the MITAS (Ref 1) thermal analyzer. The 
models were used for both steady state and transient model cor- 
relations. The effective emittance of the baseplate as a func- 
tion of blade angle was calculated using a computer program (Ref 
2) which characterizes a covered louver system with a diffuse 
base and specular blades. The effective emittance is the radia- 
tion interchange factor from the baseplate to the cover and takes 
into account the blockage of the louver blades. In the case of 
the bimetallic louver system, the effective emittance valves were 
adjusted to account for the holes in the frame to lighten the 
assembly. 

The louvered area of each baseplate was modified to account for 
the blockage of the actuator housing and the frame. The model 
also simulated radiation from the baseplate to the cover through 
the unlouvered portion of each baseplate. Conduction through the 
teflon standoffs as well as heat leaks through the multilayer 
insulation were modeled. 

Test data input to both models consisted of the power applied to 
the component simulator heater and the cover temperature. The 
cover temperature was used as a boundary condition. However, 
the heat input necessary to hold the cover at the test data tem- 
perature was calculated and compared to the cover power test data. 
An array of blade angle versus baseplate temperature taken from 
Northrop data for the bimetallic louver system was input into the 
model. A blade angle, measured directly from the potentiometer, 
versus time array was put in the model for the fluid actuated 
test article. 
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VIII. TEST RESULTS AND MODEL CORRELATION 


The results of the testing and model correlation for the steady 
state test using the bimetallic louver system test article are 
summarized in Table A-6. The table gives the heater power ap- 
plied to the simulated Component the final test data for the 
cover* power and temperature, the final baseplate temperature test 
data reading, the model baseplate temperature correlation, as well 
as the figure number which shows the model correlation and 
data for the baseplate temperature versus time. As shown in the 
plots of the baseplate temperature history, the actual test data 
was tabulated from approximately two to five hours depending on 
the duration of the particular test. The last test data value for 
the cover temperature was used to project the transient simula- 
tion out to 10 hours at which time a ttde- steady state solution 
was simulated. The extended transient simulation indicates the 
steady state temperature which would have been obtained had the 
actual test been allowed to continue. 

The baseplate temperature results of the transient tests using 
the bimetallic louver system are shown in Figures A-21 through 
A-23 for 0, 20, and 45 W of component power respectively. The 
power supplied to the cover corresponds to an absorbed environ- 
mental flux (Ref 3) given in Table A-7. \ 


A- 2 7 


Table A-6 

Siomavy of Steady State Tests Using the Bimetatlia Actuated 
Test Article 


Component 

Power 


W 


15.0 

15.0 

30.0 
30.0 

30.0 

45.0 
45.0 


W 


43.3 

57.7 

14.3 
28.2 

43.7 

0.0 

13.1 


“K 


273.3 

287.6 

254.4 
272.2 

288.6 
254.4 
272.1 


Baseplate Temperature 


Model 


°K 


(“F) 


297.6 ( 75.7) 

304.6 ( 88.2) 

302.7 ( 84.8) 

310.5 ( 98.9) 

318.8 (113.8) 

313.2 (103.7) 

321.3 (118.4) 


Baseplate 

Temperature 


Figure No. 


Table A-7 Cover Heatvng Environment Used for Transient Cases 


Time, hr 

Absorbed Flux 

W/ cm^ 

(Btu/hr-ft^) 

0.000 

1.695E-02 

(53.767) 

0.135 

1.079E-02 

(34.203) 

0.269 

5.139E-03 

(16.294) 

0.404 

2.002E-03 

( 6.348) 

0.538 

7.951E-04 

( 2.521) 

0.592 

6.069E-04 

( 1.924) 

0.592 

6.069E-04 

( 1.924) 

0.673 

6.135E-04 

( 1.945) 

0.807 

1.712E-03 

( 5.429) 

0.942 

4.269E-03 

(13.537) 

1.022 

6.557E-03 

(20.792) 

1.022 

6.557E-03 

(20.792) 

1.076 

1 

8.396E-03 

(26.620) 

1.211 

1.299E-02 

(41.214) 

1.345 

1.811E-02 

(57.407) 

1.480 

1.998E-02 

(63.349) 

1.608 

1.696E-02 

(53.767) 


Average 8 . 467E-03 

(26.842) 
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Figure A-14 

Bimetallia Actuated Test Article Baseplate Temperature for 
IS W Power, 273. 3° K Cover 
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figure A-15 

Bimetallic Actuated Te^t Article Baseplate Temperature for 
15 ^ Power j 28?.6''K Cover 
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Figure A-16 

BimetdlliG Actuated 'lest Article Baseplate Temperature for 
ZO W Bower ^ 264, 4°K Cover 
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Figure A-19 ^ ’ I 
Bimetallic Actuated Test Article Baseplate Temperature for ^ ] 
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Figure A-20 

Bimetaltic Actuated Test AvticZe Baseplate Temperature for 
45 W Power^ 272. VK Cover 
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The results of the steady state cases using the fluid actuated 
louver system test article are summarized in Table A-8. The base- 
plate temperature model correlation and test data plots extend 
only to the end of the actual test. Further simulation using the 
model (past the duration of the test) could not be performed be- 
cause test data only was used for the blade position. An analy- 
tical model of the actuator position as a function of baseplate 
temperature was attempted based on the steady state curve fit of 
the actuator temperature versus stroke relationship shown in Fig- 
ure A-10. Hox^ever, in the covered louver system test, even though 
the actuator housing temperature indicated a particular blade 
position based on Figure A-10, the measured blade position would 
not correspond to this value. This was probably due to a dif- 
ference between the Freon vapor temperature and the housing tem- 
perature because of the relatively inefficient gas conduction 
mechanism necessary to head the Freon vapor. Once the vapor tem- 
perature started to respond to the change in housing temperature, 
the housing temperature had exceeded the control range and re- 
sulted in driving the actuator from one extreme to the other. 

The influence of cover temperature on the baseplate was minimized 
by not applying power to the cover in an attempt to gain control 
of the actuator. 


Table A-8 Simnary of Steady State Tests Using the Fluid Actuated Test Article 


Component 

Power 

Cover 

Power 

Cover 

Temperature 

Baseplate 

Temperature 

Baseplate 

Temperature 





Est 

Model 


w 

w 

°K 

(° F ) 

“K 

(°F) 

“K 

(»F) 

Figure No. 

0.0 

45.0 

261.8 

( 11.3) 

264.2 

( 15.5) 

265.3 

(17.6) 

A-24 

5.0 

0.0 

170.1 

(-153.8) 

249.6 

(-10.8) 

251.7 

(-7.0). 

A-25 

7.5 

0.0 

172.1 

(-150.3) 

256.8 

( 2.3) 

258.0 

( 4.4) 

A-26 

10.0 

0.0 

190.6 

9-117.0) 

260 . 6 

( 9.0) 

262.6 

(12.6) 

A-27 

15.0 

0.0 

199.6 

(-100.8) 

273.5 

( 32.3) 

272.3 

(30,2) 

A-28 

20.0 

0.0 

210.3 

(- 81.5) 

277.4 

( 39.3) 

277.5 

(39.5) 

A-29 

30.0 

0.0 

228,9 

(- 48.0) 

283.3 

( 50.0) 

281.8 

(47.2) 

A-30 

35.0 

0.0 

238.9 

(- 30.0) 

292.8 

( 67.0) 

289.2 

(60.5) 

A-31 


A-39 
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Fluid Actuated Test Avticle Ba 
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Two transient tests were simulated for the fluid-activated louver 
system. Case 1 involved continually varying the input power to 
the simulated component with no power to the cover* The baseplate 
temperature model prediction and test data for this case is shown 
in Figure A-32. The power to the base, blade angle, and cover 
temperature all versus time are shown in Figures A-33, A-34, and 
A-35 respectively. Transient Case 2 involved applying a constant 
10 W of power to the simulated component and imposing the time 
varying environmental flux given in Table A-7 to the cover. Base- 
plate temperatures for Case 2 are given in Figure A-36. The blade 
angle and cover temperature test data for this case are given in 
Figures A-37 and A-38 respectively. 
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Figure A-33 

Fluid Actuated Test Article Transient Component Power for Case 1 



Figure A-34 

Fluid Actuated Test Article Transient Blade Angle for Case 1 
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Figure A-36 

Ftuid Actuated Test Article Transient Baseplate 
Temperature for Case 2 
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IX. CONCLUSION 


The results of both louver system steady state and transient tests 
indicate the condition when component thermal control can be main- 
tained. The analytical models developed to correlate the test 
data are capable of predicting the operational characteristics of 
an interior louver system. The bimetallic actuated louver system 
operated as expected. It was able to maintain component tempera- 
tures below 305 (90°F) while dissipating 15 W of power to a 
cover at a temperature of approximately 289°K (60°F) or dissipating 
30 W to a cover at approximately 256°K (0°F) . Also there was no 
detect'able time lag between component baseplate temperature and 
the blade response. 

The operation of the fluid actuated louver system test article 
was affected by the actuator design. The desirable operating tem- 
perature range of the louver system had to be lowered to be com- 
patable with the actuator temperature range. This had an adverse 
affect on the capability of the louver system to reject heat 
because the radiation heat transfer mechanism depends on the ab- 
solute temperature I'^vel of the radiating bodies. This can be 
easily corrected in the actuator design by either using a different 
fluid or by pressurizing the existing Freon actuator with about 
10 psi of GN 2 or He, A more critical characteristic of the ac- 
tuator was the significant time lag involved for the actuator to 
respond to a temperature change of the component. This was at- 
tributed to the relatively slow heating of the Freon actuator vapor. 
This did not effect the original Viking thermal switch application 
of the actuator because the response to temperature change was not 
a critical factor. However, when response to temperature change 
is critical, such as the louver system, then this type of actuator 
may not be adequate. 
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I . INTRODUCTION 


The application of heat pipes to thermal control of satellites 
has grown significantly in recent years. In most cases the heat 
pipe was used only as a transport medium. This test was intended 
to demonstrate the ability of an integrated heat pipe system to 
collect, transport, distribute, and control heat from the source 
to the sink. This method can result in significant weight reduc- 
tions while improving the system performance and reliability if 
the heat pipes perform as required. 

The ability to control a system of this nature requires the use 
of variable conductance heat pipes (VCHP) . The VCHP has been, 
and still remains, the pacing item for the development and flight 
of an integrated heat pipe system. Configuration (a) of Figure B-1 
was tested. 



Figure B-1 Ir^ egrated Heat Pi-pe Systems 


Each of the systems shown in Figure B-1 uses a collector panel 
called a thermal conditioning panel. This panel has integral 
heat pipes enabling the application of nonuniform heating while 
controlling the temperature differential over the panel to a 
small value. For low heat loads, this panel would provide three 
of the four functions presented above. 
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Additional rejection capability can be achieved with large heat 
loads by adding a separate radiator with integral control capa- 
bility in the radiator feeder pipes, in the header pipe, or both. 

Our original plan was to test configuration (c) (Figure B-1) ; 
however, performance problems were experienced with the government 
furnished property (GFP) header pipe. Hence, configuration (a) 
was tested. We accomplished this by procuring a fixed conductance 
header pipe. To facilitate installation into the chamber, the 
header pipe had a 90° bend in the adiabatic section of the heat 
pipe. 



VuA 1 A. .1.1 \ - >1 iL4.(. 


SUMMARY 


Testing of the heat pipe radiator/ thermal conditioning panel was 
conducted in the 6x15 thermal vacuum chamber at ^rtin Marietta’s 
Space Simulation Laboratory. Three series of tests were conducted 
to demonstrate the system capability. This was the first series 
of vacuum testing of the radiator which performed flawlessly 
throughout the test program. The design of the radiator and the 
thermal conditioning panel was limited to 324. 8 °K because ammonia 
was used as the working fluid in each unit. This resulted in the 
test plan being modified to perform an additional series of test'- 
ing to further demonstrate the radiator capabilities. 
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III. 


OBJECTIVE 


The overall objective of the Space Tug Thermal Control Follow-On 
contract was to experimentally demonstrate that the thermal con- 
trol subsystem concepts selected and defined in Phase I of the 
contract were viable, to obtain data for refining analytical 
models, and to optimize concept configurations and performance. 

The objective of this test was to demonstrate and measure the 
performance of an integrated heat pipe thermal control system. 
The special test objectives \<feTe to: 

1) Measure performance of the system under a constant thermal 
environment on the radiator at constant heat loads to the 
thermal conditioning panel to obtain steady-state results. 
Hot and cold radiator environments were tested. 

2) Measure the performance of the system under a cyclic thermal 
environment with a constant heat load to the thermal condi- 
tioning panel. 
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IV. TEST DESCRIPTION 


These objectives (Sec. Ill) were accomplished by conducting the 
testing consistent with the radiator and thermal conditioning 
panel design requirements. The upper temperature limit of 
324. 8°K required a modification to the original test plan to 
avoid over temperature of the thermal conditioning panel. 

Table B-1 summarizes the tests that were conducted. Table B-2 
presents the original test plan summary. 


Table B-1 

Thevmal Conditioning Panel/Eeat Bi-pe Radiator Test Siamary 




Testing Conducted 



1 

2 

3 

4 

Type of Test 

Transient 

Transient 

Cyclical 

Transient 


to Steady- 

to Steady- 

1.5 hr 

to Steady- 


State 

State 

Period 

State 

Conditioning 

100 

200 

200 

600 

Panel Heating 

200 

400 

400 

400 

Level, W 

300 

500 


500 


380 





200 



. 

Radiator Sink 

Constant 

Constant 

Cyclical 

Cons tant 

Temperature, °K 

238.66 

194.29 

238.7 

144.29 




to 194 


Simulated Flight 

Hot 

Cold 

Orbital 

Colder 

Environment 

^ i 

Environment 

Environment 

Cyclical 

Environment 





Table B-2 

'Planned Thevmal Condit'ioning Panel/Eeat Pipe Radiator Test Sumary 



Test 


1 

2 

3 

1. Type of Test 

Transient to 

Transient to 

Cyclic, with 


Steady-State 

Steady-State 

1. 5-hr’ Period 

2 . Conditioning 

100 

200 

200 

Panel Heating 

200 

400 

400 

Level, W 

300 

600 

600 


400 

800 



200 

1000 


3. Radiator Sink 

Constant 

Constant 

Cyclic, 

Temperature, °K 

238.66 

194.29 

238.7 to 194 

4. Simulated Flight 

Hot 

Cold 

Variable 

Environment 

Environment 

Environment 

Environment 


Orbit 

Orbit 

Orbit 


The first two series of tests were conducted at steady-state 
shroud temperatures of 238. 7 °K and 194. 3 °K while varying the power 
to the levels shown until the temperatures achieved steady-state. 
The third series of testing was conducted at constant power levels 
while varying the shroud temperature in a sine wave manner between 
238.7 and 194. 3®K with a 90 min period. The fourth series of 
testing was conducted at a shroud temperature of 144 -3°K at 
three power levels. 


TEST ARTICLE DESCRIPTION 


Heat Pipe Radiator 

The radiator panel (Figures B-2 and B-3) is an aluminum panel com- 
prised of six variable conductance feeder heat pipes. The. evap- 
orator section of each pipe is bent 90° from the condenser to 
permit attachment to the header pipe. The radiator is 1.22 by 
1.83 m (4 by 6 ft) with a total radiating area of 2.23 m^ (24 ft^) 
excluding the heat pipe and reservoir areas which extend beyond 
the radiator panel. The feeder pipes used ammonia as the working 
fluid and nitrogen as the noncondensable control gas. 
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Figure B-2 
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The radiator panel was developed for NASA-MSFC under contract 
NAS8-29905 by Grumman Aerospace Corporation, The radiator was 
then supplied GFP to Martin Marietta for this test program. 
Reference 1 describes the design^ build > and test of the radial- 
tor. 

The radiator panel was received with 38 copper-cons tantan thermo- 
couples installed on the radiator. An attempt was made to use 
these thermocouples but proved impractical due to the facility 
instrumentation capabilities and cost of supporting hardware. 

The laboratory is setup to use chromel-constantan thermocouples. 
The use of copper cons tan tan would have required the modification 
of the recorders, chamber feedthroughs, connecting wiring, etc. 

The modifications and unnecessary cost items associated with using 
the existing thermocouples were avoided by installing thermo- 
couples compatible with the existing instrumentation system. 

The radiator was instrumented with chromel-constantan thermo- 
couples (Figure B-4) . The thermocouple locations agree with 
the present thermocouple locations with some deletions and addi- 
tions. A total of 34 thermocouples were installed. Table B-3 
presents, for information purposes, a comparison of the thermo- 
couples installed by Grumman versus those installed for this 
series of tests. 


Table B-3 

Thermocouple Comparison 


GAC 

MM 

GAC 

MM 

GAC 

MM 

GAC 

MM 

1 

NI 

11 

5 

. 21 

15 

31 

27 

2 

NI 

12 

NI 

22 

17 

32 

28 

3 

1 

13 

7 

23 

18 

33 

29 

4 

2 

14 

8 

24 

19 

34 

30 

5 

NI 

15 

9 

25 

20 

35 

31 

6 

3 

16 

10 

26 

21 

36 

32 

7 

NI 

17 

11 

27 

22 

37 

33 

8 

NI 

18 

12 

28 

23 

38 

34 

9 

4 

19 

13 

29 

24 



10 

NI 

20 

14 

30 

25 



Notes 

: 1) 

NI - 

Not installed 




2) 

Thermocouples 6, 

16, and 26 

do 



not have corresponding 

GAC 




thermocouples . 
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After installing the thermocouples, an insulation blanket was 
assembled from 20 layers of vented aluminized Mylar with nylon 
netting separators. The blanket was sized to drape over the 
radiator sides and the reservoir end to form an enclosure with 
the flat plate shroud. A second 20-layer blanket was fabricated 
for completing the enclosure at the header pipe end of the radia- 
tor. 

Header Heat Pipe 

The header heat pipe was purchased from Isothermics Incorporated. 
This pipe (Figures B-5, B-6, B-7 and B-8) was designed to carry 
a 1-kW load using Genetron 31-114 es the working fluid. The 
pipe was made of copper with mounting provisions integral to the 
design. The pipe was designed for use in a 1 g environment and 
does not include a wick internal to the pipe. 

Fifteen thermocouples were installed on the header pipe as shown 
in Figure B-8. Thermocouples 10 through 15 were matched with 
radiator thermocouples 1, 6, 11, 16, 21, and 26 respectively. 

The header pipe was checked out at Isothermics before shipment 
to Denver. Attachment A presents the data resulting from that 
checkout. 

Thermal Conditioning Panel (TCP) 

The thermal conditioning panel was developed for NASA-MSFC by 
Donald W. Douglas Laboratories under contract NAS8-28639. The 
panel was government furnished property (GFP) supplied to Martin 
Marietta by MSFC. References 2 and 3 describe the design, build, 
and test of the panel. The panel was designed as an isothermal 
panel with 12 heat pipes between 2 aluminum panels. Ammonia was 
used as the working fluid. 

The thermal conditioning panel (TCP) was instrumented (Figure B-9) 
with 11 thermocouples. In addition, five heater zones of 200 W 
capability each were added (See Figure B-10) . All heater zones 
were controllable over the 200 W range. A heater tape 0,95 cm 
(0.375 in.) wide was applied to the TCP covering the entire area 
of each zone. 

The application of heat to the panel agreed with the recommended 
0.31 W/cm^ (2 W/in.^). Removal of the applied heat load at the 
header pipe interface was accomplished at a flux density of 0.66 
W/cm^ (5.56 W/in.^). Tliis appears to be consistent with previous 
testing conducted under the original build and test contract to 
Donald W. Douglas Laboratories. 

A 20-layer aluminized Mylar blanket was fabricated to provide an 
insulated enclosure for the TCP and the header pipe evaporator 
section. This blanket was installed after chamber installation 
of the TCP. 
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Figure B~6 Header Heat Pipe and Thermal Conditioning Panel 
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SIDE B 

Figure B-9 Thermal Conditioning Panel Thermocouple and Heater Zone Locations 














Environmental Shroud 


A flat shroud (Figure B-11) was used in this test as the radiator 
radiation sink. The shroud was 122x213 cm (48x84 in.) and was 
instrumented with 10 thermocouples (Figure B-12) . The shroud was 
required to enable meeting the 1.5«hr environmental temperature 
cycle. The chamber shroud would not respond fast enough over the 
required temperature range in the 1.5-hr cycle time. The flat 
shroud was insulated on the chamber side to avoid cooling down the 
present chamber shroud. This minimized the heat load on the gas 
heat exchanger system used to control the shroud temperature. 


TEST CONTROL AND INSTRUMENTATION 


The heater control of the test was accomplished by using an 
existing regulated dc power supply and control panel shown sche- 
matically in Figure B-13. Figure B-14 shows the control panel 
including digital potentiometers for controlling each zone. 
Voltage and current measurements were obtained by selecting the 
desired channel and recording the appropriate meter outputs. 
Figure B-15 shows the overtemperature controller on top of the 
control panel x^hich was wired to interrupt power to the heaters 
on the TCP when it exceeded 324. 8°K. Also shown are the multi- 
point recorders used in recording the data. The data were re- 
corded manually on summary sheets during the test. 

The number of thermocouples and heaters was quantized to be con- 
sistent with previous testing and to provide data consistent with 
the overall test objectives. Table B-4 summarizes the instrumen- 
tation requirements for each test article and the chamber shroud. 

Table B-^4 

Instrimentation Summary 


Item 

Thermo coup les 

^ ^ 

Heaters 

Radiator 

34 


Header Pipe 

15 


Thermal Conditioning Panel 

11 

5 

Shroud 

10 

— 

Total 

70 

1 

5 
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The thermocouple data were recorded on 24-channel recorders. The 
radiator data required two recorders and permitted recording the 
TCP temperatures on the second recorder. The header pipe and 
shroud data were recorded on separate recorders. The heater power 
was recorded manually on separate data sheets. The accuracy of 
thermocouple recording was ±1.4°K (±2.5°F). 


CHAMBER INSTALLATION 


The three subassemblies were installed in the chamber as shown 
schematically in Figure B-16. 

The radiator and thermal conditioning panel were suspended from 
the aluminum frame shown in Figure B-17. Figures B-2 and B-3 
show the teflon blocks used to suspend the radiator. The header 
heat pipe condenser mounting surface was greased with a DC 29 
heat sink compound before clamping to the radiator feeder pipes. 
The same compound was used to interface between the header pipe 
and the TCP before bolting the two units together (See Figure 
B-18) . 

The radiator, header pipe, and TCP were then leveled to ensure 
proper heat pipe operation. The shroud was installed 6.5 cm 
(4 in.) below the radiator completing the piping installation. 
After conducting a controls and instrumentation checkout, the 
radiator and shroud were enclosed by an insulation blanket com- 
prised of 20 layers of 1/2 mil-vented double aluminized Mylar 
with nylon netting separators. This created an enclosure 
between the radiator and the shroud and minimized the influence 
of the chamber upon the test. 

Figure B-19 shows the final installation with the TCP and line 
insulation installed. The radiator /shroud blanket is not visible 
in the figure. 
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Figure B-16 
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Figure B-16 Cl'icmber Installation Schematio 
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PveZiTniyioot*ij Chamber Installation 










V, 


TEST FACILITY 




A 6xl5“*ft ion/titanium sublimatioii'-puinped chamber capable of 
attaining a vacuum level of 10"® torr with a pumping speed of 
50,000 liters/sec was used during the conditioning panel/heat 
pipe radiator tests (See Figure B-20) . 

The operating characteristics of the 6xl5-ft solar-- thermal/vacuum 
chamber, based on both design specifications and actual test data, 
are shown in Table B-5. 

Table 5-5 

Operaixing ChaTaotevistlos of the SolaT-TheTmal/VaouiM Chamber 


Nominal Chamber Size 

Six ft O.D, X 15 ft 6 in. long 

Ultimate Pressure 

lxl0~® torr 

Pumping System 

Roughing 

Portable cart capable of pnmping the 
chamber to 5x10"^ torr. Chart contains 
a cold trap and 650-cfm booster pump 
backed up with a 60*-cfm compound pump. 


High Vacuum 

Bulk titanium subliminator with LN 2 
cooled substrate with an IntTinsto 
pumping speed of 50,000 liters/sec at 
1x10" torr. Ion pumps are capable 
of pumping 1400 liters/sec of air with 
a noble gas load capability of 20% 
relative to air for argon. 

Chamber Feedthroughs 

lus trumentation 

300 chromel/constantan thermocouples. 

' 

Power 

24 channels, each with a capability of 
28 Vdc at a maximum current of 10 amps. 
Eight feedthroughs with a 100-amp capa- 
city. 


Additional capabilities of the 6xl5-ft solar-thermal/vacuum 
chamber, which are available, but were not required to support 
the proposed test program, are shown in Table B-6, 
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Table B-6 

Additional Capabilit'ies of the Solar--Thermal/Vaout07i Chambev 


Solar Simulation 


Intensity 

One to 2.5 solar constants (1400 to 4000 
W/m^) . 

Uniformity 

±5% over a 4-cm^ area in a plane per- 
pendicular to the axis of the beam. 

Spectral Range 

0.5 to 2.7 p 

Spectral Match 

Filtered zenon with a olose match to the 
Johnson curve. 

Residual Gas Analysis 

Continuous scanning of the mass spectrum 


from 1 to 500 AMU or any portion thereof. 


Temperature control of the shroud was obtained by using the por- 
table gas heat exchanger shown in Figures B-21 and B-22. The 
unit uses gaseous nitrogen forced through a closed loop which is 
controlled to the desired temperature. Liquid nitrogen is used 
for cooling and a 10 kW heater is used for heating. An elec- 
tronic controller is used to control the system by manual or 
automatic setpoint inputs. The system is capable of controlling 
the output gas temperature between 116 and 339 °K (-250 and 150°F) 
while subjected to heat loads up to 10 kW. 

Standard Operating Procedures (SOPs) govern the operation and 
protection of both the test specimen and the test facility. 
Specific test requirements are integrated into these SOPs through 
a controlled release system. The facility is operated by certi- 
fied operators who have a thorough knowledge of vacuum technology 
and test techniques. 
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VI. 


TESTING 


The data presentation for all of the testing used the following 
nomenclature shown in Table B~7. 


Table B-7 Nomenolatuve Used for All Data Tresentation 


Thermocouple 

Designation 

Test 

Hardware 

Reference, 

Figure 

T 

Thermal Conditioning 
Panel 

B-IO 

H 

Header Pipe 

B-9 

R 

Radiator 

B-4 

S 

Shroud 

B-13 

Note: The shroud temperatures were presented as average values 

since the variation in temperatuj 
5°K. 

ce over the shroud was 


TEST 1, HOT CASE, 238.66°K SHROUD 


Testing of the radiator /thermal conditioning panel was performed 
from 16 April to 24 April 1975 beginning with the hot case. For 
this series of testing the shroud was held at 238. 7°K (-30°F) 
and heat was applied to the thermal conditioning panel (Table B-7) 
until steady-state was achieved. The planned 400 W test was 
reduced to the 380 W level to avoid over temperature of the thermal 
conditioning panel. The steady-state and transient data obtained 
from the tests are presented in Figures B-23 through B-52. 

Figures B-23 through B-32 summarize the test data by location on 
the test hardware for the five power levels at steady state condi- 
tions. Figures B-33 through B‘-52 present temperature histories 
from the same runs shox^ing the transient response of the system 
to varying power loads for constant environment' temperatures . 

The thermal conditioning panel exhibited a maximum of 8°K AT 
(Table B-7) over the panel. The header pipe exhibited ATs ranging 
from 3°K at 100 W to 8°K at 380 W. The two 200 W runs both ex- 
hibited a 5.5°K AT. 
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Note; Teinperature, °K 














Note: Temperature, ""K 












Figure B-27 





































Figure B-31 



Legend : 

• Thermocouple Locations 

Figure B-31 Test 1 Hot Case ^ 200.23 W 
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Figure B^36 
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Figure 8^43 
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Figure 8-43 Teet 1, Hot Case fheimal Conditioning Panel 
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Figure B-51 


Teet 1, Hot Case Radiator Skin Temperature 
















The radiator data clearly shows the front movement throughout 
this test series as exemplified by Pigures B-35 and B‘-45 by 
thermocouples R3 (Table B^-?) and R4. Thermocouple R2 was ob- 
served during the test to read lower than R3* This thermocouple 
(Figure B-4) was located on the adiabatic section of feeder tube 
1. After removing the radiator from the chamber it appeared 
that this thermocouple had become detached from the tube and was 
very likely influenced by the shroud temperature. The radiator 
skin temperature shown in Figures B-41 through B-51 show the 
expected trends including a nonuniform front location in the 
feeder tubes. This was anticipated since the evaporator section 
of feeder pipe 6 is only 15.24 cm (6 in.) long compared with the 
30.48 cm (12 in.) length of the other 5 pipes. 

The shroud temperatures are presented as average values in the 
transient plots and as individual temperatures by location in the 
steady-state summaries. Thermocouple S-10 (Table B-7) failed to 
yield data comparable with the other shroud temperature and is not 
presented in any of the data summaries. The average temperature 
plots are based upon the nine thermocouples for each plot. No 
attempt was made to correct the problem on the S-10 thermocouple 
due to the time required to remove the insulation blankets and 
the potential requirement of removing the shroud from the chamber 
before corrective action could be performed. After completing 
the test, this thermocouple was found to be detached from the 
panel. 


TEST 2 COLD CASE 194.29'^K SHROUD 


Cold case testing of the system was continued after completing the 
hot case by reducing the shroud temperature to (*"110 ° F) . 

The 200 and 400 W power levels were accomplished without difficulty 
after increasing the power to 600 W over the temperature of the 
thermal conditioning panel was achieved. At that point the test 
plan was modified to three power levels; 200, 400 and 500 W. The 
600, 800, and 1000 W tests x^ere deleted. In place of the two 
deleted tests a colder shroud temperature test was added by 
agreement with the COR, Mr. Jim Wise. 

The data derived from the cold case testing are presented in 
Figures B-53 through B-68. Figures B-5 3 through B-5 8 present 
the steady-state results by location on the test hardware for the 
three power level tests. Figures B-59 through B-68 present 
temperature histories from these runs showing the transient 
response from a changing power load at a constant shroud tempera- 
ture . 
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Figure B-5S 











Figia?e B-S6 Test 2 Cold Case, 398.27 W 










Figure B-57 







Note: Temperature, 
















Figure B-60 










Figure B^Sl Teet 2, Cold Case Radiator Beat Pipe Bo* 1 













reproducibility of the 
ORIGINAL PAGE IS POOR 


















Hgure B-66 Teet 2, CoU Com Radiator Beat Pipe No. 6 
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Figure B-55 
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The thermal conditioning panel exhibited a maximum of 9°K AT over 
the panel at 200 W. The 400 and 500 W runs both resulted in 8°K 
AT. This agreed with hot case results for the 400 W condition, 
however the 200 W results were 1®K higher. The reason for this 
discrepancy is unexplained. 

The header heat pipe exhibited ATs of 4, 7, and 9°K respectively 
for the power levels of 200, 400, and 500 W. These data are in 
agreement with the hot case results. 

The radiator temperatures show a pronounced movement of the vapor/ 
gas interface as shown in each of the feeder tube temperature 
plots (Figures B-61 through B--66) . 


TRANSIENT TESTS 


The transient tests were conducted by cycling the shroud tempera- 
ture between the hot and cold case shroud temperatures 238. 7 °K 
and 194. 2°K, in a cosine wave form with a 90 min period. Two 
transient tests were conducted as a result of the over tempera- 
ture problem experienced during the cold case test series. A 
minor problem was encountered in the gas heat exchanger which 
was used to control the shroud temperature. This problem was 
corrected by increasing the vent area of the unit eliminating 
the choked flow conditions in the vent* Two sets of transient 
runs were conducted as a result of the gas heat exchanger 
problem, however only the second set of runs was reduced and 
presented. 

Figures B-69 through B-78 present the temperature histories of 
the data that resulted from the 200 and 400 W tests. Figure B-78 
presents the average shroud temperature. The unconnected points 
indicate the desired temperature curve. The comparison between 
the desired curve and the actual curve was within a few degrees 
Kelvin throughout the transient test series. 
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Figure B^70^ Teat Z, Transient Environment Header Heat Pipe 
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Figure B^73 Teet 3, Transient Environment Radiator Heat Pipe Ho. 3 






























Ftgure B~7? 














Temperature, *K 








144'^K SHROUD TESTING 


The 144°K shroud testing was conducted to obtain additional cold 
case performance data instead of the testing that could not be 
conducted under the original plan. This test exceeded the cold 
case design environment for the radiator and was compensated for 
by increasing the load to be dissipated. 

Three power levels 400, 500, and 600 W were tested, however the 
lower load was terminated before achieving steady-state. The 
400 W load test was terminated early as the reservoir temperature 
continued to drop below the freezing point of ammonia, the feeder 
pipe working fluid. The power level was increased to 500 W and 
the test continued until steady-state was achieved. 

This test demonstrated the capability of the radiator to operate 
in a significantly colder environment (144°K versus 194°K) . The 
response of the radiator to changing load for this environment 
was also obtained. Figures B-79 through B-84 show the tempera- 
tures obtained for this series of testing by location on the test 
hardware. Figures B-85 through B-94 show the temperature his- 
tories of the data. 












Figul’e B-81 



Legend ; 

• Thermocouple Locations 


Figure B-^81 Test 4 144'^K Shroud^ 400 W 



Note : Temperature , 












Figure B-83 
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Figure 8^85 Teat 4, Thermal Conditioning Panel 


























KRPRODUCffilLITY 
OltiGlNAL PAGE II 




























REPRODUCIBn.nT 
ORIGINAL PAGE I! 















VII. 


CONCLUSIONS 


The heat pipe radiator/ thermal conditioning panel demonstrated 
that an all heat pipe thermal control system can satisfy future 
spacecraft thermal control needs. Wiile all of the heat pipes 
performed as expected, one conclusion drawn from the test was 
significant. 

The use of heat pipes in a thermal control system should minimize 
the number of physical interfaces between pipes and other equip- 
ment. Each interface introduces a AT into the system which will 
result in a lower system performance per unit of weight. This 
test was intended to integrate two existing heat pipe systems and 
demonstrate their operational capability. Three interface AT^s 
resulted in connecting the tx^o major test articles. Specific 
thermal control applications should consider making the header 
pipe or its equivalent an integral part of the collector device. 
One or two of the system interfaces can be avoided improving 
overall system performance. 
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ISOTHERMICS INCORPORATED • P.O. BOX 86* AUOUSTA* NEW JERSEY*07822 (201) 383*3300 


CS>ISBtherniics,inB. 

March 4, 1975 


I 


Martin Marietta Aerospace 
Denver Division 
P.O. Box 2371 
Denver, Colorado 80201 


Attention : Terry Ward 


Dear Terry: 

Enclosed please find a copy of our test report and lab 
data sheet concerning the heat pipe we manufactured for you. 

We did not know exactly how the pipe would be inter- 
faced with heat source and sink. We assumed that the heat 
input would be uniform over the 30 inch evaporator and that 
the heat out put would be uniform over the 48 inch condenser. 
Tests were run level and horizontal. It was decided to verify 
the upper load requirement and assume that the lower loads 
would be acceptable. 

If you have any questions regarding the report or data, 
please feel free to call me at any time. 


Very truly yours, 
ISOTHERMICS, INC. 

Robin Rhodes 
Project Engineer 

RR/dss 
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ATTACHMENT A 


March 3, 1975 


TO: R. Rhodes 

FROM: C. W3mian 

SUBJECT: Testing of Martin Marietta Pipe MO 5114 


1. Test Set Up 

A. Thermocouples No, 1 to 8 ere mounted on the evaporator 
plane in the first and third; bolt holes from each end. 

Thermocouples No. 9 and 10 are mounted on the heat pipe 
in the evaporator section. 

Thermocouples No. 11 and 12 are mounted on the heat pipe 
on the 90° bend. 

Thermocouples No. 13 to 16 are mounted on the condenser 
plate. 

B. Heaters 

Two heater wires were attached to the evaporator plate 
outside the bolt holes. Each heater measured 20 ohms 
resistance. 

C . Insulation and Cooling 

Evaporator section was covered with two layers of 2 in. 
fiberglass. The condenser section was cooled by evapora 
tion of a refrigerant. 

Tests 1 and 2, R-113 refrigerant used. 

Tests 3 to 10, R-11 refrigerant used. 

D. Pipe Set Up Level and Horizontal 

E. Load Evaporator - each heater had 20 ohms resistance. 
Voltage varied current measured. 

Watts calculated using I^R. 
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Results 

_ . , Watts Watts 

Over all pipe conductance T evaporator - T condenser 


0.0 44.4 

8.4 36.1 

34.96 
36.76 
44.6 


T Condenser 

I 

130 . 1 

3 

91.4 

3.6 

89.1 

4.0 

93.8 

5.0 

94.3 

5.5 


Note: Temperatures in °F 


10 Tests at five sets of conditions 

One and tx^o used R-113 cooling - insufficisnit cooling capacity 

Three changed to R-11 for cooling - uneven temperature on condenser 
indicated that R-113 was insulating portions of the condenser. 

Ran approximately hr to stability. R-11 attached Silicone in- 
sulation on the t/c's (condenser). 

Four and five R-11 uniform condenser temperatures restored. 

Six, seven and eight R-11 cooling excellent operation at 100% 
pox^er requirement. 

Nine and ten R-11 cooling - 210 W beyond requirement, tests 
terminated. 

Note: Due to time factor and lack of equipment, the loxxr limit 

(~30°F) could not be verified. 

Conclusions : 

1. Pipe is capable of meeting high end requirements . 

2. Overall pipe conductance based on the five conditions tested 
is 38,9 W per °F rise. 
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INTRODUCTION 


The analysis of the heat pipe radiator/ thermal conditioning panel 
test results was performed to demonstrate the ability of predict- 
ing multi-heat pipe system thermal performance. The test data 
from Appendix B was compared with the model results for both steady- 
state and transient conditions. To accomplish this comparison a 
112 node thermal model of the header pipe and radiator was de- 
veloped. Steady-state and transient runs were made and the results 
compared \‘7±th the test data. 

The predicted results agreed with the test results both in heat 
rejection and temperature. The most difficult task in the analysis 
was matching the front location in the individual feeder pipes on 
the radiator. 

The complexity of this analysis is significant when compared with 
the analysis of previous radiator systems using a pumped fluid. 

The number of variables introduced by the variable conductance 
heat pipes over-shadows those of the pumped fluid systems. While 
the pumped fluid radiators have been designed and analyzed for 
over ten years, heat pipe radiators are relatively new and their 
analytical techniques have not achieved the maturity of other 
thermal control systems and techniques. 

Over the past two years, Martin Marietta has devoted an IR&D task 
to the development of these analytical techniques and has verified 
the techniques for several fixed and variable conductance heat 
pipe designs including one radiator system. The analysis described 
in this report was the first data correlation attempted for an 
axial grooved variable conductance heat pipe. 
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II . SUMMARY 


The heat pipe radiator thermal model correlations with the test 
data demonstrated that heat pipe systems can be modelled as a 
part of an overall system thermal model. The predicted results 
agreed very well with the test data for all steady-state and 
transient runs . 

The added complexity of modelling heat pipes and including them 
in a thermal model is significant. The predictability of the 
front location in a variable conductance pipe represents the most 
difficult task. In general the predicted front locations appeared 
to be within 10 cm (approximately 5% of total length) of the 
actual location on the average with the maximum error being 30 
cm (approximately 15% of total length) . 

Future applications of these techniques and improvements in heat 
pipe modelling will provide even better results as heat pipe 
applications to spacecraft increase. The results indicate that 
heat pipe analytical techniques are now achieving the maturity of 
techniques developed for analyzing other thermal control devices 
such as pumped fluid systems. 
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III. MODEL DESCRIPTION 


The thermal analysis of the radiator/TCP test was accomplished 
using a 112 node thermal model of the system. The radiator con- 
tained 60 nodes, the header heat pipe 16 nodes, and the shroud 36 
nodes • 

The radiator model was divided into 10 nodes per feeder heat pipe. 
Each pipe contained one evaporator node, one transport node, six 
condenser nodes, one reservoir node and one vapor node. The 
header pipe consisted of four evaporator plate nodes j four evapo- 
rator nodes, one transport node, six condenser nodes, and one 
vapor node. Figure C-1 presents the nodal breakdown of the radi- 
ator and header pipes . Each of the heat pipes were modelled in 
the same general manner as shown in Figure C-2 for feeder pipe 1 
on the radiator. 

The evaporator of each feeder pipe was connected with the appro- 
priate header condenser node via a contact conductance. The con- 
tact conductance was calculated based upon the test AT and the 
predicted heat load per pipe Five steady-state runs (three hot 
case and two cold case) were used to determine the values shown 
in Table C-1. 

The average values in Table C-1 were used in the model to obtain 
the predicted performance presented later in the test. 

Heat x>ras applied to the header evaporator plate uniformly over 
the four-plate nodes based upon the measured heat loads from each 
test. An attemp was made to model the thermal conditioning panel 
(TCP) as a single node and apply the heat load at this point. This 
did not yield reasonable results due to a variation in the apparent 
contact conductance between the plate TCP. This variation 
is attributed primarily to variation in ^jemperature on the TCP 
under the plate. Since this temperature distribution was not 
measured in the test, the proper conductance value was unobtain- 
able from the model or test data. 

Since this was the first attempt to checkout the axial groove pipe 
simulation, problems with the axial groove logic in the subroutines 
were anticipated. Initially low ATs were obtained for the feeder 
pipes. A match of feeder pipe ATs was obtained by reducing the 
tube to vapor conductances by a factor of 18.6. The original con- 
ductance calculation did not account for conducting through a liquid 
layer in the area of vaporization or condensation. A brief liter- 
ature search revealed that a liquid layer of approximately 0.00254 
cm (0.001 in.) yields the appropriate conductance in series with 
the radial conductance of the tube (ref 1). Furthermore, Reference 
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1 states that the vaporization in the evaporator occurs in the 
vicinity of the meniscus attachment line while condensation in 
the condenser occurs primarily on the fin ends. For the results 
discussed later, a constant radial conductance for both the evap- 
orator and condenser sections was used. 
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Figure C-2 Heat Pipe Model 


Table C-1 Contact Conductances 


‘g ^ 
G 

Pipe 1 

Pipe 2 

Pipe 3 

Pipe 4 

Pipe 5 

Pipe 6 

10.26 

10.73 

9.38 

7.07 

8.69 

7.29 

Standard ^ 

Deviation ®K 

1.69 

1.90 

1.99 

1.65 

0.79 

1.97 


The early runs further indicated that using the gas load of 0.0135 
Ib^ of nitrogen (Ref 2) did not yield the proper gas front loca- 
tion. The model assumes flat front theory with tube conduction. 
This was not sufficient in locating the front until the gas load 
was reduced by 6%. The reduced load was within the loading ac- 
curacy accordin^to Grumman Aerospace Corporation (±10%) . The 
actual gas weighW is very difficult to weigh and the load was 
determined by the perfect gas equation. Some diffusion effects 
are also likely; however, this effect was considered to be small. 
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IV 


MODEL CORRELATIONS 


Model correlations with the test data yielded good agreement for 
the steady-state and transient runs both in temperature and front 
location. The discussion of the model results are presented first 
for the steady-state cases followed by the transient runs. 


A. STEADY-STATE CORRELATIONS 


1. Hot Case 

All of the hot case test runs were simulated by the model and 
the results are presented in Figures G-3 through C-12. Each of 
the six feeder heat pipes’ temperatures are plotted versus node 
location superimposing the test temperatures for direct comparison 
with the model results. The 100 W hot case is presented in Fig- 
ures C-3 and C-4. The front locations for pipes 3 and 4 were in 
excellent agreement with test while pipes 5 and 6 predicted fronts 
10 and 15 cm, respectively beyond the locations indicated by test. 
Pipes 1 and 2 predicted fronts 25 and 15 cm inside those in- 
dicated by test. This run yielded the most significant front 
location errors of the hot case analyses. 

In reviewing all of the data one should bear in mind that the 
reservoir temperatures were set to the measured test temperature 
with one exception. The thermocouple on reservoir 1 was loose 
and did not yield accurate test data. Its reading was adjusted 
to a more reasonable value. 

The 200 W case is presented in Figures C-5 and C-6 . The tempera- 
ture predictions were in excellent agreement with test as were 
the front locations. The 300 W case is shown in Figures C-7 and 
C-8. Pipes 1, 4 and 6 were in excellent agreement with test x<rhile 
pipes 2, 3, and 5 showed errors in front local^on 20, 12, and 17 
cm, respectively, beyond that indicated by test. The 380 W case 
is shown in Figures C-9 and G-10. This case resulted in excellent 
agreement for the first 5 pipes with pipe 6 predicting the front 
inside of that indicated by test by 32 cm. The 200 W test was 
repeated to show repeatability of results which also showed re- 
peatability of the model. The data is presented in Figures C-11 
and C-12. 
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Cold Case 


The cold gas correlations are summarized In Figures C-13 thru C-18. 
Figures C-13 and C-14 present the 200 W case. The model predicted 
temperature and front locations very well for this case. Pipe 5 
presented approximately 10 cm error In front location. The 400 W 
case. Figures C-15 and C-16, resulted In slightly larger front 
location errors, 17 cm for pipe 2, 10 cm for pipe 3, 5 cm for 
pipe 5, and 15 cm for pipe 6. The 500 W case. Figures C-17 and 
C-18, gave very good results for the first 5 pipes with the sixth 
pipe showing a significant error In front location. 

3 . 144"K Shroud 


The 144"K shroud 500 and 600 W cases showed nearly excellent agree- 
ment with test. The 500 W case Is shown In Figures C-19 and C-20. 
Temperature agreement was excellent for the first 5 pipes. The 
sixth pipe predicted the front 22 cm Inside that Indicated by 
test. Figures C-21 and C-22 present the 600 W data and show ex- 
cellent agreement In temperature and front location for the first 
5 pipes. The sixth pipe appeared to have at least a 30 cm error 
In front location. 

4. Discussion 

As shown In the above data general agreement In predicted versus 
measured temperatures was achieved by the model. Due to the lack 
of sufficient number of thermocouples from the test the precise 
location of the vapor gas front can not be determined. Assuming 
that the temperature profile Indicated by the model Is reasonable 
the apparent errors In front location were determined by shifting 
the model prediction to agree with test. Feeder pipe number 6 
gave the most Inconsistent set of predicted results. A further 
reduction In gas load below the -6% load used In the model would 
have resulted In poorer agreement. In all of the runs the energy 
balance of the model at steady-state was less than 1 W maximum 


Simulation of a variable conductance heat pipe Is strongly dependent 
on reservoir vapor pressure, which In the case of ammonia. Is an 
exceptionally strong function of reservoir temperature. The com- 
plex geometry In the vicinity of the reservoir made an accurate 
estimate of the radiation couplings to the shroud difficult. Con- 
sequently, the reservoir temperatures were Imposed In this analy- | 

sis In order to minimize these unknowns. All the steady-state 
runs contained logic to compute a radiation coupling from reservoir 
to shroud based on the heat leaking Into the reservoir from the 
condenser. Several cases resulted in apparent front movement Into 
the reservoirs, a condition which the model is not capable of 
simulating. When the average radiation couplings for the runs with 
the front not Influencing the reservoir were used between the res- 
ervoirs and shroud, predicted temperatures within 2^K resulted for 
the 500 W 144"K shroud case. 
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Figure C^IZ Cold Case Stoady-Stata, 200.39 W 
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TRANSIENT CORRELATIONS 


The transient data is presented in Figures C-23 through C-29 . 

In general the model predicted transient responses in very good 
agreement with the test data . The reservoir temperatures were 
controlled to the measured test values except for pipe 1 which 
was controlled to the same curve as reservoir 2. All of the 
capacitance of the system was limited to that of the four header 
pipe evaporator nodes Including the thermal conditioning panel . 

All of the remaining nodes were input as arlthemetic nodes. The 
reason for this was that the resultant time step would be extremely 
small indicating that the chosen method would yield a close 
approximation . Time step problems were encountered resulting in 
a sluggish response of the problem at an interval of 0.05 hr. 
Reduction of that interval to 0 . 02 hr yielded the presented data. 
The data presented here were derived from Appendix B superimposing 
the predicted results for comparison. 

In the header pipe comparison the upper set of model points rep- 
resents the average of thermocouples HI, H3, H5 and H7. The next 
cooler set represents H2, H4, H6 and H8. The lowest two sets 
are the hottest and coldest condenser results representing thermo- 
couples HIO through H15 . Figure C-23 presents the data comparison 
for the header pipe which resulted in very good comparison in 
temperature . The most severe transient conditions were chosen 
for this correlation representing the end of the 200 W run and 
the beginning of the 400 W run. 

In the feeder pipe comparisons the upper set of model points 
represents the evaporator, the middle set the first condenser 
thermocouple from the evaporator. Feeder pipe 1, Figure C-24, 
yielded excellent agreement with the test. Pipe 2, Figure C~25, 
gave results approximately 2°K warmer on the evaporator and dis- 
played faster response of the front movement compared with test. 
Pipe 3, Figure C-26, compared very well with the test and showed 
slightly faster response in front movement than the test. Pipe 
4, Figure C-27, showed good agreement in the evaporator and active 
condenser section. The agreement with the second condenser thermo- 
couple was mixed. The front did not respond out far enough during 
the Initial front movement resulting from the power change and 
later showed the front moving back into the pipe beyond that in- 
dicated by test. Pipe 5 , Figure C-28, yielded very good agree- 
ment except for the front tending to move back into the pipe be- 
yond that indicated by itest during the cold part of the orbit while 
at 400 W. Pipe 6, Figure C-29, gave good agreement in the active 
sections of the pipe but also showed a front location closer to 
the evaporator than indicated in the test result during the 400 
W portion. 
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Figure C-23 








Figure C-2S Radiator Heat Pipe Ho. 2 





















REPRODUCIBILITY OF TH 
c-34 ORIGINAL PAGE IS POOR 



Figure C-28 Radiator Heat Pipe No. S 













Discussion 

The transient results were in good agreement : with the test data 
except for pipes 4 and 6. A slightly sluggl^ih response to the 
load change in the model caused the front in pipe 4 to be too far 
toward the evaporator when the shroud temperature began to drop. 
Consequently, a much greater temperature excursion resulted in 
the model than in test. Preliminary runs in the 540 to 630-min 
region indicated that this excursion would be eliminated by that 
time. Pipe 6 showed the most inconsistent comparison with test 
throughout the analysis, and in the transient again a slightly 
further decrease in noncondensible mass would have resulted in 
better results. The error indicated in the evaporator of pipe 2 * 
is uniform throughout the analysis and indicates an instrumenta- 
tion problem especially since in some cases this thermocouple . 
registered temperatures below the condenser temperatures of the 
same pipe. 


Although the capacitance of the radiator was ignored in the tran- 
sient analysis, the response rate results indicated that the ther- 
mal conditioning panel was the dominant mass in the system. 




V. CONCLUSIONS 


The model correlations with the test data indicate the techniques 
of analyzing heat pipe thermal control systems are maturing and 
can be incorporated into the design analyses of spacecraft. Since 
future spacecraft will very likely depend upon heat pipes for 
controlling the vehicle thermal balance and equipment temperatures, 
this analytical tool will become more useful and important to the 
thermal design. 


Specific pipe configurations and designs still require empirical 
data for the model such as the determination of the characteristics 
of liquid layer in the axial groove pipe through which the heat 
is transferred. This can be evaluated through test and analysis 


yielding an effective ~ for use in modelling.. 


Mechanical joints 


are inherent problems in heat pipe systems and particularly in a 
radiator system. Considerable data is available for design appli- 
cations, however the performance of a specific system is still 
dependent upon the values and the repeatability obtained. 


The techniques employed in this analysis did not incorporate the 
front movement into the transition tube and reservoir. Generalized 
approaches to handle this problem would aid ones ability to pre- 
dict the reservoir temperatures with more accuracy during full 
open pipe operations. 
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.002365, REPEAT, 

36, 

• u 3 6 5 , 

REPEAT, 

6, 

ENDS 

5 84, 

.CJ3333, REPEAT, 

35, 

•COOOl, 

REPEAT, 

6, 

ENDS 


20. , PHFEAT, 6, 
1.31, REPEAT, 6, 


592, 1.31, REPEAT, 6, 
OPERATlNr, CONDITIONS 

602, 32. 17*,, REPEAT, 

603, 90., REPEAT, 7, 
6 C A 


ENOS 


90 

“ 1 , f 


, REPEAT, 
REPEAT, 


9-J., E 
1.23, 

7, ENOS BOOT FORCE 
ENHS ANGLE 3ETHEEN 

“ ti., IF 

-1., IF 


ENOS 


;tting angle (degrees) 

SPECIFIC HEAT RATIO 

FORCE ACCELERATION (IBM-FT/LBF- 
BOOY FORCE AND PIPE AXIS 


:rs 


EVAP 

EVAP 


H U. rv M » JH- 

BODY FORCE AND PIPE A> 
HIGHER THAN CONDENSER 
BELOW CONOFNSER 


610, 

0., REPEAT, 42, ENOS ABSORBED FLUX 



611, 

• u0^4l67, 

.084167, 

. GO 2667, 

•002667, 

.004167, 

.G04167 


• C u 41 6^ , 

. Q C4167 , 

. 0G2667, 

.002667, 

.00 4167, 

.004167 


• U C *4 1 6 ? , 

• C Ch 167 f 

. G02667, 

.002667, 

. CQ4167, 

.004167 


• 00*1 167, 

.'uuh167, 

. 0 G 26 6 7, 

.002667, 

.0X4167, 

.004167 


• 4.J4167 , 

• u 04167 1 

. w 0 26 6 7 , 

. JG2667, 

.004167, 

.004167 


.JG4167, 

.004167, 

. 002667, 

.002667, 

.004167, 

.004167 


612, 

613, 


j., REPEAT^ 6, ENDS FINTHIGKN 
.333333, REPEAT, 36, 0., REPEAT, 

7f^.- ^ DCDCAT-'?t;_ 1.. D CD TAT- C.. 


6, 

6, END: 


614, 

.85 . 


REPEAT, 36, u., REPEAT, 6 

, ENOS 

615, 

-26. 


-25. , - 2 4 • , -26., 

-28. , 

-31. 


-26. 


-25., -24., -26., 

-28., 

-31. 


-26. 


-25., -25., -27., 

-28., 

-31. 


-27. 


-25.J -25«, ^ 27 » f 

-29., 

-31. 


-27 . 


-26. , -26. , -28 ., 

-29., 

-31. 


-27. 


-25., —26., -28., 

-29. , 

-31. 


u. , 

REPEAT, 6, END 



6 16, 

169 . 

♦ 

REPEAT, 36, 0,, REPEAT, 

6, ENDS 


WICTH (FT) 

K (BTU/FT-HR-F) 
EMTSSIVITY 


617, 

618, 


6, ENDS 


C-...ARTHRY PARA»^ETERS 
G. ...EVAPORATOR 


TIE TO ADJACENT 
/ TTcn rn An.iArcK 
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c 

c 

c 

c 


c 


c 

c 

c • • * • 


C • • t • 


r 

i 


633, a 


621 , a. 

622, 0. 
623 , 0 * 
6 24, u • 

625, C. 

626, 3. 
6 27, a. 
6 2 3 , 3 « 
629-1 J. 
6 3 u , u « 
6 31, 0. 

632, 0. 


.REPEAT 


>R6P£AT 

, REPEAT 
, REPEAT 
, REPEAT 

, REPEAT 
, REPEAT 
, REPEAT 
, REPEAT 
^REPEAT 
, REPEAT 
, REPEAT 

, REPEAT 


,13, ENDS TYPE OF ARTERY (INTEGER) — 
u = NONE 

1 = PEDcSTAl 

2 - SPIRAL WITH SPACER WIRE‘S 

,13, ENDS OUTER OIAM (FT) — USE EFF OIAH IF 
NOT CIRCULAR? EFF D = (4.^ X-S AOEA / PT)'^».5 
10, ENOS INNER DIAM <FT) — USE EFF Q 
10, ENOS OIAH OF SPACER WIRES IF M620+I = 2 (FT) 
iJ,FND$ ARTERY SPIRAL SCREEN THICKNESS (FT) 

IF M620+I = 2 (FT) 

iJjENnS ARTERY SCREEN PORCSTTV (01 ME NSIC NL E SS) 
13, ENOS NC CP SPIRAL TURKS IF M620+T=2 (REAL) 
10,ENOSART£RY SCREEN OENSITY^GP FROOUCT 
IJrENOS SPACER WIRES DENSITY *CP PRODUCT 
13 ,ENDSARTERY SCREEN THPML CCNOUCTIVTTY 
1 J ,£Nn$SFACER WIRE THERMAL CCNCUCTiVlTY 
1 j,EM 03 THICKNESS C^ OUTER ARTERY W^AP 
M62wH = 2 (FT) 

la.ENOS EFF RADIUS OF LARGEST ARTERY HOLE 
EXCLUCING THE CCNOEMSER END OF THE ARTERY (FT) 


(aTU/FT»»3-F) 
It=' M62C + I = 2 
(^TU/HR-FT-F) 
TF M62C+I = 2 
IF 


ADIABATIC SECTION 


6 4 w , j 

, REPEAT 

,7, £NO 

641v 0. 

.REPEAT 

.7. END 

6 4 2, u • 

, REPEAT 

,7, END 

643, 0. 

, REPEAT 

,7, END 

6h4, 0. 

.REPEAT 

,7, END 

645, 0. 

.REPEAT 

,7, END 

6 H 6 , 3 • 

.RcPt AT 

,7. END 

£^7, 3. 

, REPEAT 

,7, END 

6 4 6, w • 

.REPEAT 

,7, END 

649, 3. 

.REPEAT 

,7. END 

650 . 3. 

.REPEAT 

.7.ENO 

6 51, 0* 
ONDENSER 

.REPEAT 

,7, END 

6 6u , u 

, REPEAT 

,42, END 

661. 0. 

.REPEAT 

.42, END 

6 6 2, Q • 

.REPEAT 

,42 ,ENO 

663, 

.REPEAT 

,42, END 

664, 3. 

, REPEAT 

,42, END 

6 6 3, 0 4 

.REPEAT 

, 4 2 , E NO 

6 6 6. 0 • 

.REPEAT 

,42. END 

667, u • 

.REPEAT 

,42, END 

6 6 6, 0. 

.REPEAT 

,42, END 

6 6 9, u • 

.REPEAT 

,42 ,END 

6 7 u . 0 • 

.REPEAT 

.42. END 

671, U4 

.REPEAT 

,-.2, END 


C *** RESEPVOl.^ PARAMETERS 


8 3 6. 

.01269/28. 016, 

REPEAT, 

6, 

0., END' 

MOLES NCG (L"MOLES) 

8 0 9. 

33.4/1728., 

REPEAT, 

6» 

C. , EN05 

RESERVOIR VOLUME (FT*»3) 

822, 

6. /1 2., 

REPEAT, 

e, 

0 ENOS 

RESERVOIR LENGTH (FT) 

823, 

2.902/12., 

REPEAT, 

6, 

0 * , K 0 3 

RESERVOIR WICK 00 AND RESERVOIR IC 

324, 

2.fi«5/12. , 

REPEAT, 

e, 

C . , ENDS 

PfSERVCIR WTCK ID (FTV 

82 5» 

3* ju/12* . 

REPEAT , 

6f 

0 . . ENDS 

RESERVOIR 00 (FT) 

6 2 6, 

2* 857 12# . 

REPEAT, 

6, 

0 . , END3 

FFEOTUEE LENGTH (FT) 

82 7. 

.381/12.. 

REPEAT, 

6, 

0. .ENDS 

FEEOTUEE WICK 00 ANO FEEOTORE 10 

8 23, 

.380/12., 

REPEAT, 

6, 

0 ..ENDS 

FEEnTUEE MICK 10 (FT) 

829. 

.437/12., 

REPEAT, 

6, 

0 . , ENOI 

FEEOTUEE 0 0 (FT) 

8 3a, 

8., 

PEFEAT, 

6, 

0.,END3 

ffedtuee thermal CONO (BTU/HR-FT-F 

8 31. 

0. . 

REPEAT, 

6, 

0 . , E NO 3 

IN'SULATICN CCNO (RTU/HR-FT-F) 

837, 

8 • , 

REPEAT, 

6, 

u . .ENDS 

RESERVOIP THERMAL CCND (SAME) 

638. 

0 . . 

REPEAT, 

6 , 

C. ,EMD3 

INSULATICN THICKNESS (FT) 


C **• ANGULAR WIDTHS OF EVAPORATOR SPIRAL ANNULI (DEGREES) 

910, J, REPEAT, 7, ENDS NUMBER OF ANNULI IN EACH PTPE 

911, 0., REPEAT, 7, ENDi ANGULAR WIDTHS OF ALL ANNULI 

C *** ANGULAR WIDTHS CF ABIARATIC SECTION SPIRAL ARTERY ANNULI (DEGREES) 

912, J, REPEAT. 7, ENOi 
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f,'RP10BlJC]BII-rrf OF TIO* 
PAGK IS mmi 


913. a. .REPEAT, 7, END* 

C ••• ANGULAR WIDTHS CF CONDENSER SPIRAL ARTERY ANNULI (DEGREES) 

914. 0, RER-flT, 7, ENOS 

915. j., REPEAT, 7, EN'Oi 

C... ARTERY CAPACITANCE (9TU/F) (COHFUTEn IN SU3 ARTFRY) 

681,SPAC£,i:,EN0iEVAF 
682,SPACE,7,ENDS ADIA 
683, SPACE, 42, END? CONB 

C... ARTERY AXIAL CCNCUCTANCE (9TU/HR-F) (COHPUTEO IN SUB ARTERY) 
6d4,l.E-12,REPEAT,lC,END3 EVAP 

685. 1. E-12, REPEAT, 7 , ENOS AOTA 

686.1. £-12, REPEAT, 42, END* CONO 

C... EFFECTIVE VAPOR PASSAGE DIAMETER (FT) (COMPUTED IN SUB ARTERY) 

687, SPACE, 10, ENOS EVAF 
688, SPACE, 7, ENOS AOIA 
689, SPACE, 42, ENOS CONO 

C...(PERMEABILITY)»(X-SECTI0NAL AREA) PRODUCT (FT**4) (COMPUTED IN SUP PERMnD 
690, SPACE, i:, ENDS EVAF 
691, SPACE, 7, ENOS AOIA 
692,SPAC£,42,EN0$ CONO 

C...HALL AXIAL CONDUCTANCE (BTU/HR-F) (COMPUTED IN SU« AXCOND) 

693. 1. E-12, REPEAT, lo. ENOS EVAF 

694. 1. E-12, REPEAT, 7, ENOS AOIA 

695. 1. E-12, REPEAT, 42, ENOS CONO 

C...WICK AXIAL CCNOUCTANCE (BTU/HR-F) (COMPUTED IN SUB AXCOND) 

696. 1. E-12, REPEAT, 10, ENOS EVAP 

697. 1. E-12, REPEAT, 7, ENOS AOIA 

698. 1. E-12. REPEAT, h2, ENOS CONO 

C...EVAP WORKING FLUID LIQUID PROPERTIES (OETHPMINEP IN SUa LOCKUP) 

701, SPACE, i;, ENOS LIOUIO DEMSTTY {LPM/FT**3) 

732, SPACE, 10, ENOS LIOUIO SPECIFIC HEAT (BTU-'LRH-F) 

70 3, SPACE, 13, ENOS LIOUIO VISCOSITY (LBM/FT-HR) 

704, SPACE, 10, ENOS LIQUID THERMAL CONDUCTIVITY (gru/HR-FT-F) 

7C5, SPACE, 10, ENOS SUPFACE TENSION (LSF/FT) 

7Q6,RPACE,10, ENDS HEAT OF VAPORIZATION (pTU/IBM) 

73 7, SPACE, 1. , ENOS HEAT OF FUSION (OTU/LpM) 

C...AOIA working FLUID LIOUIO PROPERTIES (OETEP'^IMEO IN SUB LOOKUP) 

711, SPACE, 7, ENDS LIQUID DENSITY (LeM/FT**3) 

712, SPACF,7, ENOS LIOUIO SPECIFIC HEAT (RTU/L°M-F) 

713, SPACE, 7, ENOS LIOUIO VISCOSITY (LPM/FT-HR) 

714, SPACE, 7, ENOS LIQUID THERMAL CONDUCTIVITY (P TU /HR-FT -F) 

71S, SPACE,?, ENOS SURFACE TENSION (LBF/FT) 

716, SPACE, 7, ENOS HEAT CF VAPCRIZATICN (BTUXLPM) 

717, SPACE, 7, ENDS HEAT OF FUSION (BTU/LBM) 

C...CONO WORKING FLUID LIQUID PROPERTIES (DETERMINEO IN SUP LOOKUP) 

721, '•PACE, 42, end? LIOUIO DENSITY ( UPP/FTT’S) 

^22,SPACE,42,£N0i LIOUIO SPECIFIC HEAT (PTU/Lpm-F) 

723, SPACE, 42, ENO:- LIOUIO VISCOSITY (LBM/FT-HR) 

724, SPACE ,42, ENOS LIQUID THERMAL CCNOUCTIVITY (BTU/HR-fT-F) 

725, SPACE, 42, ENCE SURFACE TENSION (LRF/FT) 

726, SPACd, 42, ENOS HEAT OF VAPORIZATION (PTU/LPM) 

727, SPACE, 42, ENC? HEAT OF FUSION <aTU/LPM) 

C... WORKING FLUID VAPOR PROPERTIES (DETERMINED IN SUE ClOEGl) 

731, SPACE ,7,ENnS VAPCR DENSITY (LBK/FT»»T) 

732, SPACE, 7,EN0S VAPOR SPECIFIC HEAT (BTU/LBM-F) 

733, SPACE, 7, ENOS VAPOR VISCOSITY UBM/FT-HP) 

734,SPACE,7,ENDI VAPOR THERMAL CCNOUCTIV TTY ( oju/HR-FT-F) 

C... OPERATING CCNOITIONS (CALCULATEC IN SUP LIMITS) 

751, SPACE, -,£NCS MA'^S FLOWRATE (IBM/HR) 

752*SPACF,7,ENni LATENT HEAT TPANSFFR (PTU/HR) 

753. SPACE, 7, ENDS HICKING LIMIT (PTU/HR) 

754, SPACE, 7, ENOS SONIC LIMIT (BTU/HR) 

755, SPACE, 7, ENTS ENTRAINMENT LIMIT (BTU/HR) 

756,SpACf,7,.FN0S BOILING LIMIT (BTU/HR) 

757, Space, 7, ENOS ovefall limit (lowest of hzsj thru a25b) (pru/HR) 

758, SPACE, 13, ENOS EVAPORATOR VAPOR PRESSURE (PSIA) 
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759,SPACEf7,EN0J A0IA9A7IC SECTION VAPOP PRESSURE (PSIA) 

760, SPACE, 42, ENOS CONDENSER VAPOR PRESSURE (PSIA) 

761, SPACE, 7, ENOS RESERVOIR VAPOR PRESSURE (PSIA) 

763,SPACE,7,ENOi INJECTION REYNOLDS NUMPER ( OINENSIONLESS) 
76',,SPAC£,7,£N0j SUCTICN REYNCLOS NUMDER (DIMENSIONLESS) 
C...LIQUin/HICK MATRIX EFFECTIVE RACIAL THERMAL CONDUCTIVITY ( pTU/HR -FT -F ) 
C.... (COMPUTED IN SUB EFFKMK) 

765, SPACE, 13, ENOS EVAP 
766,SPACE,7,EN03 AOIA 
767, SPACE, 42, ENOI CONO 

C, ..PANEL ENVIRONMENTAL PARAMETERS (COMPUTED IN SUB RADFIN) 

76 e, SPACE, 42, ENOS NODAL £N VI RONNEN T AL TEMPERATURE (F) 

770, SPACE, 42, £NOS NOOAL HEAT REJECTION (“TU/HR) 

771,SPACF,7,ENDS PIPE ENERGY BALANCE (DTU/HR) 

811, SPACE, 7, ENOS PIPE HEAT REJECTION (RTU/HR) 

C... FLAGS 

781, SPACE, 7, ENDS NC OF ITERATIONS IN SURRCUTINE FRONT 
783, SPACE, 7, ENOS BACKUP SWITCH IN SUBROUTINE FRONT 
C,,.VCHP PARAMETERS 

810. . 6. .6. .5. .5. .4. .4.1., ENDS ACTIWE7T0TAL CCNCENSER LENGTH (C-LESS) 

820. 3. 6.3. 6. 3. . 3.. 2, H, 2. 4, 4. , ENOS NEW ACTIVE CONDENSER LENGTH (FT» 
832, SPACE, 7, ENOS CL 0 ACTIVE CONDENSFR LENGTH (FT) 

836, SPACE, 7, ENOS NUMBER OF INACTIVE CONDENSER SECTIONS 
841, SPACE, 7, ENDS OSCILLATION TEST 
842,SPACE,7,EN01 OSCILLATION TEST 
843, SPACE,?, ENDS OSCILLATION TEST 

844.2., Rc.PFAT,7,EN03 INITIAL DAMPING 
C...FIN capacitancf and axial CCNOUCTAMCE 

853, SPACE, 42, ENOS CAPACITANCE (COMPUTED IN SUB CAPFTN) 

851, SPACE, 42, ENDS AXIAL CCNOLCTANCE (CCMOUTEO IN SUP RACFIN) 

C... WORKING FLUID THERMOPHYSICAL DATA 

IflSl? LIOUIP DENSITY (Lpm/FT*» 3) VS TEMP (F) - AMMCNTA 
-107.9,46., -g..,45.3, -50., 43. 5, C.,41.34, 50. ,39.1, IGO.,36.4 
15C.,33.6, 23o.,2S.5, 274.3,16., ENC 
13021 VAPOR DENSITY (LBM/FT*''3) VS TEMP (F) - AMMONIA 

-1C7.9, . CC57, -E4.,.3t,7, -50. ,. 0332, 3. ,. 1397, 50. ,.3036 

101.. ..7,48, 124, ,1.01, 150. ,1.54, 200 . , 3. 45 , 251 . , ° .1 

270.3.16., END 

10035 LIQUID SPECIFIC HEAT (PTU/LBK-F) VS T£MP (F) - AMMONIA 
■-137.P, 1.032, -75., 1.3h5, -58. ,1.066, -22. ,1.369, -4. ,1.077 

14.. 1,C9, 32.,l.lj7, 5i-., 1.126, 68. ,1.146, 86. , 1.168, 104. ,1.194 

122. . 1.222, 15u., 1.29, 175. ,1.365, 20-., 1.45, 225. ,1.555, 250. ,1.7 
270.3,1.86, END 

10043 VA“CR SPECIFIC HEAT (BTU/LBM-F) VS TEMF (F) - AMMONIA 
-107.9,0., ENDS UNNECESSARY 

10055 LICUIO VISCOSITY (LBM/FT-HR) VS TEMP (F) - AMMCNIA 
-1-7,9,1.415, -92.2,1.1484, -58., . 7668, -40 . ,. 666, 20. ,.3888 

43.. .3384, 63., .2916, ICO., .2197, 150. ..153. 200 ., .0^8 

250. . .0533, 27C, 3, .03384, END 

10365 VAPCR VISCOSITY (LBM/FT-HP) VS TEMP (F) - AMMONIA 
*107.9, .u1d2o, — 5-.,.C1863, 32. ,.02221, 68. ,.02376, 122., .02642 
153 . ,.32755, 20C., .33325, 270.3, . 03384, END 
10375 LIQUID THERMAL CONDUCTIVITY (BTU7HP-FT-F) VS TEMP (F) - AMMONIA 
*lL7.9,.32tJ, — 2-#,. ^17, 3<i.,.3l2, 53. ,.307, 8—., .293, 120.,. 275 
l5o».«254, 4u0.t«2-5, 225. ,.171. 250., .11, 270. 3, .068, END 
10Q8S VAPOR THERMAL CONDUCTIVITY (8TU/HR-FT-F) VS TEMP (F) - AMMONIA 
■* 0.1. 7, 9, .0 492 , -75. ,. 33955, -43. ,* 0105, C‘.,.0117, 40., .0129 
8 j«, .0142, ISO... -149. 120. ,*3156. 150*. .0173. 2C3...D223 

253.. , 0 36, 270. 3,. 368, END 

10395 SUPFAC? TFNSICN (LDF/FT) VS TEMF (P) - AMMCNIA 
-13 7.9, , 3-345, -6-. ,. 30278, -2-.,.uj281, 20., .00189, 52. ,,001576 

93. 4. . 03124, 1 36 . 2 . . 0 -- 891 . 150.. . 03-81. 17 5 . . * - CO 6 35 , 200., .00047 

225. . .03029, 250. ,.30313, 270.3,0., END 

1313? HPAT OF VAPORI7ATION (pTU/LBM) VS TEMP (F) - AMMONIA 
■“1-7*9,639., — 5'-. ,6-4, 3, -.,:>68.9, P-.,527. 3, US. ,477.8 

124. , 453.1, 153.. 4:g.. 175..365., 20 0 . . 328, , 2 25 , , 2 49 . 
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250.fl53.» 27C.3,J., tND 

10113 HEAT CF FUSICN (ETU/IBK) VS TEMP (F) - AMMONIA 
-lu7.9,j. » ENBS UNNECESSARY 

1J123 SATURATION PRESSURE (PSIA) VS TEMP (F) - AMMONIA 
-lu7»9f«85» -75.»3.6(, -50»»7.67j ~25«»15»5» fl«j3C»A2» 25r»53«2 

50. . 89. 1E, 75. ,143., 100. ,211.9, 12A. ,303.4, IEO.,425., 175. ,580. 

200. . 765., 225. ,1008., 250. ,1330., 270.3,1650., END 

1013$ THERMAL CCNOUCTIVITY < 0TU7HR-FT- F) VS TEMP (F) - 6061 ALUMINUM 
~26C.,65., ”160 . ,76., “60., 84., 4ii.,91., 140., 96., 240. ,101., END 
1021$ LIQUID DENSITY tLfiM/FT»»31 VS TEMP (FI - GENETPON 31/lli, 
-1C0.,9I..6, “5C., 90.6, 0., 86.4, 50. , 81. 9, IOC. ,77. 

15C.,71.5, 23J.,64.8, END 

1022S VAPOR DENSITY (LBH/FT**3» VS TEMP (FI - GENETRON 31/114 
-ICO. ,. 01224, -53. ,.0647, 0.,.2298, 55. ,.627, 100. ,1.4339 

150.. 2.93u6, 220., 5.6593, END 

1023$ LIQUID SPECIFIC HEAT (PTU/LSM-F) VS TEMP (F) - GENETRON 31/114 
“ICO. ,.21, -50. ,. 22, u»,. 237, 50. ,.256, ICO. ,.276, 150.,. 303 
20:.,.348, END 

1024$ VAPOR SPECIFIC HEAT (BTU/LBM-FI VS TEMP (FI - GENETRON 31/114 
-ICC.,0., ENOS UNNECESSARY 

1025$ LIQUID VISCOSITY (i.8M/FT-HR) VS TEMP (F) - GENETRON 31/114 
j.,. 9676, 25. ,. 8346, 50.,.7257, 75. ,.64, IOC., .5685, 125., .5 


i5C.,.H32,£ND 

1026$ VAPOR VISCCSTTY (LBM/FT-HR) VS TEMP (F) - GENETRON 31/114 
-ICC. ,.0195, 2C0.,.O345, END 

1027S LIQUID THERMAL CCNO (8TU/HR-FT-FI VS TEMP (FI - GENETRON 31/114 
— 1C>*.,.C755, 2^0. ,.036 9, END 

1026! VAPOR THERMAL CONQ (BTL/HR-RT-F) VS TEMP (FI - GENETRON 31/114 
-ICC., .00356, 202. ,.012 3, END 

1j29$ SURFACE TEN’SIQN (LBF/FT) VS TEMP (F) - GENETRON 31/114 
- 1C C .,. 0 03391 , 200 . , . 00 0415 , END 

1030$ heat OF VAPORIZATION (PTU/L9M) VS TEMF (FI - GENETRON 31/114 
-ICO., 112. 72, -50. ,108.77, Q.,104., 50. ,98. 06, 100. ,90. 66 


C...FIN 
600$ LTV 
13, 


151.. 81.22, 2C0.,68.'4E, END 

1031$ HEAT CF FUSION (BTU/LPH) VS TEMP (F) - GENETRON 31/114 
-ICG.,0., ENOS UNNECESSARY 

1032$ SATURATION PRESSURE (PSTA) VS TEMP (FI - GENETRON 31/114 
-100. ,.503, -50. ,3. 001, -20. ,7. 085, 3., 11.73, 20. ,18. 55 

43. . 28.18, 60. ,41. 35, 80. ,58.84, 10 0 . , 81 . 5 2 , 120 . , 110 . 3 

140. . 146. 2, 160. ,190. 2, 180 . , 243, 5, 2u3 . , 30 7.4 , END 
GRADIENT TA"LE 

DATA- Dt^'£NSIC^'LESS TEMP GRADIENT FOR FLAT EXTENDED SURFACES 


12, . 00, .6030, .70 30, .8003, .9300,1.003, 

1.10130 9 

1.2uOG ,1 

•3000, 

1.4000,1.5300 

. u C, . 2303, . 1503, .IGCC,. 0503,. 000 3, 

-#u5Q0, 

-.1000 

.1500, 

-.2000,- 

.2500 

. 13, .2566, .2133, .1708, .1293, .0890, 

.0499, 

.0120,- 

.0247, 

-.0597,- 

. 0939 

. 2 5.. 331 4, .2948,. 2598, .2265, .1949, 

.1649, 

.1336, 

.1100, 

.0 849, 

.0613 

. 7 i. , . 436 9, * 41»79, • 38c4, .355 2, . 3317 , 

.3099, 

.2898, 

.2713, 

.2542, 

. 2335 

. 7 5,, 5275, .5.- 32, .4 80 6, .45 99, .4411, 

.4239, 

.4383, 

.3940, 

.3810 , 

.36 92 

1. 00,. 607 4,. 58 63, .5671,. 5497, .5340, 

.5196, 

.5 070, 

.4954, 

.48 49, 

.4755 

1. 50,. 7458,. 7292,. 71 43, .7011, .6893, 

.6783. 

.6 694. 

.6610, 

.6535, 

• 6468 

2. 00, .8650, .8513, .8392, .3285, .8191, 

.8109, 

.8034 , 

.7969, 

.7911, 

.^860 

2. 5 J, .9709, .9593, .9491, .9402, .9324, 

.92 55, 

.9195, 

.9142, 

.9095 , 

. 90 54 

3. 00,1.067,1.357,1. 048, 1,341, 1,034, 

1.028, 

1. J23 , 

1,019, 

1.015, 

1.011 

3. 50,1.156,1.147,1.139,1.133,1.127, 

1.122. 

1.118, 

1.114, 

1. 110, 

1.107 

4. 00,1.239, 1.231,1.224, 1. 218,1.213, 

1.2u8, 

1.205, 

1.2 Cl, 

1.193, 

1.196 


10,12, .53, .6000,.7COO, .80 0 J,.9QQ0 
. Cu,.2U03,.15j0,.ljuut*u500 
.li.. 2508,, 20 7 5,. 1651,. 1237 
, 2 5,. 31 8:,. 2 81?-. 24 67,. 2 13 9 
. sc, .413 5,. 38 46, .3578, .3 ■'30 
.75,. 495 3, .4712, .4495, .4295 

1. 00. . 5670.. 5469. .52 85.. 5 120 

1. 50. . 6916. .6763. .6626. .6504 

2. 00. . 7993. .7868. .7759. .7664 

2. 50. . 8941. .8842. .87 54.. 8 67 7 

3. 03.. 9807. .9723. .9649. .9566 


1»OC3,1.100C,1.2COO 
.0000, -.0530, -.1000 


.3834, 
.1824, 
.3399, 
.4112, 
.4970, 
.6395, 
.7580, 
. 860 9, 
.9531, 


.3443, 
.1527, 
. 2835, 
.3946, 
.4834 , 
. 62 98 , 
.7505, 
.8550, 
.94 82, 


.0 065 
.1246 
.2687 
.3793 
.4712 
.6210 
.7439 
.8498 
.9439 


1.3000,1.4000,1. 5000 
-.1500, -.2000, -.2500 
•. 0301, -.0651,-. 0992 
.0981, .0732, . 0498 
.2536, .2338, .2184 

.3656, .3530, .3416 
.4601, .4501, .4411 
.6133, .6064, .6001 
.3381, .7328, .7282 
.8452, .8411, .8374 

.9402, .9369, .9339 
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lUfiZ 


C774, 

1691, 



10 , 12 , 


• JOf *2uuij| 

• Ut aZldOl 

.25,.245wj 

r» f ^ 


lUfic 


1.033, 1.029, 

1.C26, 

1.024» 

1.021 

1.111, 1.103, 

1.106, 

1.103, 

1.101 

.. 1000, 1.2u00 ,1 

.3000, J 

..4000,1 

.5000 

. 0500, -.1000,- 

.15QO>-t200iJ»- 

.2500 

.0384. .0006,- 

.0359, -.0709,- 

.10 49 

.1396, .1117, 

.0854, 

.0607, 

.0375 

.2657, .2463, 

•2205, 

.2121, 

.19 70 

.3635, .3488, 

.3355, 

.3232, 

.3122 

.4449, .4333, 

.4227, 

.4133, 

.40 46 

,5784, .5703, 

.5632, 

.5567, 

.55 09 

.6878, .6319, 

.6766, 

.6720, 

.6678 

.7821, .7776, 

.7736, 

.7700, 

.7669 

.8660. .8625, 

.8593, 

.8565. 

. 8540 

.9422, .9393, 

.9368, 

.9345, 

.9326 

1.312, 1. UlO , 

1.GC8, 

1.006, 

1.00 5 

..1000,1.2000,1 

..3wC0,1 

L.4000,1 

. 5000 

•.0500,-.1000,- 

’. It:lj0.*i2000.* 

-. 2500 

.Q286,-.C C9C,- 

•o 3454, -.3804,- 

••1142 

.1181, .0905, 

.0646, 

.0402, 

.0175 

.2284, .2096, 

.1925, 

.1765, 

.1619 

.3130, .2992, 

•2866. 

.2751, 

.2645 

.3829, .3721, 

.3624, 

.3536, 

.3455 

.4963, .4892, 

.4830, 

.4772, 

. 4721 

.5884, .5836, 

.5791, 

.5752, 

.5717 

.6676, .6640, 

.6608, 

.6579, 

. 6554 

.7378, .7350, 

.7327, 

.7305, 

. 7286 

.8313, .7993, 

.7974, 

.7958, 

• 7943 

. 85 99 , . 8 p 8 3 , 

.8568 , 

.8556, 

. P544 

L. 1020,1. 2003, i 

L.3C00,; 

1.4000«1.^000 

3500, -.1 COO,- 

-.15;q, 

-,2000,-. 2500 

.0175, -,0238, -.0588, 

-.0950,- 

-.1285 

•Q853, .0583, 

.0329, 

.0091,- 

-.0133 

.1721. .1543. 

.1380, 

.1228, 

.1089 

.2375, .2248, 

.2131, 

.2026, 

.1929 

.29J7f .2812, 

. 27 24, 

.2647, 

. 2575 

.3758, .370C, 

. 3648 , 

•3600, 

.3556 

.4444. .4406. 

.4371, 

.4341, 

.4313 

.5028, .5003, 

•4979, 

.4958, 

.494C 

.5546, .5528, 

.5511, 

.5497, 

.5484 

•6015, .6C02, 

.5q9Q, 

.5979, 

. 5970 

.6447, .6437, 

.6428, 

.6421, 

.6414 

1.1000,1.2050,1.3000, 

1*40 00 ,; 

1.5Q0Q 

*.0500 ,'*.1400 

-.15 00, 

-.2000, • 

-. 2500 

.0038, “.0336, 

-. 0697, 

-.1047,« 

-.13 79 

•0637. .0370, 

.0118. 

-.0117, 

-. 0336 

.1353, .11P2, 

.1024, 

.087^, 

. 0743 

.1886, .1765, 

.1655, 

.1555, 

• 1464 

.2314, .2225, 

.2145, 

.2072, 

. 2007 

.2992i .2941, 

.2894, 

.2852, 

.2814 

.3534, .3502, 

.3473, 

.3447, 

.3423 

.3996, .3974, 

.3955, 

.3936, 

• 3923 

.4403, .4389, 

.4576, 

.4365, 

.4354 

.4773. .4763. 

• 4754, 

.4746, 

. 4739 

.5113, .5106, 

.5099, 

.5094, 

. 50 89 

1.1000, 1.2000, 

1.3GC0, 

1.4000, 

1. 5000 

" • 0 5 0 0 , . i -C 0 u , 

-.1500, 

2U00, 

-, 2500 

. C081.-.0453, 

-.0812, 

-.1159, 

-.1492 

.0381, .0119, 

-.01 29 , 

-.0362, 

-.0577 

.0921, .C756, 

.C6C4, 

.0463, 

.2334 

.1312, .1200, 

.10 98i 

.10 03 , 

. 0918 

.1623. .1543. 

.1469» 

.1403, 

.1342 

.2108, .2063, 

.20 23, 

.1986, 

.1953 

.2491, .2464, 

.2440, 

.2418, 

.2398 

.2815, .2798, 

.2783, 

.2769, 

• 2757 

.3101. .3093, 

. 3C80, 

.3071, 

. 3063 

.3360, i3352, 

.3346, 

.3340, 

. 3334 
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E^RODUdigiLfH OF THE 
ORIGINAL PAGE IS 


4. aa 

iOflZ, . 95 

• u 0 

• u 

• 25 
.5'j 
.75 

1. 00 
i .50 

2. JO 
2. 50 

3. 00 

3.50 

4.00 
10tl2,i.C 

.00 
. 10 
.25 
.50 
. 75 

1.00 

1.50 
2.00 

2.5 0 
3.00 

3.5 0 
4. 00 

lOf 6 y 1. It 
. y 0 
.10 
‘ .25 

.50 
.75 
i. 00 
10y6yi.25 
.00 
.10 
.25 
.50 
.75 
l.OQ: 
10f4yl. 4C 
.00 
.10 
.25 
. 5 0 

10,3,1.50 
.00 
.10 
.25 

10,5,1 .75 

.00 
.10 

• 2 5 
.50 
.75 

10,3,2.00 
.DC 
.10 
.25 
2ND 
2QC0. 
ENC 

3 

MTEST 


. 3 63 4 
. 6000 
I 2CCC 
► iai3 
.1612 
1 141 9 
.1336 
.1314 
• 1359 
.1456 
. 1571 
. 1689 
.1806 
, 1919,. 

• 6 0 0 u , 
.2Q0C, 
.1645, 
.1235, 
.0787, 

• 0521, 
.0356, 
.0179, 
.Cu97, 
.0056, 
.0034, 
. 0021 , 
.0013, 

. 6 0 0 g , 

. 2 0 J G , 
.1213, 
.0295, 
-.0740, 
-.1413, 
-.1892, 

• 6 y 0 J , 
« 24 u u , 
•0322, 

-.16G6, 
-.371C, 
-.50 82, 
-.6107, 
.6CC0, 
.200 C, 
-.0941, 
-.4179, 
-.7517, 

. 6 w j c , 

. 2 0 0 C , 
-.2018, 
-.6292, 
•6000 , 
. 2 u 0 0 , 
57-^2, 
•1.295, 
•1.927, 
-2.384, 
. 6 0 u 0 , 
. 20 u D , 
- 1 . 1 . 2 , 
■2.147, 


3625 
70 00 
15C0 
1385 
1273 
1187 
1174 
1198 
1296 
142u 
1549 
1675 
17 97 
1914 


.36 l7 , .3610, .36u4, .3598, 
.8000 , .9000,1.000, 1.1000, 
• IGOu, .0500, . uuuO.'*.C50C. 
.0 967, .0 56 G, . 0 164, - . D219 , 
.0923, .y 644, .0354, .0082, 
.0973, ,0773, .0588, .3418, 
.1025, .3890, .0765. .0651. 
.1092,. 0996,. 0909, .0830, 
.1239, .1188, .1142, .HOC, 
.1387, .1358, .1332, . 1308, 
.1529, .1511. .1495, .1481, 
.1663, .1652, .1642, .1633, 
.1789, .1782, .1775, .1770, 
.190 8,. 1903, .1899, .1895, 


.7000. .8000, .9000 
.1500, .IGCO, .0500 
.1219, .u3G2f .'j-o97 
.u9o2, .0586, .0286 
.0568. .0364. .0176 
.0374, .0238, .0114 
.0254, .0162, .0077 
.0128, .ooei, .0039 
.gC69. .GO 44 , . u 0 2 1 
.gUnC, .0J25, .G012 
.0024, .0015, .0037 

• 0 G 1 5 f • u 0 u 9 , • 0 0 0 4 
.0009, .0006. .0003 
.7GwO, .6030, .9000 
*15wj, alOuc, .0500 
•0797, .u364,“.uC17 

-•CG17,*.G3l8,*.C 603 
-.0929, -.1107, •.1270 
•.1529,-. 1637.-. 1736 
-.1965, -.2034, -.2097 
. 7 L 0 u , • 8 J U 'j , .9000 

• 15uii, .10 Ou, .0 500 
-.0094.-. 0494. -.C884 
-.1890, -.2159, -.2414 
-.3853, -.3986, -.4111 
-.5157, -.5227, -.5291 

£148, -.6187, -.6222 
.7000, .8000, .9000 
.15C0, .1003, .0500 
-.1336, -.1726, -.2105 
-. 4421. -.4652, -.4fi68 
-. 7616, -.7705, -.7793 
. 7 u u C f • 8 J y C , • ^ w 3 i 
. 15 Cl, .IJjw, .053u 
-.24X9, -.2787, -.3152 
6503, -.67 C5,-. 6892 
.70X0, .8300, .9000 

• 153L, .Igiiu, «L5GL) 
-. 6074, -.6409 , -.6732 
-1.308,-1.320,-1. 331 
-1.931,-1.935,-1.938 
— 2.38s, — 2.3 87 , — 2 . J 8 8 

. 700 0 , . 6 3 WL , .900 3 

• 1 5 u C , .13 3 3, . 0500 
-1.131,-1.158,-1.164 
-2.154,-2.160,-2.166 


i.GOO, 1 

0 . 000 ,- 

3 • C , — . 
J . u , 
0.0 , 

3.0 
3.G 
3.0 
0 • 0 
3.C 
0.0 
G.O 
0.0 • — . 
1.000, 1 
• 0 000 
0437,-. 
-.0872,- 
-.1423,- 
-.1627,- 
-.2154,-. 
1.030, 1 
.0003,-. 
-.1263,-. 
-.2653,-. 
-.4226,-. 

5353.- 
- . 6254,-. 

i.oca, 1 
.0000 . 
-.2471,- 
-.5072,- 
-. 7869,-. 
1.033, 1 
.GGOO,-. 
-.3 53 7,- 
-.7069,-. 
l.ODC, 1 
.0003,-. 
-.7043, 
-1.342^-1 
-1.942,-1 
-2.389,-2 
1 . 000 , 1 
. 0 C u 3 , — 
-1.209,-1 
-2.171,-2 


.3598 
1.2300 
-.1000 
-.059C 
-.0 179 
. G26C 
.0 547 
.0758 
.1C62 
.1286 
.1468 
.1625 
.1 764 
.1892 
ICG 
. u5 0 
0379 
3 266 
016 0 
Qiu 3 
3069 
0 035 
0 019 
0011 
COO 6 
000 4 
0GC3 
.ICC 
G53 3 
3784 
1127 
1563 
1911 

22t:6 

1C C 
0 500 
163 2 
2880 
4331 
5 43 4 
6284 
130 
05CG 
2825 
5264 
^94 2 
. ic: 

0 5 C 3 
3848 
7 235 
10 3 
0 50 3 
7341 
352 
945 
591 
10 C 
0 5 33 
.23 3 
177 


.3589, .3585, .3581 

1.3GCQ,1.40CQ ,1. 5000 
-.1530.-. 2000, -.2500 
-.0948, -.1293, -.1625 
■.0420, -.0646, -.0860 
.0115,-. 0019,-. 0143 
.0452, .0365, ,0285 
• 0692, .0632, .Q57fi 
.1028, .0997, -0968 
.1267, .1249, ,1233 
•1456, .1446, .1436 
.1617, .1611, . 1605 
.1760, .1755, .1751 


i0021 C^2'9,-.0u37 , — « 

• 0012, — . 0018,— .00 22 , — . 
.0008,-. 0011, -.3014,-. 
.0005,-. 0007, -.0009,-. 
\ 1 « 3n n . 1 .lp n . i 


.20C, 1.700, 1.4:0, 
1000 , -.150 0, -.2000 J 
3165,-. 3493, -.3808, 
5445,-. 5615, -.5776, 
-.8013,-. 8074, -,8i33j 
L.2G0, 1.300, 1.4C0, 
.1G3C,-. 150 0, -.20 00 1 
.4177,-. 4494, -.4796, 
.7391, -.7538, -.7676, 
L.20D, 1.30C, 1.4Ca j 
.1000,-. 1500, -.2030, 
7628,-. 7903,-. 81 68, 
. 362,-1.771,-1.380, 
.948,-1. 951,-1.954, 
.391,-2.392,-2.3 93, 
..200, 1.300, 1.400 , 
iG:c,-.i50G,-.2caa, 
,.256,-1.277,-1.298, 
'. 182,-2.186 ,-2.191, 


SPACE, 6, END 


EXECUTIGK 

= li 
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o o 




N = H50:^ 

DO iu I=i(,N 



IF(A83A4l .ME 

« C • ) 

MTEST 


iJ 

COM’TINUE 
IFCJTEST ,EQ. 

1) 

GO 

TO 

IG 

.SOlVt^ STEACY STATE FRCRLEM 





VARPLl 

IFfMTEST .EQ. 

Cl 

GO 

TO 

70 

2ii 

LTEST = KOUM 

= u 





lOUn = JOUM 

= KOUM = 1 





DO 4 ; ITEST=1 

>N 





IF(A8aB+ITEST 

0 £ G • 0 « ) 

GO 

Tr. 

30 


MOUH = HOUM 

FRCNT 

+ 1 





IF(M783+ITEST 

• EG.:) 

LTEST 

= 

3fl 

lOUM = inUN 

+ (M5:0+ITEST) 





JOUM = JCIM 

+ .(M51G + ITEST) 





KOUM = KOUM 

+ (M52:flTEST) 




40 

CONTINUE 
IF(LTEST .EO. 

C) 

GO 

TO 

5u 


STDSTL 
GO TO 2. 

5o WRIT^I (b,3i) (H7ait J, J-1,N) 

F 51 FCRMAK^u Fi?CKT ITERATION COUNTERS -- 2 0 ( I 3, 2 X r ) 

write C6,52) J? N) 

f 52 FORMATC^ FRONT DAMPING FACTORS -- * , 13 ( ^6 .0 y 2X ) > 

DO 6^ I=lfN 
M7ai+r = h 
60 ABhA+I - 1. 

GO TO a: 

70 CONTINUE 
RTOSTt 

80 CONTINUE 

L IMITS (A5J., A5Ib 5G1 .-.l,Tl2C, TIC 1 yRl ,P2 ,PT , R4 , R6 , A 1 U2, Al.j 32) 
C nMAF(5HOOAHS) 


II^C 

00 liJ TTT=lt6 

A2.LufIII~(G'318+III-l) (TlGT^fll) -{T1154-III-1)) /X (T 119+ ITT- 1 + h6C . ) 
+ »^4-(R2w>46:'. ) 


llu * 11=11+6 

PRINTA (6HSINCCN, G3i8,6yi) 

PRI^ TA(6HFASGMA, A2u j J + 1, 6, 1) 

RETUPN 

C,-. SOLVE TRANSIENT FRCfiLEM 
IGO continue 
FWQBCK 
HETPMP 
END 

BCD 3VARIABLES 1 
N = M5CC 

" 1 1 = j . 

t = LL - u 

Ml = Me = Ml = 1 

DO 1030 1=1, M 

IF(ASwfl+I .ME. u.) 1=14 i 

C.. .LIQUID PROPERTIES 

- . IF (I,r,T, 5) GO TO 100 

L 0CKUFfM3C3+I,Ti01+M-l,aiCul,A73l4-Ml) S LinUIO OE'MSITY 
1 OCKUP (••l5i;4l,Tl,e + N2-l, AlJOi, A7114M.e) 

L OCKUP(M52j+I,T10S+N3-l, A1 vu1,A721+M?) 

Lr30KUP(M50C4T,Tiul + M'l-l, A10C3,A7324NH i LIQUin Sr-ECTPIC HEAT 
1 OOKUPt'ISiutl, Ti;24-K2-i, Ai:0 3,A7124N2) 

LOOKUP [”5214 1,T 1:34- N3-1,A10. 03, A722+N3) 

lCCKUP(”5Ll + I,Tli,14Nl-l,Ai005,A7334M) ; LIOUIC VISCrslTY 
L 00KUP(‘f51-4i,Tij24-\'2-i,Alo:?, A713 4M2) 


C-51 




C « • t 


c • • 

1C 


L0CKUP(H52C*I,TiC3+N3-i, AlOO 5 , A723^N3> 
LOOKUP(M5GO<-I.TlQl^Mi-l, AIGD 7,A7 04fNU 
LOOKUP0151C + I,T1J2 + N2-1, A10 u 7,A714 4-N2) 
L00KUP(M52J^I,Tllj3+N3-l, AiCO 7 , A724^N3) 
LOCKUPtM5 00>I fTlOl+Nl-lf AlOO c,A7e5^Nl) 

I 00KUP<tl51C+I, T1^2 + N2-lt AiCC 9 , A715 f N2 ) 

L OOKUPCM52Q + I,T10 34-N3-1» AIGC 9, A725>N3) 
L00KUP<M5CG^-I f T101*N1-1, At 01 Q , A7 06 + Nl) S 
L OOKUPCM51G^I,riC2frN2-l,AiaiO»A7i6+N2) 
LOOKUP{M520^I ,T iO 3 f N3-1, Alyl Q , A726+N3) 
VAPOR PPOPEPTTES 

DlDEGlt T12C + I-1, AlC02,A73H-I) 5 VAPOR 
DlOEGlt T12C^I-1 ,A1Q04|A732+I)S VAPOR 
OlOEGK T12a + I-lt Al006,A7334l) 5 VAPOR 
OiOEGK ri2y+I-i,AluJB,A734fI)5 VAPOR 
GO TO lul 

. HEAOER PIPE PPOFERTIES 
0 L OOKUP(M5aO>I»TlCl^Nl-l,A1021f A7U1 + ND 
T102fN2-l, A1G21» A7llf N2> 
TlC3fN3-i, A1021,A721^N3> 
Ti:ilfNl-l, AlJ2 3,A7a2fNl) $ 
Tlu2fN2-l, A1U23, A712^N2) 
T1C3 + N3-1, A102 3, A722^-N3) 
TiCl+Nl-i,A1025,A7u3fNi) 
T1C2+N2-1, Ai:i2 5,A713 + N2) 
TlC34-N3-lrA102 5, A723 + N3) 
7101 + N1-1,A102 7, A7 04<^N1» t 
T1C24N2-1, AiC27,A7i4+N2) 
TlG3f K3-lf A1027, A7244‘K31 
TlCl^Nl-1, AlQ2c,A7C5+Nl ) $ 
Tl02+N2-irAlJ29,A715+N2) 
Tlu3+N3-lf A132 9, A72P^-N3) 
TlDlfrNl-l» AlC3QtA7C6+Nl) 
Tlw2f N2-1, AiaSG, A716VN2) 
T1C3^N3-1| A103Gf A726 + N3) 


I LIQUID THERHAL CONHUCTIVITY 


9 SURFACE TENSION 


HEAT OF VAPOFI7ATION 


DENSITY 
SPECIFIC HEAT 
VISCOSITY 

THERMAL CCNCUCTIVITY 


LIQUID DENSITY 


LIQUID SPECIFIC HEAT 


3 LIQUID VISCOSITY 


LIQUIO THERMAL CONDUCTIVITY 


SURFACE TENSION 


3 HEAT CF VAPORIZATION 


C • • • 


DENSITY 
SPECIFIC HFAT 
VISCOSITY 

THERMAL CCNDUrTiVlTY 


iJ 

C • • • 


c.. 


L00Kt)P(M5iu^I 
LOOKUPtM520 + I 
lCCKUP(M5oC+I 
L OCKUP<M510+I 
L 00KUP(M52 j^*I 
LOCKUP<M50Cf I 
L OOKUPt M51j + I 
L00KUP(M52C+I 
LOOKUP(M5CCi-I 
L 00KUPCM5lGfI 
LOOKUP(M52:+I 
LOGKUP<M5Cw^H 
L OOKUPIM5ia+I 
LOOKUP{M520+I 
LOOKUP(M500+I 
L 00KUP(M31G + I 
LOOKUP(N520^I 
VAPOR PROPERTIES 

DIOEGK T12C4l-l,Alu22,A73l4-I) ? VAPOR 
OlOEGU T12C4l-lf A1G24,A732+I)S VAPOR 
OIOEGK T120 + I-1, A1026/A7334-I) 3 VAPOR 
DIOEGK T12KI-i,Ai:23jA734f 1)3 VAPOR 
1 CONTINUE 

COMPUTE ARTERY CAPACITANCES AND CONDUCTANCES 

ARTERY (M50 0 + 1,1191 04^ I, H62l^N 1,A621+N1, A 6224N1,A6 31 + Nl , A623+N1 , 

A624 + N1,A6 2 5+N1,R1, A6 26 + N1, ABEZ^-Nl, A628+K1, A6 29+K1, 

A63£l + Ni^A70 4 + Ni-#A53 7+Nl,A701+Nl,A?a<l+Nl» A5 35+N1, 

A6814N1 ,A6 84 4-N1 ,A6d7+Nl) 

ARTERY (M51C+I,M9i24l,M64y+N2, A64i+N2,A642+N2,A651+N2,A64T+N2, 
A6444N2, A645 + N2 ,R1,A6 46 + N2, A647 + M2 , A646 + N2 , A649 + N2, 
A65D4N2,A7144N2,A557+N2tA711+N2,A712+N2t A955+K2, 

A6 82 + N2, A6 85+N2 ,A6a8+N2) 

ARTERY(M5aG+I,M9i4+I,M66C+N3 ,A661+N3,A6624N3,A671+N3,A6634N3, 

Ao 64 + N3, A665+N3 ,Ri, A666+N3, A6674N3, A6684N3 , A6694N3, 
A67C+N3t A7244N3,A577+N3,A721+N3,A722+N3, A5754K3, 
A6ft3+N3,A686+NM,A689+N3) 

COMPUTE CAPACITANCES 

CAPFIN(M52l4N3, A612+N3,A611+N3, A 52 6+N3 , A6 164 N3, A617+N3, A850 +N3) 
CAPWCK(M5 Cu+I, A5-ai4M,A5Q2 + Nl,Rl, A5044Nl ,A5C5+Nl,A5G64Nlt A5384Nl, 
+ A7D14N1, a7u24Ni,A53l4Nl,A5324Nl,A534+Nl, A535+N1, A5374N1, 

^ A6814N1,0*,CAPE) 

CAPWCKCMpiOfi, A5114N2, A5124N2,R1, A514+N2, A5154N2, A5164N2, A5584N2, 
♦ A7U+N2, A712+N2 ,A55i+N2,A55 2+N2,A5544N2,A5 55+N2, A557+N2, 

+ A6824N2^G.,CAPA) 

CAPWCK(M52u 4I, fl52l4N3, A522+N3,R1,A5244N3 ,A5254N3, A526+N3,A57B4N3, 
4 A7214N3, A7224N3>A571+N3, A572+N3,A5744N3, A5754K3, A577+N3, 

+ A683+N3, A8534M3,C1034N3-1) 

C1^.34N3-1 = C1G34N3-1 4 CAPE 4 CAPA 
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C... COMPUTE AXIAL CCNCUCTANCiS 

15 AXC0N0(‘T5u^ + I ,A5 J3»N1, A5o4+Mlf A5C5<-N1, A536+N1,R1, A5374-N1, 

♦ A533 + M. A533+N1,A5344-N1,A53 5 + N1 » A653+M, A696+K1) 
AXCCND(H51C+ I. A513+N2, A514 + N?.,A515fN2, A516+M2,R1» A557«-N2, 

• A556«•^>^, A55 3+N2,A55Af N2, A555+N2, A594«-N2, A6 97+N21 
AXCONDtMScJ+IjASETtNS, A5 24+M 3 , A525+N 3 , A526tN3,Ri , A577+N3, 

♦ A5 78+N3, A 57 3«-N3 ,A574«-N3» A57 5+N3 ,A695<- N3, A698+N3) 

C 

53 GO TO (2Cl,2j2,2G3,2aH,2J5,2G6,221) , I 

2ul RADFINCM5 2u+I, ti.3,Ti^3,T2J3,A613 + N3,Ao12 + N3, A6114-N3,A526<-N3, 

* P5, A6„ j, A6J.‘t+N3,Abl3<-N3, A:'b0 + N3,A613«-N3, GIO, A851 + N3, 

* T13»A773*N3, A311+I) 

GO TO 213 

20 2 RAOFIN(M5 20+I,T20 3, T1 03tT303, A613I-N3, A612 + M3, A611+M3,A526tN3, 

* R5,A5.J,A&lH4-N3,A615fN3,A7b8tN3,A61C*N3,G23,AR5i+N3, 

» T2C jA77G+N3, A411+I) 


GO TO 210 

2C3 RAOFIN X t'5 23 + 1 , T 3 v3, T 2 0 3 , T4 03 , Aol 3+ N3 , A612 t-NS , A6 11 +N3 , A5 26<-N3, 

* 9 5, A6uC, A614+N3, A615+N3, A7b8 + N3,A&lu+N3,G?Gsi A85H-N3, 

• T33,A77C + K3, A81H-I) 

GO TO 210 

2C<* RA0FINtM52j4l ,T40 3,T3 J3,T5 03 ,A613^N3, Aei2X-N3, A611 + N3,A52f,+N3, 

♦ R5, A6C 0, A614»-N3,A615+N3, A76 8 + N3,A61C»-N3,G4 0, A85x+N3, 


* T4C ,A770i-N3,A81i+I) 

GO TO 210 

2C5 RAOFIN(M52o + I ,T 50 3 , T43 3 , T6u 3 , A61 3+N3 , A61 2 + N3 , A6 11 tNT , A5 26 + N3, 

* R5,A6C«,A614«-N3,A615<-N3,A7fi8+N3,A61C + N3,G50, AP51I-N3, 

* T53,A773+N3, A811+I) 

GO TC 213 


23 6 PA0FIN<M5 2j + I,T6C3,T5-3,T603,A6l3<-N3,A612fN3,A611 + N3, A5 26+N3, 

* R5, A63 3, A614+N3 , A615+ N3 , A76 8 + N3 , A 610 + N3, G6 C , A851+N3 , 

* T63, A770+N3, Abll+I) 

213 CONTINUE 

221 CONTINUE 
C 

AGCCM9(M5J: + I ,«510+T, "523 + 1, A684 + N1 , A6 85 + N2r A68 6 + N-?, A69 3+N1, 

+ A69A+N2, A695+N3,A696+N1,A697+N2,A698+N3, A618+N3, 

■* A851+M3,G311+LL) 

C... COMPUTE LIOUn/KICK EFFECTIVE RAOIAL THERMAL CONOUCTIVITIES 
eFFKWK(M5G3+I,M533+Nl, A533+N1, A734+N1, A5 0 3 +N1, A53 8+Nl, 

* A534+M1, A535+M,Ri, A542+N1, A54 3+N1, A5 44+N1, A765 + K11 
EFFKWK(i5i; + I, "55iJ + N2 , A5E3+N2, A714 + H2, A513+N2,A558 + N2, 

* A554+N2, A555+N2,R1, A562+N2, A563+N2, A564+N2,A766+N2) 
EFFKWKtM523+I,K57a+N3, A573+M3,A724+N3, A523+N3,A578+N3, 

* A574+N3,A575+N3rRl, A562+N3,A583+N3, A584+N3,A767+N3) 
C... COMPUTE RAOIAL CONDUCTANCES 

RACOND(H5t>i + I,Rl, A7 65 + M , A53 7 + N1, A534 + N1, A535+N1, A5 l’+N1, 

* A 5 3 6+K‘l , A 5 j s+N i, A 50 5+ Nl , G 1 0 u 1 + N !■ 1 ) 

PAC0N0<VI513 + I,R1, A7 66+K2, A557 + N2, A554 + N2, A555+N2, A513+N2, 

’ A516+N2, A5l4+N2,A515+N2,G1002+N2“l) 

RACOND(M52u+I,Rl , A767+N3, A577+N3, A574+N3, A575+N3, AS23+N3, 

* A526+N3, A524+N3, A525+N3,GiaC3+N3-l) 

GO TO (32 : .,321, 32 2, 3 2?, 32 4, 325,30 0) ,T 

320 FAC = 18.6 
GO TO 330 

321 FAC = 18.6 
GO TC 333 

322 FAC = 18 .6 
GO TO 33o‘ 

323 FAC = 18.6 
GO TO 330 

324 FAC = 18.6 
GO TO 33ii 

•325 FAG = 18.6 

330 Glj£l+Nl-1 = GlCOl+Ml-l /FAC 
G1CC2+N2-1 = G1002+N2-1 / FAQ 
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00 3ul 11=1,6 


301 

G10C-3+H3TII-2 = GlJ03*N3 + II-2 / 
GO TO 310 

FAC 

3G0 

00 303 li=l,A 


3C3 

G1001*Nl*lI-2 = GlOuH-M + II-2 / 
G1002TN2-1 = Gl3u24-N2-1 / 14A.2 
DO 304 11=1.6 

144.2 

304 

G1003»N3+II-2 = Gia03+N3TlI-2 / 

144.2 

310 

CONTINUE 



C.«.C0HPUfE FEEOTUBE CONDUCTANCES 


JJ = K52CfI 


'qSCONO( A8GR+I,Rl,A5 26^JJ,A52 3f JJ, A524^-JJ, A52 5+JJ, A8 26<*I, 

* A82G^I, A825^-I, Aa23 + I, Afl22fI,A85 7^I, A8 2S + I,A8 27^I, 

♦ A577 + JJ, A57 8^JJ,A57 3+ JJ, A57 4+JJ,A57 5f JJ, A82 84-I, 

> A824+I,r,3i8 + L-1) 

C. •.ADJUST CONOENSEP COUPLINGS 8ASE0 ON LA/LC 

AO JG{ A3C8f I,M52Cf I, A526+N3 ,A^e2Q4-I , Gig J 8 + N 5-1 ) 


C.,.C0MPUr£ HEAT LCAOS 

Rll = Rll ^ A8114I 

N1 = K1 <H50G + n 

N2 N2 + (M51G+I) 

N3 = N3 f (N52:+I) 

LL = LL f (H50a^-r) f (H51C + I) 4- (M52G4-I) - 1 

lOJO continue 

0731 = 0703 = G7D5 = 0707 = 683.21/4. 
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RCO 3 VARIABLES 2 

IF(JT£ST .EO. a .OR. H1EST .EO. j) RETURN 
N = M5:C 

MOUM = 3 

lOUH = JGUN = KOUN = 1 

DC 5u ITEST=1,N 

IF tA3u8+IT£ST .EC. u. ) GO TO 40 

MOUM = NQUK 4-1 

FRONT 

IF(N783+IT£ST .EC. 0) BACKUP = 1. 

40 ICUM = lOUK 4- (H500 + ITEST) 

JDUM = jnUM ♦ (MEiOfITEST) 

K.DUM = KOUK f (M5 234 -ITcST> 

50 CONTINUE 

IF(BACKUP .NE. C.) t^ETURN 

ICO WRITE (6,51) (M7814J, J-i,N) 

F 51 FORt1AT(^J FRONT ITERATION COUNTERS — * , 2u ( I o , 2 X ) ) 
WRITE (6,52) (A8444-J,J=1,N) 

F 52 FORMAT(*- FRONT OAKFING FACTTRS -- ♦ , 1 3 (F6 . 0 , 2X ) ) 
DO 11 J 1 = 1, N 

M7 811*1 = i 

110 A8441-I = i. 


STGREP 

ENC 




t »¥ »¥ *¥* 


BCnSOUTPUT CALLS 
TPPINT 

PRIMTAISHTEMVIR, A7&BH-i ,M768t 1) 

PRINTA(6HQLATKT , A7S2<i 1) 

PRINTAUHGPvJrTN, A77 j + l,M77G, 1) 

PRINTA(6HQRJCTP, ABlH-l,Mfill, 1) 

WRITE(6,1G) Pli 

F 10 FORf'ATl*: TCTAL PANEL h£AT R E JECTI CN =♦, F15. 3,* 0TU/HR*) 
IF(JTEST .EO. 1) GO TO 20 

STOREP 
RETURN 

C... COMPUTE HEAT TRANSFER LIMITS 


23 L IFITST A5'.;, A5:C , G 1 OG 1 , T 1 20 t T 1 J1 ,Ri, R2,R3 , PA , R6 , Ai: 12, Al.) 32) 
FNO 
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RCO 3 sun ROUTINES 

C...IN several cf- the fcllcming suhrcutines, multiple cases can re ccnstderfr 

C....WITH ONE CALL, SINCE ALL VARIAELES ARE SURSCRIPThO. THIS OPTION IS USED 
C....BV INPUTTING THE VARIARLES IN ARRAYS WHERE EACH ARRAY roNTAINS MULTIPLE 
C.... VALUES, REPRESENTING EITHER SEVERAL °IPES CR SEVERAL NCDES IN ONE PIPE. 
C...,IF ONLY ONE CASE IS OESIREO, INPUT N=1 IN THE SUBROUTINE GALL. IF SEVERAL 
C... .CASES ARE DESIRED, ENTER THE INTEGER COUNT OF GKE CF THE ARRAYS. THE 
C....SUPRCUTINES CONTAINING THIS CPTION ARE EFFkrWN. ARTERY, RACONO, AXCPND, 
C....CAPWCK, CAPFIN, RAOFIN, AND LOCKUP. 

FSTART 

* SUPPOUriNE EFFKWKLN,NWICK,KMICK,KL , KW ALL , P CR OS , CD , IC , PI , 

GNQ, WIDTH, DEPTH, EFFK) 


THIS SUBROUT I 
OF THE LlOUin 
AND HALL TS G 
WHF.RE POOR CC 
RACONO, THUS 
N 

NHICK 

KHICK 

KL 

KHALL 

POROS 

on, ID 

GNO 

WIDTH 

DEPTH 

EFFK 


NE COMPUTES THE EFFECTIVE RADIAL THERMAL CONDUCTIVITY 
-SATURATED WICK. IT ASSUMES THE CCNTACT BETWEEN WICK 
COD, WITH THE EXCEPTION OF THE AXIAL SPACFR WIRES 
NTACT IS ASSUMED. EFFK IS USED IK SUDOCUTINE 
EFFKHK MUST BE CALLED BEFORE RACCKO IN VARPLl. 

= NU'MBER OF CASES TO BE TREATED IN THIS CALL 
= TYPE OF WICK 

= THERMAL CCNQUCriVITY CF WICK "ATEPIAL 
= THERMAL CONDUCTIVITY OF LICUID 
= THERMAL CCNDUCTIUITY CF WALL MATERIAI 
= HICK POROSITY 

= WICK CUTER AND INNER OTAMETERS 
= NUMBER OF GROOVES, WIRES, OR THREADS/FT 
= GROOVE OP. THREAD MEAN WIDTH 
= C-ROCVc OP THREAD DEPTH, OR WIRE DTAM 
= effective thermal CCNnuCTTVTTY 


REAL KWICK, KL , KWALL, ID 

DIMENSION NWICK fl) , KWICK (1> , KL U) , KWALL (1) , POROS ( 1) , CO ( 1» , ID ( 1» 
DIMENSION GNO (i),HIDTHU) ,0EFTHU),EFFK(1) 

DO 2 L u i, 1= i , N 

IF(NHICK{I) .LE. J .OR. NWICK(I) .GF. 6) GO TO ICOO 
J = NHICK(I) 

GO TO (iu’O.Si J.30'J,4Gu,5GO), J 

...WRAPPED SCREEN, SIKTEREO POWDER, FELT, CR FCAM CN IC OF PIPE 
IJO EFFK(I)= KL (I)^(KL(I) fKHICK(I) - ( 1 . -POROS ( I ) ) » (KL ( T ) -KWICK ( I » )» / 

* (KL<I) #^KWICKtIJ ♦ (l.-PORCS(Tn»(KL(I)-KWICK(IV)) 

IFCNWICKd) .EC. 3) GO TO 310 

IF (NWICK (II .fq. 4) GO TO itlJ 

GO TO 2400 

...OPEN AXIAL GROOVES ON ID OF PIPE — APPROXIMATE COD-SHAPED 
....GROOVE X-SECTICNS WITH ECUIVALENT SDUARE LENGTH AN^ DEPTH 
20 0 EFFK(I)= GNCdl »WIOTH (I) * ( KL ( II -KW AL L ( I) ) / (R I* ( CD ( I) -Q EPTH ( T) )l 

* > KWALL (I) 

IF(NWICK(I) .EQ. 3 ) GO TO 32Q 

GO TO 2000 

C...SCREEN-COVEREO AXIAL GROOVES ON ID OF PIPE 
300 GO TO iOj 
3lC EFFKl = EFFKtll 

GO TO 200 

320 EFFK2 = EFFKdl 

FFFK(1I= 1. /(I. /EFFKl + 1./EFFK2) 

GO TO 20C0 

C. ..SCREHN-CCVERED AXIAL SFACEP WIPES 
400 GO Tn 100 

410 EFFKW = EFFKdl » 2.* PI * ID(I)/(OC(I) - ID(I) - 2.» DEPTHdH 

EFFKL = KL(II * (PI * (OD(Tl/DEPTH(I)-i.) - GNCdl I 

EFFKS = KLd» » KMALL(I1 /(. 1738 * KHALL (I) + .8262 * KLdll 
EFFKd)= 2.*(DD(I)-ID(I))/(PI’f(ODd) + l0d» )* (1 . / ( EFFKSi-EFFKL I 

* * i./EFPKW) ) 

G 0 T 0 2 0 4 L 

C...C1RCUMFEPENTIA1, GROOVES ON TO OF PIPE (THREADS) 
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WIDTh/I) ♦ (KL(II-KHALL<in ^ KWALUI) 


SCO EFFK<I1= GMCa) • 

GO TO 20QC 
IQuO WRITE(6,lu0i) 

iOOl FORMATf^J + + IN\)ALI0 WICK TYP£,^» 

• ♦ RAOIAl KICK CCNCUCTCRS ERRONEOUS^) 

ZQOu CONTINUE 

RETURN 

="NO 

SURPOUTINE PERHRL(^WICK,00,Ir♦CN0t WinTH,OEPTH,KART,nrART> 

* CIAR7, WTfA ART, K,NT,SO,TfPERW,PFRMTA,PI,DR) 


THIS SUPROUTINE COMPOTES THE PERMEABILITY ♦ CROSS-SECTIONAL AREA 
PRODUCTS CF THE ARTERY <IF ANYI ANO WICK. SINCE THE ARTERY *ANO 
WICK FORM PARALLEL LIQUID FlCW PATHS, THEIR PERMEABILITIES ARE 
SUM*ME0 LINEARLY. PERMTA IS USED TY SUBROUTINE LIMITS AND PERMBL 


IS CALLED RV 

limit: 

NWICK 

= 

OD,IO 

= 

GNC 

■= 

WTOTH 

- 

DEPTH 

= 

NAPT 

= 

OOART 


DIART 

= 

WT 

= 

AAPT 

= 

K 

= 

NT 


SO 

=■ 

T 

= 

PERM 

- 

PERMTA 


OP 

= 

PIKENSION AART(l) 

REAL ID, NT, 

KAC 

PERMA = Q. 


IF (MART .£0. 

0) 

IFCNART .IT. 

j • CR 

GO TO (i3,2i) 

,NART 


TYPE OF WICK 

WICK OUTER ANO INNER DIAMETERS 

NUMBER OF GROOVES, WIRES, CR THREADS/FT 

GROOVE OR THREAD MEAN WIDTH 

GROOVE OR THREAD DEPTH, OR WIPE OIAM 

TYPE OF APTEFY 

ARTEPY OUTER DIAMETER 

ARTERY INKER OIAKETER 

THICKNESS OF OUTER ARTERY WRAP 

ARRAY OF ANGULAR WIDTHS OF SPIRAL ARTERY ANNULI 

NUMBER OF SPIRAL ARTERY ANNULI 

NUMPER OF SPIRAL ARTEPY TURNS 

SPIRAL ARTERY SPACER WIRE DIAPETEP 

SPIRAL ARTERY SCREEN THICKNESS 

HOKCGENEOUS WTOK PERMEABILITY 

PERMEABILITY » CROSS-SECTIONAL AREA 

OEGPEES-RAOTAKS CCNVERSICN (57.3) 


GC TO 9C 

5£. 3) GO TO 10 u2 


128 


G, •.PEDESTAL ARTERY 

10 PERMA = PI ^ 0IART»»4 / 

GO TO 90 

C... SPIRAL ARTERY WITH SPACER WIRES BETWEEN SPIRALS 


20 


DART 

PI 


= OOART - Z.* WT 

= niART/2. 4^ T + (0ART-0IART-'2. n) AART(i) / 

(fl.» PI ♦ NT DP) 

= PI + SO . 

= AART (1)V (2.*OR) - SO/fRO + RI) 

CALL ANNSPCCRI, PC, THETA, KAO 
PERMA = KAC 
DO 25 1^2, K 

= RT 4- {DART-OIART-2. »T) > (AART(I-l) ^ AflRT(T)) f 
(6.^ FI » NT ^ DR) 

= RI + SO 

= AART(I)/(2.»OR) - ST>/(R04-RI) 

CALL ANNSPC(RI, RC,THETA,KAC) 

FERMA =: PERMA + KAC 
COMINUE 

PERMA ^ PEPKA + FI + (DIART / 2 . ) 4 
IFCNWICK .LE. 0 .OR. NWICK .GE. 6) 

GO TO (ICO, 200,3.jQ,4Ga,500) rNWrCK 
C...WPAPPEO SCREEN, SIKTEREO PQWDEn, FELT, OR FOAM ON ID OF PIPE 
ICO PERMTA = PERt" ♦ FI ♦ (CD»^2 - ID^^2) / 4. 

RETURN 


RO 

THETA 


RI 

PO 

THETA 


25 


90 


/ 8 . 

GO TO IDCO 


+ PERMA 
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...OPEN AXIAL r,RCGVES ON 10 CF PIPE — APPROXIMATE COO^SHAPED 
....GROOVE X-SECTTCNS WITH EGUIVALENT SQUARE LENGTH ANC OEPTM 
2 j- PERMTA = 5.* GNO * MI0TM**3 * DEPT H** 3/ ( 16 . ♦ (WIDTH**? 

* t w.* 0FPTH**2)> * PERMA 
PETUPN 

...SCREEN-CnVEREO AXTAL GROOVES ON ID OF PIPE 
3C0 PEPHTA = 5.* GNC * WI0TH**3 * DEPTH**3/ C’2.* (WIDTH**2 

* + CEPTH**2n * PERKA 
RETURN 

...SCREEN-rOVEREO AXIAL SPACER HIRES 
4u3 RO = on / 2. 

PI = PO - DEPTH 

THETA = (PI * (RO*'’I) , GNO*DEPTH)/f GNO * (RO^RD) 

CALL ANN^Pn(RI,RC, THETA, KAC) 

PERHTA = GNO * KAC <■ PERHA 
RETURN 

...CIPrUMFERENTIAL GRCCVES ON TO OF PIPE (THREADS) 

5jJ PERHTA = 5.* WI0TH**3 » 0EFTH**3 ,'(GNC * PI * (CD + IO) * 

* 9.* (W10TH**2 + 4.* DE°TH**2)) f PERHA 
PETURN 

luJu WRITF(6,luul> 

lOJl FQP.PAT(*u ++*+t£RRCRiM-'M-* INVALID WICK TYPE,*, 

* * PERHEAEILITY ERPONE CUR*) 

RETURN 

1002 WRITE(6,iai3) 

1003 F0RHAT(*3 <-+«■ ttC AUT ION ++ T++ INVALID ARTERY TYPE,*, 

* * NO ARTERY ASSUMED*) 

GO TO 9C 

END 

♦ **•*♦**♦**♦<(1 •»♦«»*»•»-**♦*♦♦*****•#»*•*»«»**»». ♦*»->M,***¥»»»»*** 

SUPROUTINE ANNSFC (R 1 , R C, THET A , KAC) 


THIS SURROUTINE COMPUTES THE PERMEABILITY * CROSS-SECTIONAL AREA 
PRODUCT OF AN ANNULAR SPACE. IT IS US^O PY SURPOUTINF PERROL AND 
IS CALLED PY FERKSL INTERNALLY. 

PI,RO = ANNULAR SPACE INNER AND OUTER RADII 

THETA = half-angle WIDTH OF ANNULAR SPACE 

KAC = FFRHE ABILITY * CROSS-SECTIONAL AREA 


o£AL KAC 

KAC = 5.* THETA**3 ♦ (RO**A - PI**4 - 

* 2.* (PC**3 » RT - PO * RI**3))**2 / 

* (Ifi.* ((4.* TH£TA**2 * 5.) * (R0**4 - RI**4) - 

* 6.* (THETA**2 + 5.) * (R0**3 * RI - RC * ^1**3) *■ 

* 6^.* RC**2 * RI**2 * ALOG (RO/RI) ) ) 

RETURN 

END 

SUBROUTINE ARTERY (N,NA1,NART ,OOART,OIART, WT, SD,T, POPCS,PI,NT, 

* RHC0FA,RH0CFS,KA,KS,KL ,L ,D ENL ,CPL, I0,CAP, 

* CONOjEFFVO) 


THIS SUBROUTINE COMPUTES THE CAQACITANCE AND AXIAL CCNOUrTANCE 
OF AN APTFRY. RADIAL CCNDUCTANCE IS UNNECESSARY. BOTH THE 
ARTFRY MATERIAL AND THE LIQUID FILLING IT ARE CCNSIDERED. THE 
CAPACITANCE AND THE CGNOUCTANCE ARE AnOEO TO THOSE OF THE HICK 
IN SUBROUTINES CAPHCK AND AXCCNO, THUS ARTERY MUST BE CALLED 
BEFORE THOSE ROUTINES in VAFBLl. TH£ EFFECTIVE VAPOS SPACE 
DlAMfTER RESULTING FRCP THE PRESENCE CF AN ARTERY IS ALSO COM- 
PUTED FOR USE IN SU4R0UTINP LIMITS. 


N 

= 

N’UMPER OF cases TC BE TREATED IN THIS CALL 

NAi 

=: 

NUMBER CF spIpaL ARTEPY ANNULI 

NART 


TYPE OF APTERY 

00 APT 

=: 

ARTERY OUTER OTAMETER 

OIART 


ARTERY INNER DIAMETER 

WT 

= 

thickness OF OUTER ART£py WRAP 
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SD 

T 

PCROS 

MT 

RHOCPA 

PHOCrS 

KA 

KS 

KL 

L 

OEM 

CPI 

ID 

CAP 

CONO 

EFrvO 


SPIRAL ARTERY SPACER WIRE DIAMETER 

SPIRAL ARTERY SCREEN THICKNESS 

APTERY SCREEN POROSITY 

NUiBER OF SPIRAL ARTERY TURNS 

ARTERY SCREEN DENSITY * SPECIFIC HEAT 

SPACER WIRE OrNSITV ♦ SPECIFTC HEAT 

ARTERY SCREEN THERHAL CCNCUCTIVITY 

SPACER WIRE thermal CONDUCTIVITY 

LIQUID THERMAL CONDUCTIVITY 

NODAL LENGTH ALONG PIPE AXIS 

tlQUin DENSITY 

LIQUID SPECIFTC HEAT 

WICK INNER DIAMETER 

ARTERY CAPACITANCE 

ARTERY AXIAL CONDUCTANCE 

EFFECTIVE VAPOR SPACE DIAMETER IN PIPE 


PEAL L, LSCRNi 10, <A , KS, KL, NT 

DIMENSION NART(l) ,DQART(1> ,0IART(1) , WT (1 ) . SO ( i) , T < 1 » ,*>OPOS <1) 
DIMENSION NT(1) , RHOCPA (1 ), RHOCPS tiJ , KA (1) , KS ( 1) ,KLUI,t (1) 
DIMENSION DENL (1> ,CPL (il ,10(1) ,CAP (1) ,COND (1) ,EFFVD(1) 

DO 2CJ0 1=1, N 

IF (NARK ]> .NE. o) GO TO 5C 

C...NO ARTERY 

10 CAP(I) =0. 

CONDtI)= l.E-12 
ARTXSA = 0. 

£FFVD(I)= ID(I) 

GO TO 20.: 

.50 IF(NART(I) .LT, u .09. NART(I) .GE, 3) GO TO 1000 
J = NARK II 

GO TO CiaG,2C0),J 
C. ..PEDESTAL ARTERY 

IDd VSCPN = PI^(0CAKT(I)»»2 - 0 I APT (I ) 2 ) »L ( T ) ^ (i . -P CRCS ( D) / A . 
VLIO = PI ^ OOARKII^^Z ^ t(I) / 4. - VSCPN 
CAP(I) = VSCRN ^ RHCCPA(I) ^ VLIO ^ OENL(I) CFLtl) 

COND(I)= KA(I) ^ VSCRN / (2. ’ LCD^-^a) 4- 

♦ KLd) ♦ PI ♦ 0IART(I)^^2 / (4.^ L(D) 

ARTXSA = FI » D0ART(II»»2 / 4. 

EFFVD(I)= KFI IC(I)»^2/4. - ARTXSA) ♦ 4./ FT)^».5 
GO TO 20 CO 

C. ..SPIRAL ARTERY WITH SPACER WIRES BETWEEN SPIRALS 
200 DART = DOARKI) - 2.» WKI) 

A = (DART « DIAPT(I) - 2.^T(I)> / PI ♦ NKIT) 

B = BIARKI)/2. + TCI) 

LSCPN = (2.-^ PI ^ NT(I) ^ B/A) » ((2.^ FI ♦ NT ( I) » A + P) *^2 ^ 

♦ A^»2)^».5 - P/A ♦ *■ + A/2.^ 

» ALCG(C<(2.» PI * NKI) ^ A > 8)^*2 f > 

♦ 2 .* PI ^ NKII » A 4 0) / ((B^^2 + A^^2)*^^.5 + B> ) 
VSCRN = LSCRN T( I) * L (I) ♦ ( i . -PCRQS ( I ) ) f PJ tOOARKI)»^2 

^ - OART»*23 ♦ L(I) * (i.-POROS(I) ) / 4. 

VSFAC = FLOAT (NAi^'il ^ PI ♦ Sn(i)»*f2 ^ L(I) / 4. 

VLIO ■= FI * CCART(I)^»? * id) / 4. - VSCPN - VSPAC 

CAPd) = VSCPN • RHOCFA(^) + VSPAC » RHGCPSd) + VLIO ^ DENLd) 

^ CPL (I) 

CONO(I)= (VSCRN/2.» KA(I) + VSPAC ♦ KS(I) 4 <PI * 0flRT»*2/4. - 
» LSCRN ^ Td)) ^ LCD * KLCIM / L(I)*»2 

ARTXSA = PI » DCARKI)»^2 / 4. 

£FFV0(I)= ((FI » I0(I)^/^2/h. - APTXSA) ^ 4./ FI)^*,E 
GO TO 2000 
lOOfl WRITE{6,i0ul) 

1001 FOPMAK^u +44 + 4GAUTICN4+4 + 4- INVALID ARTERY TYPE , 

♦ * NO ARTERY ASSU^^ED*) 

GO TC 10 

2000 CONTINUE 
RETURN 
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END 

SUBROUTINE LOCKUP C N |T , A, p ) 


THIS SURROUTINF OBTAINS PROPERTIFS AS A FUNCTION OF TEMPERATURE 
ON AN AUTCMATFO BASIS. THE TKOEPENOENT AND nEPENOENT yAPIARlES 
MUST BE INPUT SE CUE KT I ALL Y . 

N = NUMBER CF LOCKUPS IN THIS CALL 

T = INOEPENOENT VARIABLE ( TEMPERATURE) 

A - ARRAY OF T,P VALUES 

P “ OEPENOENT VARIABLE fP^OPERTY) 


niHENSlON T(l) fPil) 

DO 10 1=1, N 

. 10 CALL OIDEGKT(I) ,A^P(D) 

RETURN 

ENC 

^ if. If. Jf. if ififif ififw.ifJfifJHf.9^ifif If if if if if if if if.. if if if m if if if .if if if H. if if ifififififif 

subroutine CAPWCK(N,RHCWA,CPWA,PI,0DHA,I0WA,LHA , 

» POROSfRHOL ,CPt ,RHCWK,COWK,OOWK,IOWK, 

• LWK, CAPA, CAPF, CAP) 


THIS SUBRCUTINc COI'PUTES THc CAPACITANCE OF THE L tOUIO-SATlIPATEO 
HICK AND PIPE WALL. THEY APE SUHHEO ALONG WITH THE AR-^£RV CAPA- 
CITANCE FROH SURPQUTINE ARTERY ANO THE FIN CAPACITANCE FROH SUB- 
ROUTINE CAPFIN (IF ANY) TO OBTAIN THE TOTAL CAPACITANCE CF EACH 
PIPE NODE. THE CAPACITANCE CF THE VAPOR IS TGKCRED. 


N» 

- 

NUMBER OF CASES TO PE TPEATEO IN TMjc CALI 

PHCWA 


PIPE WAIL DENSITY 

rPWA 

= 

PIPE WALL SPECIFIC HEAT 

OnWA. TDNA 


PIPE wall outer ANO INNER DIAMETEPS 

1 WA 

=; 

PIPE WALL NOOAL LENGTH 

PC^OS 


WICK POROSITY 

RHCL 

- 

LIQUIO DENSITY 

CPL 

- 

LIGLIO SPECIFIC HEAT 

PHCWK 

- 

WICK DENSITY 

CPWK 


WICK SPECIFIC HEAT 

OOKK, ICWK 


WICK OtlTER ANO INNER DIAMETERS 

L WK 

= 

WICK NODAL LENGTH 

CAPA 


ARTERY CAPACITANCE 

CAFF 


FIN CAPAriTANCE 

CAP 


TOTAL CAPACITANCE 


REAL IOWA, LWA,ICWK, LWK 

niKFNSION RHOKA (1) ,CPWA( 1) ,OCWA(i) , IOWA (1) ,L WA U) ,PCROS (1) 
OIMFNSIOK PHCL(l) ,CPL(1) , RHOHK (l) , CPWK (1 ) , OOWK( 1 ) ,IOWK(l) 
OIMENSrON LWK (1) , CAPA (1 ) , C Ap F ( i) . C AF (1 ) 

00 10 1=1, N 

CAPWA = RHOWA(I) * CPWA(I) ^ PI ^ (CDWA(I)**2 - I0WA(I)»*2) 

• < * LWAa)/4. 

CAPWK = (PCROS(I) » RHOL(I) * CPL(I) *■ ( 1 . -pCRCS ( I ) ) 

* * FHCHK(I) * CPWKfl)) » PI * (OOWK(I)»*2 - IDHK(I)»»2) 

♦ LWK(I)/4. 

Ij CAPCII = CAPWA f CAFMK ♦ CAPA(I) + CAPF(T> 

RETURN 

ENC 

SUBROUTINE CAPFIK(N,X, Y, 7,0EN,CP,CAP) 


THIS SUBROUTI 
FIN (IF ANY) , 
N 
X 
Y 
Z 

OFN 


NE POHPUTES THE CAPACITANCE OF A SPACE RAOTATOR 
ATTACHED TO A HEAT PIPE CONDENSER. 

= NUiPER OF CASES TC "E TREATEG TN "^HTS CALL 
= FIN LENGTH (HALF PIP£-TC-PIPF DISTANCE) 

= PIN THICKNESS 
= fin width CCCMOFNSFR LENGTH) 

= FIN nFNSlTY 
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CP = PIN SPECIFIC HEAT 

CAP = FIN CAPACITANCE 


DIMENSION X(U,Y(1» ,2(1) ,0£N(l),rP(ll ,CAP(1) 

DO lu 1=1, N 

10 CAP(I) = 2.* X(l) • Yll) • 7(1) » 0£N(I) ^ CP(I) 

RETURN 

subroutine AXCONC(N,KHA,nOMA ,I0HA,LWA,PI ,L HK , POROS , KWK , ODWK , 
* IDWK,CONDWA,CONOWK) 


THIS SUPROUTINE COMPUTES THE TOTAL AXIAL CONDUCTANCE OP THE, 
HICK AND HALL. NEGLIGIBLE CONTACT CONDUCTANCE IS ASSUMED 
0ETHEEN CROSSEC MEMBERS IN HCHOGF.NECUS H1CKS--THE WICK VOLUME 
IS CONSEOUENTLY HALVED. THF CONDUCTANCES ARE CCMBINED WITH 
THOSE FROM SLDRCUTINE ARTERY IN SUBROUTINE AGCCVR, THUS ARTERY 
AND AXCOND MUST BE CALLED BEFORE AGCOMB IN VARBLl. 

N = NUMBER OF CASES TO BE TREATED IN THIS CALL 

KWA = WALL THERMAL CONDUCTIVITY 

OOWA.TOHA = HALL OUTER AND INNER DIAMETERS 

LWA = WALL NODAL LENGTH 

LWK = WICK NODAL LENGTH 

POROS = WICK POROSITY 

KWK = HICK thermal CONDUCTIVITY 

0nWK,I0WK = WICK OUTER AND INNER DIAMETERS 

CONDKA = HALL AXIAL CONDUCTANCE 

CONOWK = WICK/LIQUIO AXIAL OOMDUCTANCE 


PEAL KWA, IOWA, L)«A,LHK, KWK, IDMK 

DIMENSION KWA(1> ,0CHA(1) ,I0WA(1) ,L WA (1) , I WK ( 1) , POROS (1) 
DIMENSION KWK(l) ,0DWK(1) ,I0WK(1) ,CONDWAtl) ,C0N0WK(1) 

00 10 1=1, N 

CONOWA(I) = KWA(I) * PI ♦ (ODWA(I)»»2 - IQ WA 1 1) *• 2) / (4 . »L WA (I ) ) 
13 CONOWK(I) = tl.-PORCS(I) ) * KWK(I) ♦ PI * (ODWK(I)*»2 - 

• rOWK(I) *«2)/ (4.* LWK(D) 

RE TUP N 

END 

« **«*•*»««««•**«« ««****««*«**«* «« V 

SUBROUTINE RSCCND (NCG , PI ,L WA ,KWA , OOWA , IDWA ,L FT, KFT, OOP , TOR, 

• LR, KR,0DFT,IDFT,LWK, POROS, KWK, ODWK, TOWK, 

• lOFTWtlOPW.GFTTnT) 


THIS SUBROUTINE COMPUTES THE TOTAL AXIAL CONDUCTANCE FROM THE 
CONOENSER TO THE RESERVOIR OF A VCHP. THE SAME ASSUMPTIONS 
USED IN SLBRCUTTNE AXCCNO APPLY. THE WICKS USED IN THE CON- 
DENSER, FEEDTUBE, AND RESERVCIR ARE ASSUHED TO BE lOENTICAL. 
NCG = FLAG FOR VCHP 

LWA = CONDENSER WALL NODAL LENGTH 

KHA = CONDENSER WALL THERMAL CONDUCTIVITY ‘ 

OOHA,IOWA = CONDENSER WALL OUTER AND INNER DIAMETERS 

LFT = FEED TUBE NODAL LENGTH 

KFT = feed TUBF THERMAL CONDUCTTVITY 

ODR.IOR = RESERVOIR CUTER AND INNER OIAmETE°S 


LR = RESERVOIR LENGTH i 

KR = RESERVOIR THERMAL CONDUCTIVITY 

00FT,10FT = FEED TUBE OUTER AND INNER DIAMETERS 

LWK = CONDENSER WICK NODAL LENGTH 

FCROS = WICK POROSITY 

KWK = WICK THERMAL CONDUCTIVITY 

OOWK,IOWK = CONDENSER HICK OUTER AND INKER DIAMETERS 

lOFTW = FEED TUBE WICK INNER DIAMETER (ODFTH = IDFT) 

lORW = RESERVOIR WICK INNER DIAMETER (OORH = IDR) 

GFTTOT = TOTAL AXIAL FEED TUBE CONDUCTANCE 


REAL NCG, LWA, KWA, LFT, KFT, IOWA, TDFT, IDR, LR , KR 
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I 




ii 

i,; 


REAL 

LWK 

, KWK 

, IGHK 

, lOFTW, lORW 


IF (NOG . 

£0. 

u 

• ) 

RETURN 


GWA 


PI 

•9 

KWA ♦ 

(00WA**2 - IDWA**2>/(2 

LWA) 

GFT 

■;= 

PI 

¥ 

KPT ♦ 

(OOFT»»2 - TDFT-**2)/CA 

.* LFT) 

GR 

= 

PI 


KP » 

(OOS**2 - IDR**2)/(2.Y 

LF) 

GWK 


PI 

¥ 

KWK ♦ 

(1.- POROS) ♦ (0DWK*»2 - 

lOWK^'^a)/ (4.^ LWK^ 

GFTW 


PI 

¥ 

KWK » 

U.- P0R0S)*(I0FT**2 - 

I0FTW^^2)/(8.» LFT) 

GRW 

= 

PI 

¥ 

KWK » 

(1,- P0B0S)*tI0R»*2 - 

I0RW»»2)/'(A.* LR) 


GFTTOT = l./U./GHA ^ l./GFT <• l./GP) + 

♦ l./(l./GWK ♦ l./GFTW <■ l./GRW) 

RETURN 

ENC 

********** 

SUBROUTINE A6C0K8{NE,NA,NC,GAE,GAA,GflC,GWAE,GWAA,GWAC, 

* GHKE.GWKA,GWKC,NCTRAP,GFIN,GTCT) 


THIS SURROUTINF. CONPllT'-S THE TOTAL AXIAL CONOUCTANCE OF THE 
ARTERY, L lOUIO-SATURATEO WICK, ANO MALL. THF RESULTS OF 
SUBROUTINES ARTPRY, AXCONO, ANO FADFIN ARE REQUIRED, THUS 
THOSE ROUTINES MUST «£ CALLED BEFORE AGCOH" IN VARRLl. 

N = NLIH8ER OF NODES IN EACH PIPE SECTION 

r-A = ARTERY rONDUOTANCF IN EACH SECTION 

GWA = WALL CONDUCTANCE IN EACH SECTION 

GWK = WICK/LIOUIO CONDUCTANCE IN EACH SECTION 

NCTRAP = COLD TRAP FLAG 

GFIN = AXIAL FIN CONDUCTANCE 

GTCT = TOTAL AXIAL CONDUCTANCE 

E,A,C = °IPE SECTION!? EVA®, ACIA, OR rCN.D 


01 N't NS ION GAE (1) , GAA (i» ,GAC( U ,GWAE tlJ , GM A A ( 1 ) , GH AC 1 1) 
niKENSICN GWKF(l) ,GWKA(1) , GWKC (i) , NCTRAP ( 1) , GFIN ( IT , GTOT (1) 
00 I=i,NE 

IFLNE ,EC. IJ GO TO 20 

10 GTOT(I) = 2./(1./GAE(1U1./GAE(In-1)) * 2 . / 1 1 . /GW AE L U 

♦ H./GKAF <lf 1)) «• 2./(l./GMKE(T)+i./GWKF(I+l)) 

2J GTOTtn = 2./tl./GA£(I) + l./GAA(lH ♦ 2 ./ (1 . / GMAE ( Tl 

'» ti./GMAA (1) ) 2,/(l./GHKE<I)+l./GHKAtl)) 

II = NE 4- 1 

III = NE t NA 

no 31 I=TI,III 

J = I - NF 

IFCIII .£0. T) GO TO kC 

30 GTOT(I) = 2./(l./GAA(J>H,/GAA(J + U) + 2 . / (1 . /G F A A ( J) 

* fl./GHAA(J+l) ) f 2,/(l./GWKA(J)4-l./GHKAtJ4-l)) 

AO GTCTLU = 2./(l./GAA(jH-l./GAC(in + 2./ (i./GWAA ( J) 

» fl./GWAC(l) ) 4- ?./(l,/GWKA(J) 4-1./GWKC (in 


II = N E 4- NA 4- 1 

III = NE 4- NA + NIC 
DO 5C I=ri,III 

J = I - NE - NA 

IFdll .FC. I) RETURN 

GTCT(I) •= 2./ (l./GAC( J)4-1./GAC(J4-1) ) 4- 2./(l./GWAC(J> 

* N 1 ./GHAC (Jfin P 2./(l./GWKC( })P1./GWKC(J4-1») 

IF (NCTRAO ( Jl .EC. i) GTOT(I> = GTOT(I) p 2./ 

* (l./GFiN(J) P1./GFIN(JP1) ) 

50 CONTINUE 

RETURN 


END 


* ******* I ***** < 


SURROUTINS RACCNO(N,°I,EFFKWK,LWK, OOWKjIDMK, KKL,LWL, 
GDHL ,IOWL,CCNO) 


THIS SUBROUTINE CCMPUTFS THE TOTAL RADIAL CONDUCTANCE OF THE 
LIOUIO-SArURATPO WICK AND THE OTRE WALL. THE EFFECTTVP 
THERMAL CCNOUCTIYITY CF THE L IQUID-SATURAT£n WICK P^CM SUB- 
ROUTINE EFPKWK IS RfcvIUIoEO, THUS EFFKWK MUST P£ CALLED BE- 


/ 


;7-.rr 3£,TlVS<<£p' 




bJL~ 
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FORt RACOND IN VARBll. 

N = NUM9HR OF' CASES TO 9E TP.EATEO TN THIS CALL 

EFFKHK = LIOUin/WICK EFFECTIVE THEOHAL COMniJCTTVITY 

LMK = HICK NODAL LENGTH 

OOHKjIOWK = HICK OUTER AND INNER DIAMETER': 

KML = HALL THERMAL CONDUCTIVITY 

LHL = HALL NODAL LENGTH 

ODWL»IOWL = WALL OUTER AND INNER DIAMETERS 

CONO = TOTAL RADIAL CONDUCTANCE 


REAL LHK,IDHKtKKL,LHL,IOHL 

niKENSION EFFKHKIl) ,LHK(1» ,ODHK(l) ,IDHK(1) ,KHL(1) ,LHL(1) 
OIKENSION ODWLtl) ,IOHL (1) ,CONOlll 
DO lu 1=1, N 

CONDWK = 2.* PI * EFFKHK(T) * LHK ( T) / ALOG ( OD WK t I )< IDHK ( I) ) 
CONDWL = 2.* FI * KWL(I) » L HU I ) / ALOG <ODWL < I) / IOHL ( H ) 

10 CONO(I)= 1./(1./CGN0HK l./CONOWL) 

RETURN 

END 

SUBROUTINE '’SLEAK (NCG, KINS, PI, OOP, KINS, LR, GLEAK) 


THIS SU3RCUTINF COMPUTES THE RESERVOIR-TC-ST NK HEAT LEAK CCNPUrTANCES 


NCG 


FLAG FOP VCHP 

KINS 


INSULATION 

THER^*AL CONDUCTIVITY 

ODR 


RESERVOIR 

GUI EH DIAMETER 

XINS 


INSULATION 

THICKNESS 

Lf? 


RESERVOIR 

LENGTH 

GLEAK 


RESERVOIR 

HEAT LEAK CONCUCTANCE 


REAL NCG, KINS, LR 

IF(NCG.£Q.C. .OR. XINS.LE.Q.) RETURN 

GLEAK = 2.* FT * KINS * LR / ALOG UCOR +2.* XINS1/P0R) 

RETURN 

ENC 

SUBROUTINE RAOFI N (N, TR,TRL,TRR,CCNO , X ,T, 7, SI G, ASLCPE .EMTSS, TEX, 
* TENV ,OABSA,GSINK, AXFINC-,TSINK,QNOOE»aTOT) 


•> 


THIS SUPROUTINE COMPUTES THF EFFECTIVE STNK TEMcERATURE (BCUNOAPY 
NOC£» FOR A SPACE RADIATOR. ACCOUNTING FOR FTN GEOMETRY, ABSORBED 
EXTERNAL ENERGY, AND CCNCUCTICN TO ADJACENT RADIATOR PANELS. IT 
CONSIDERS END EFFECTS--FIN RCOT TEMPERATURES MUST BE INPUT SE- 


QUENTIALLY. 

N 

TR 

TRL 

TRR 

CONO 

X 

Y 

7. 

SIG 
A SLOP: 


EMISS 

TEX 


RAOFIN SHCULtJ BE CALLED FROM VARRLl. 

= NUMBER OF CASES TC BE TREATED IN THIS CALL 
= FIN ROOT temperature 

= LEFT SIDE adjacent FIN ROOT TEMPERATURE 
= RIGHT siot ADJACENT FIN ROOT TEMPERATURE 
= FIN THERMAL CONDUCTIVITY 
= FIN LENGTH (HALF PIPE-TO-PIRE DISTANCE) 

= FIN THICKNESS 

= FIN HIOTH <CCNOENSER NODAL LENGTH) 

= STERHAN-BCLTZMANN CONSTANT (.1713E-8) 

: = ARRAY OF NCNDIMcNSIQMAL FIN TEMPEFATURE 

GRADIENTS AS A FUNCTION CF ENVIRONMENT, ADJA- 
CENT PANEL TEMPERATURE, ANO FIN GEOMETRY 
= panel EMISSIVITY 

= EXTERNAL SURROUNDINGS TEMPERATURE 


TENV 

QABSA 

GSINK 

AXFING 

TSINK 

ONODE 

OTCT 


ENVIRCNHENT TEMPERATURE 
ABSORBED EXTERNAL FLU > 

RADIATION CONDUCTOR FROM RCOT TO EFFECTIVE SINK 

AXIAL FIN CONOUCTANCE 

EFFECTIVE SINK TEMPERATURE 

HEAT rejected FROM EACH HEAT PIR£ NODE 

total heat rejected from all nodes in THIS CALL 


nit'EKSIOK TP(1) ,TRL (1» ,T«R(1) ,CONO ( 1 ) , X ( 1) , Y (1) , 7 (H ,EMISS(H 
DIMENSION TEX(l) ,TENV(1» , QAD SA (1 ) , GSINK( 1) ,AXFING(1» ,TSINK(1) 


DIMENSION ONODE(l) 

OTOT = 0. 

DO li, I=lfN 

IF(X(II .ED. J.J GO TO 1C 

C.. .CONVERT TEMPERATURES TO RANKINE 
TEXT = TEX(I) ♦ 46u. 

TRQOT = TR(I) ♦ 460. 

TROOTL = TRL<I) + 46-j. 

TPOCTR = TRRtI) + 46a. 

C. ..COMPUTE FIN PRCFllE NUMPER 

PRCFNO = SIG * EMISS(I» * X(I>**2 * TP00T*-»3 / (CONDdI * Y(I>) 
TFIPROFNO .LT. 0 .OR. PPOFNC .GT.4.) WRTT£{6,1) PROFNO 

1 FORMAT(*0 «-<-*»-«-£RROR»** + <- PRCFNO (*E12. 5*) EXCEEDS ROUNDS*) 

C.. .COMPUTE ADJACENT fiCCT TEMFEPATIIRE RATIOS 

THETAL = TRCCTL/trcOT 

IFdHETAL .LT. .6 .OR. THETAL .GT. 1.5) WPITE(6,2) thETAL 

2 FORMAT(*0+ ++++ERROR++t+'i- THETAL (*£12.5*) EXCEEDS ROUNDS*) 
THETAR = TROOTR/TROOT 

IFtTHETAfi .LT. .6 .OR. THETAR .GT. 1.5) WRITE(6,3) THETAR 

3 FORMAT(*j+* + +*-ERPOR*«- +++ THETAR (*E12.5*) EXCEEDS ROUNDS*) 

C. ..COMPUTE EFFECTIVE ENVIRONMENT TEMPERATURE AND TENVTR RATIO 

TENVIR = ((.P2fl * TEXT**4 + .122 ♦ TR00T**4 
* + OARSA(I)/SIG)/EMISS(I))**.25 

THETAE = TENVIR/TROOT 

IFdHETAE .LT. 0. .OR. THETAE .Gt, 2,3) WRITE(6,4) theTAE 
H FORHAT(* ;tt + *«-FRROF++**+ THETAE (*F12.5*) EXCEEDS ROUNDS*) 


C...LOOK UP FTN GRADIENT AT ROOT 

CALL P30EG1 (THETAL, FRCFNO, THETAE. ASLOPE, DTOXL) 

CALL 030EGl(TH£TAR,Ph0FNC, THETAE, ASLOPE, QTOXR) 

C.. .COMPUTE HEAT REJECTION FROM EACH SIDE CF ROCT, AND TOTAL 
OL = -COND(I) * 7(1) * V(I) * TROOT * OTOXL / X(I) 

OR = -CCNO(I) * 7(1) * Y(T) » TROOT * pTOXR / X(T) 

ONOOE (I) = OL * OR 

C...CCMPUTF SINK TEMPERATURE REQUIRED TO GIVE ONODE 
TSINK(l)- QNOOETD/GSINKd) + TR(I) 

C. ..CONVERT ENVIRONMENT TEMPERATURE TO FAHRENHEIT 
T£NV(I)= TENVIR - 46C. 

C...SUM HEAT REJECTION FOR EACH NODE INTO TOTAL PIPE HEAT ppJECTION 
OTOT = QTOT f GNCO-Ed) 

C. ..COMPUTE FIN AXIAL CONDUCTANCE 

AXFINGd) = CCNOd) * 2.* Xd) * Yd) / 

10 CONTINUE 
RETURN 
END 


lao 


SUBPOUTINr 

niMENSION 

NT 

N1 

N2 

V (1) 

V (2) 

Mi 


03DEG1 (X, y, r. A, W) 

A(i) , V<5) 

= A Cl) . C:^ . G 
= A(2) ,OR. G 
“ A ( 3) 

*= i+.OR.: 

*= A(4) 

= 2 

IF(Z.GT,A(4>) GO TO 2GD 

XI *= A(Ml+i) 

A(Hi+l) = (Nl^N2¥NlfM2fi) ,OR, J 

X2 = A 01 + 2) 

A (Ml+2) ~ ,N1.0P* J 

CALI D20EGi (X, Y, A(^l+i), W) 

A CHl + 1) ^ XI 
ACMl+2) = X2 


FETUPN 


2Cu ^2 ~ M14 Ni»N2f Mf N2 + 3 

IF(M2.GT • NT) GO TO ICQ 
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GO TO 4u0 


N3 = A(K2>.OR.O 

M4 = AtM2+i)*OP.O 

V(4) = A(M2+2) 

IF (2.NE.A 1M242) ) 

N1 - N3 

N2 -Hh 

Ml ♦= M2 

GO TO 100 

400 IF (7.GT,A(M242) ) GO TO 5j0 

XI = ACMl+l) 

AOi + 1) = (M4N2*Nl4N24l).0P.a 

X2 = AIMli2) 

A(Ml+2) = Ni.OR.a 

X3 = A<M2+1) 

A<M2+1) = N34N441I .OF.O 

X4 r A(M2 + 2) 

A<M2>2) ^ N3.0R.G 

CALL 02DEG1 (X, Y, A(Ml4l>, Vt3)l 

CALL 02DEG1 (X, A(M24l>, V(5)) 

CALL OlOEGl (Z, Vtl), Kl 
A(M141) = XI 
A <Mlf2» X2 
A (M2+1) = X3 
A «M2 + 2) = Xa 
RETURN 

5:0 Ml = M2 

N1 - N3 

N2 r N4 

V(2) = V(4) 

GO TO 2u0 

END 

SUBROUTINE FL FRNT C P VA , PVR , Rf TR ,VF # VC ,N ,K , P VS , TS , NN ,L At C ) 


THIS SURRCUTINE COMPUTES THE LOCATION OF THE INTERFACE BETWEEN 
WORKING VAPOR AND NCNCGNOENS IRLE GAS IN A VARIABLE CCNOUCTANCE 
HEAT PIPE. THE PRINCIPLES OF CONSERVATION OF MASS AND THE PER- 
FECT GAS LAWS ARE EMPLOYED# PRESSURES ARE DERIVED FROM TEMPERA- 
TURES ACCORDING TO THE CLAPEYRON EQUATION, DIFFUSION EFFECTS 
ARE IGNOPEO. WHILE AXIAL CONCUCTION EFFECTS ARE CQNSTOEREO, THE 
RATIO OF ACTIVE CONDENSER LENGTH TO TOTAL CONDENSER LENGTH IS 
USED TO MODIFY THE VAPCR-TO-CONDENSER COUPLINGS TO SIMULATE THE 
SHUT-OFF PORTION OF THE CONDENSER. FLFRNT IS CALLED «Y FRONT 
INTERNALLY. 

PVA 

PVR 

P 


= VAPOR PRESSURE IN THE ACTIVE PORTION OF THE PIPE 
= VAPOR PARTIAL PRESSURE IN THE RESERVOTR 
= UNIVERSAL GAS CONSTANT (1S45) 


TR 

VP 

vr 

N 

K 

PVS 

TS 

NN 

LALC 


= PESERVOIR TEMPERATURE 
= RESERVOIR VOLUME 
= TOTAL CONDENSER VOLUME 

= NUMBER OF MOLES OF MONCONDENSIHLE GAS IN PIPE 
= NUMBER OF INACTIVE NODES (NOT NECESSA'^ILY AN INTEGER) 

= VAPOR PARTIAL PRESSURE IN EACH INACTIVE CCNDENSEP NODE 
= TEMPFRATURe IN EACH INACTIVE CONDENSER NODE 
TOTAL NUMBER OF CONDENSER NODES 
= RATIO OF ACTIVE TO TOTAL CONDENSER LENGTHS 


REAL N, K, LALC, KNEW 
DIMENSION PVS(l) ,TSU) 

IF<N ,EO. 0 GO TO 3G 

KNEW = K + .G2 
...ITERATE ON FRONT LCCATION 
DO 25 J=l,15 

IF(APS(K-KN£H) .LT. .01) RETURN 
K = KNEW 

— K = K 4 l.E-13 
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TRR = TR + *163. 

Kt^ = INT(K’> 

C...SUM GAS DENSITY CVcF ALL INACTIVE NODES 
SUY = 0. 

IF(KK .EC. J> GO TO 2C 

DO Iw 1=1, KK 

TSR = TS<NN+1-I) ♦ A63. 

C....NOOE TOTAILY 81ANKETED 

SUM = SUM ♦ (PVA - PVS(NN4l-in /(R • TSR) 

10 CONTINUE 

C....NOriE PARTIALLY '^LANKETEO 

20 SU*^ = SUM ♦ (K-FLOAT(KK)) *<PVA - PV S(N'N-KR ) ) / (R * 

* (TS(NN-KK) ♦•+60.) ) 

IF(ABSISUM) .LT. l.E-5u) GO TO 30 

C...COMPUTF FRONT LOCATION 

LALC = !.♦ ((((PVA - °VR)/(R * TRR)) * (VR/VC) - N/(lA4.*VO) 

* ♦ K / SUM) 

C... ADJUST TO BETWEEN 0 AND 1 

IF(LALC .GT. 1.) LALC = 1. 

IF (LALC .LT. 0.) LALC = C. 

KMc« = (l.-LALC) • FLnAT(NN) 

25 CCNTINUc 

...WRITE MESSAGE IF CONVERGENCE IS NOT SATISFIED 
DELTA = R - KNEW 
— WRIT£{6,A0) DELTA 

— -40 FCPMAT(*0^++4+ CAUTICN ♦♦♦+♦ NO LA/LC CONVERGENCE, DELTA =* , 

— • F9.4,* ACTIVE NODES*) 

RETURN 

30 LALC = 1. 

RETURN 
END 

SUPROUTINE AOJG(NCG,NC,LC,LACT,GCONC) 


THIS SURRCUTINF ADJUSTS IHF CONDENSER COUPLINGS OF A VC HP 
TO SIMULATE THE LOCATION OF THE VAPCR/GAS INTERFACE. 

NCG = FLAG FOR VCHP 

NC = NO OF CONCENSER NODES 

LC = CCNCENSEP NODAL LENGTH 

LACT = ACTIVE CONDENSER LENGTH 

GCCND = FIRST CONnENSER WAl L-TC- VAPOR COUPLING 



REAL NCG, LC, LACT 
DINENSIOK LC(1) ,GCONOf 1) 



TFtNCG .£C« U.) 

CL = G. 

DO 3C I=1;NC 
CL ^ CL 4- Lca) 

RETURN 


IF(CL .LE, LACT! 

GO TC 30 

10 

IF CCL-LCtI) .GE. LACT) 

GO TO 20 


GCOND(I) = (<LACT - (CL-LC(I))) 
GO TO 3G 

/ LC(I)) ^ GCONOd) 

2C 

GCON C(I) = u . 


3J 

CONTINUE 

RETURN 

END 



*»•»<,»**»»♦ JMMi •»¥•■** IMf »++»*»* »♦¥ »* »**■*****¥¥»»**»*»»•* I, 

SUBROUTINE LI ^ITS ( A ,m , G, T V.TE ,PI , GC .H , OR, SH, APRSS 1, APRSS2) 


THIS SUBRCUIINE COMPUTES THE MAXIMUM HEAT LCAO CAPABILITY 
OF A HEAT PIPE. THIS .“AXIMUY LOAD IS LIMITED RY THE CAPIL- 
LARY return mechanism (WICKING LIMIT) . CHOKED VAPOR FLOW 
(SCMIC LIMIT), THE INTERACTICN BETWEEN VAPOR AND LIQUID AT 
THE WICK SURFACE (ENTRAINMENT LIMIT), AND THE TENDFmcy TO 
FORK BURBLES IN THE EVAPORATOR WICK (°CILING LIMIT) . A 
FIFTH LI'^IT CCCURS WHEN AN ARTcRV IS INTROGUCEQ, SINCE THE 
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ABILITY OF CAPILLAFY FORCES TO FILL THE ARTERY IS A FUNCTION 
OF APPLIEH HEAT LOAD (ARTERY SELF-PRIMING LIMIT), SUBROUTINE 
LIMITS CALLS CN SUBROUTINE PERHPL. IF THE PIPE HEAT LOAD 
EXCFEns THE LOWEST LIMIT, THE PERGFNT OVERLOAD IS PPINTEO. 

IF NO LIMITS ARE EXCEEOEO, THE PERCENT UNDERLOAD T5 PRINTED. 

ALSO OUTPUT ARE VAPOR AND LIGUID PRESSURE DROPS, ABSOLUTE PPES- 
SUftES, AND MASS FLOWRATE. LIMITS SHOULD BE CALLED FROM VARBL2 
FOR TRANSIENT ANALYSES, EXECUTION FOR STFA OY -ST ATE* 

A - INTEGER COUNT LOCATION OF FIRST HEAT PIPE ARPAY* 

heat pipe arrays must be input as A PLCCK 

IN THE OROEP SHCWN IN THE ARRAY DATA. 

M = SAME AS A. USED TO ALLOW ACCESS TO INTEGER 

COUNTS AND INTEGER ARRAY DATA VALUES. 

G = FIRST HEAT PIPE EVAFORATOR-TO-V APCR COUPLING. 

HEAT PIPE CONDUCTORS MUST RE INPUT AS A BLOCK 
IN THE ORDER SHOWN IN THE CCNOUrTOR DATA. 

TV = FIRST HEAT PIPE VAPOR TFMPERATURE. 

TE ~ FIRST HEAT PI^E EVAPORATOR WALL TFHPERATUPE. 

HEAT PIPE NODES MUST BE INPUT AS A BLOCK 
IN THE ORCEP SHOWN IN THE NODE DATA. 

GC ^ CONSTANT RELATING MASS AND FORCE (GSUBC) 

H = MECHANICAL EDUIVALEM OF HEAT 

OR = OEGREES/RAOIAN 

SH = NO SECCNDS^'^a PEP HOUP»»2 

AFRESS ^ SATURATION PRESSURE VS TENPERATUop ARRAY 


OIHENSICN 

A(l) 

, H ( 11 

,G(1),TV(1),T£(1) 

PEAL MUL, 

MUV, 

KL, 

KV, MOOT 

N 

M(l) 



Ni 

N2 = 

N3 = 

NN ^ 1 

li 

12 = 

13 = 

1 

00 1000 I 

= 1,N 



II 

M(l) 

4 1 


J 

M(II4 

1) 4 

1 


K = 

L 

JJ 

KK 

LL 

EVAPL = 
AOIAl 
CONOL = 
QLAT 

ACTL = 

DENL 

OENV 

CPL 

CPV 

HUL 

Huy 

KL 

KV 

HFG 

HFS 

IFd.LE. 

TFd.GT. 

A OB'^IIfH 
C.*. .WIG KING I IM 
IFLAG = 
DPCAP = 


M(3m + 6»J + 1) 4- 1 

MC5^II^6»J4-6*K + 1) + 1 
M(26^Il434»Ji-33»Kf42»L4-i) 4 i 
M(27»Il434*J433^Kf42^L4jJ4l) f i 
M(284II+34'' J + 33»K442»L4JJ4KK41) 4- 1 

ACTI+6»J+1+I) 

A (3^1l4e^J46^K4-14I) 

A (5»TT464J + 6»K4-6^L414T) 

A(33^Il447^J446*K4-55»L4Jj>KK4LL4i^I) = G<NN) ^ 
(TEtNl) - TV(in 


ODCAPI = 
OP BOO = 


A(464II+49^J44d4Kf5 
A(2S*II^47»J439^K44 
A(28*II447M446'^K45 
A1284II + 47^34^u^K4-4 
A (29^II447»J445^K45 
A(28»II447»J4 J+14K+4 
A(2G^II447»J4464K45 
A(28»Il447434ii2»K4tf 
A(3l^Il447^j4-46^K4-5 
fl(E8^Il4ii7»J443^K44 
A(28»II447^J4444K44 

A(284IIf.474j4>^cA^K + 4 
6) CALL 

6) CALL 

7^j4He<X 45 5^L4 J J + KK 
IT 
D 

2.* SURFT * Cn''(A(6 
fl t6''IT + 15*J+£»><+6* 
DPCAP / IhA. 


9*L<-JJ+KK+LU+1 + I) • CCKCL 
8*L*-JJ+KK + LL + 1+N2> 

5»L+ JJ<-KK«-LL + 1+T > 
8*t«-JJ+KKtLl + l + N2) 

5*L i-JJ+KK^LL + l+I) 
8*L+JJtKK+LL'm-N2) 

5*L*-J J+KKtLL + l+I) 

8*L + JJ+KK+LU-H-N?) 
5»L+.JJ4-KK+LL+1 + I) 
B*U-JJ+KK»LL + 1+N2) 

6»Lf JJ4-KK4-LL+H-N2) 
e*L+JJ+KKftt+l+N2) 
OiO£Gl(T£(NN) , APRSSlfPEVAO) 
OlOEGl (TE (NN) , AFRS92f P£VAo) 
+ LL+1+K1) = P£\/AF 


*II*-2i»Jt2i*'<»-21*L+l<-I)/0P) / 
L+l+Nl) 


OEM » A(6*II + 2i*J4-21'*-K + 21*L«-l+I) ♦ ( A ( 10*1 14-21»J 
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c 


c 


c 

c 


c 


c 


* * (EUflPL + AOIAL+flCTL) * CCS(A(g*TI 

• +21*J«-2i*K«-21»L+l+I)/D°) + {A{II+i**J^l + Nl) +A(3*II 

♦ <-6* J+(t*K + l+N2) tA (5*II»6*Jf6*Kf4*L+H-W3) ) ^3.3 

* ■ » SIN(A(q + rH-21*J«-21*IO-21^L + l + I>/DR)l / GC 
opgooi = DPsr.o / ia4. 


CALI PERf-BL (K{6»ir + 6*Jt6*K+6*L+l + Nl) , A (6»I !♦ 10 * J + 6 *K f 6*L + 1+ N1 ) , 

* A (6*Iim*J + 6*K*6*Lm-Nl) , A( 6*11+1 B*J*6*K<-6*L>1 *N1 ) , 

* A (6* II+19*J+6*K + 6*L+1 + M) , A f 6* 11 + 20 * J + 6* K + 6*L +1 + K' 1 ) i 

* M 11*11+21* J+ 21 *K+ 30 *L+1+N1) , A ( 1 1* I I +22* J+21 *K+ 3Q *L +1+N1) , 

* A (ll*IT + 23*J+21*K + 3-*L +1 + N1) , A (1 1 * 1 1 + 3 2* J + 2 1 **< + 3C *L + 1 + N1) , 

* A (26*II+34*J*33*K+42*L+1+I1) , H (2 5*1 1+3 4* J+33 *k+42 *L +1+ I ) , 

* A (11 *11+27*3+21 *K+30*L+1+Ni) ,A (ll*II+24*J+21*K+30*L+l+Ni) , 

* A (11*1I + 25*J+21*K + 30*L + 1 + N1) , A ( 6*1 I + 16 *J * 6 *K +6*L + 1 +N1) , 

* A (28*11+3 7*3+36 *k+^5*L+JJ +KK+LL+1+N1) , PI ,DR) 

CALL PERt'BL (M (6*II + 21*J+6*k+6*L + l + N2> , A(6*II+21*J + 13 *t< + 6*L + l + N2) , 

* A (6*Il+21*J + ll*K + 6*l + l+N2) ,A ( 6*11+ 21* J+l« *n'+6*L +1 +N2 ) , 

* A (6*II+21*J + 19*J< + 6*L + 1+N2) , A ( 6*1 1+ 21*3 + 2 0* K + 6*L +1 + N2) , 

* M (11*I1 + 34*J + 21*K+30*L+1+N2) ,A (ll*II + 34*3 + 22*K+30*L + l+N2) , 

* A(11*II+34*3+23*K+3C*L+1+N2» , Atll*II+34*3+32*K+3u*L+l+N2) , 

* A {27*II+34*3+33*k+42'»L + J3+l+12) ,H{ 26*II+34*3 + 33*K+42*L+33+1*T» , 

* A (il*II+ 3tt*3+27*K + 3;*L +1 +N2) , A (1 1* 11+34*3+ 24*k+ 33 *L + 1+M2) , 

* A (11 *11+34 *3+25*4 + 3 ;*L +1 + N2) , A (6*ri + 21*3 + 16*K+6*L + l+N2) , 

* A (28*II+38*3+36*K+45*L+3J+KK+LL+l+N2) , PI, OR) 

CALL FERMPL (M { 6* 11 + 21 *3+21*K+ 6*L + 1+N3) , A ( 6* H + 2l* j + 2i »)<+10*L+l + 

* N3»,A(6*II + 21*3 + 21*K+11*L + 1+N3) , A ( 6*11 + 21*3+21*4+ 19*L+1 +N3) , 

* A (6*1I+21*3+21*K+19*L+1+N3) , A (6*TI+2l*3+21*K+20*L ♦ 1+N3) , 

* M (11*11+34* J+33*K+3w*L+l+N3) , A (1 1* 11+34*3+33 *4+31*L+l+»3) , 

* A (ll*TI+34*3+3 3*K+3 2*L +1 + N3) ,A(11 *11 + 3 4*3+ 33*4+41*L+l+’*l3) , 

* A (26*11+34 *3 + 33* k+42+L +33 + KK + 1+I3) ,M (2^* TI + 3u*3 + 3 3*4+42*1 + 

■ • 33+4K+1+I) ,A { ll*II+34*3+33*K+36*L +1+M3) , A ( 11* II + 3 4*3+33 *K+33*L + 

* 1+N3t ,A(ll*II+34*3 + 33»4+34*L+l+N3) , A (6 *II + 21 *3+ 21 *K + 16*L + 1+ «•*» , 

* A (2P*II+38*J + 37*K + 4 5*l.+33+KK+LL + i+N3) ,PI,0R) 

SDPLIQ = .5* >(UL * EVAPL / (9P * GC * OEML * HFG * A(26*T.I + 

* 37*3+36*K+45*L+33+4K+LL+l+Nl)) + 

* MUL * AQIAL f (SH * GC * DENL » HFG * A (26*11+ 

* 38*3 + 36*'<+45*L+3J+KK+LL + l+N'2) ) + 

* .5* MUL * ACTL / (SH * GC ♦ DENL * HFG * A(26*II+ 

* 3fl*3+37*K+45*L+33 + 44+LL + l+N3)) 

DPLIOI = SOFLIC * QLAT / 144. 


SOFl/PA = 126.* HU\/ * AOIAL / (PI * SH * GC * HFG * DENV 

* * A ( 28*II+37*J+35*4+44*L+3J+4K+LL+l+N2) **4) 

npVPAl = SQPVPA * QLAT / 144. 

150 IF(npBOO .G£. CPCAP) GOTO 230 

OOWL ^0. 

R£I = 40. 

RES = 40. 

DO 220 MMM =1,50 

SDPVPE = 6.* MUV * EVAPL / (FI * SH • GC * HFG * OFNV 

* * A (26*TI+36*J+35*K+44*L+33+4K+LL+l+Nl) **4 

* * (.0481+.04g4/(4,7 + P.EI)**.8)) 

CPVPEI = SDPVPE * QLAT 7 144. 


TF(A9S(RES) .LT. 4.5975) GO TO 20 3 

SDPVPC = 32.7(P1**2 • SH * GC • OENV * HFG**2 * A(26*II+ 
* 37*3 + 36*4 + 44*1 +33+4K+L L + 1+N3) **4) 

OPVFCl = SOPVFC * CLAT**2 / 144. 


OWL = (-CSOPLTQ+SOPVPA+SOFVFE) + SORT ( ( SOPL Tfi 

• +SOPVPA+SOFOFt ) **2 - 4.*snPVPO*{OPPCO- 

* OOCAP) ) )/ (2.*SDPVPr) 


c 


GO TO 210 
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2C3 SOFVPC = 8.* MUV * ACTL / IPI ♦ SH • GC ♦ HFG • OENV 

* • fl(28*II+37*J+36»Kt44*L4.JJ»KK4-LL*l+N3) **4 

* » (.3481*.u4?4/<4.7-Al=>?UES))**,8)) 

DPVPCI = SOPVPC • GLAT / 1^44. 


QWL 


= (oocAP-cP8now(SDPL iQ4-snpVPE+snPVPA+snPVPC) 


210 OQWL = APS(QWL-CQWL) 
IF(OQWL .LT. .05) 


GO TO 26C 


¥ 

¥ 


OPCAP -*f 
OPLin 
OPVPE =*r 


OQWL ~ CWl 

REI = A (40*II + 4Q»J«-47*K+56*L+JJ«-KK>LL*'1 + I) = OWt / 

» (PI * HFG • cVAPt ^ MUV) 

RES = AC^l^XIf^a^J + ^r^K+Se^L+JJ + KK + LL^H^-I) = OWL / 

♦ CPI * HFG ♦ ACTL * MUV) 

220 CONTINUE 

GO TO 24u 
233 QWL = 0. 

GO TO 263 

2 WRITEC6,250) I, DOWL, IFLAG 

250 FORMATC^C 4+f+4CAUTION+++++ HCRE THAN 50 ITERATIONS REQUIRED ^ 

♦ ^TO COMPUTE QWL FOR PIPE‘S, 12, » OQWL = ♦fFlO.I, 

» * 9TU/HR, IFLAG =^,I2) 

260 IFCIFLAG .£0. 1) GO TO 261 

QWLIM = A(34^II + 47»J4'H6^K455*L4-JJfKK+LL +1 + 1) = GWL 
GO TO 262 

261 OPRIH = QHL 
GO TO 266 

262 WRITE(6,265) I, DPCAPI, DPBOTI^ DPLIQI, QPVPAI, OPVPEI, DPVPCI 

265 FORMAT(+0 PRESSURE DRCFS FOP */+ 

E12.4,^ PSIO»/^ DP800 =+,£12.4,+ PSIO+/+ 

E12.4,+ PSIO+/+ DPVPA =+,£12.4¥+ PSID+/+ 

E12.4,+ PSIO+/+ DPVPC =+,£12.4,+ PSID+) 

C... ARTERY SELF-PRIMING LIMIT 

TF (M(li+TT+21»J+21 + k:+3w + L + l+M) .£Q. 0 .AMD. 

+ M(11+II+T4+ J+2i+K+3b+L+l+N2) .£Q. fl .AND. 

+ M(ll*II+34+J+33+K+3u+L+l+N3) .EQ. 01 GO TO 266 

OPCAP = 2.+ SURFT + CCS ( A (6+II+21+ J+ 21+K+ 21+L +1+ IT /DR> / 

+ A(ll+lI+33+J>21+K+30+L+l+Ni) 

IFLAG =1 
GO TO 15u 
C... .SONIC LIMIT 

266 OSLIM = AC35+II+47+J+46+K+55+L+JJ+KK+LL+1+I) = DENV ♦ HFG + 

+ PI + A(2ft*II+36+J+35+K+44+L+JJ+KK+LL+l+Nl)++2 + SQRTCA 

+ (7+11+21+ J+21+K+21+L+1+I) + A(3g*II+47+J+46+K+55+L 

♦ +JJ+KK+iL+l+Nl) + SH + GC + 144. /(2.+ DENV + (A 

+ C7+II+21+J+2i+K+21+L+l+T)+l. )1)/4. 

C.,.. ENTRAINMENT LIMIT 

QELEV = PI + A (28+IX+36+J+35+K+44+L+JJ+KK+LL+1+N1) *+2 

♦ + HFG + SQRT(SURFT+SH+GC+CENV/A(6+II+17+J+6+K 

+ +6+L+1 + NDT /4. 

QELAO = PI + A (28+Il+37+J+35+K+44+L+JJ+KK+LL+l+N2)^+2 
+ + HFG + SORT (SURFT +SH+GC+0ENV/A<6+TI + 21 + J + 17+K 

+ +6+L+1+N2)) /4. 

QELCO = PI + A(28+II+37+J+36+K+44+L+JJ+KK+LL+l+N3r ++2 

♦ + HFG + SQRT(SURFT+SH+GC+D£NV/A(6+II+21+J+2i+K 

+ +17+L+1+N3) )/4. 

QELIM = A(36+Il+47+J+46+K+55+L+JJ+KK+LL+l+I) = AMINICQELEV, 

+ n£LAn,0£LCG> 

C... .ROILING LIMIT 

QBLI^ = A(37+II + 47 + J+46»»^+55+L+JJ+KK+LL + i+I) = 2.+ PI + 

» (TElNl)+460.) + SURFT + EVAPL + ( A ( 6+II+iO +J+ 

+ 6+X46+L+l+Nll+A(6+II+ll + J + 6+K+6+L+l+Nin ♦ A(42 + II 

+ +48+J+47+X+56+L+JJ+KK+LL+l+Ni) / (H + HFG + DENV + 

+ A(6+II+i2+J+6+K+6+L+l+Nl) + (A(6+lX+iO+J+6*K+6+L+l 

+ +Nl)-AC6+ri + il + J+6+X+6+L + l + NUn 

WRITE (6, 2 70) I, GWLIM, OSLIM, OEtIM, ORLIM^ 
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270 FORMAT(*C LIMITS FCR P1P£*,I2,* --*/* PICKING LI«TT =», 

» F15.1,* BTU/HR*/* SCNIC LIMIT =*,F15.T,* PTU/HR*/ 

* * ENTRAINMENT LIMIT ='',F15.3,* inTU/HR*/* ROILING ♦, 

* ►LIMIT =»,F15.3,* RTU/HR*) 

IFdFLAG .EC. 1) WRITE(6,275) GPRIM 

275 FORMAK* ARTERY SEL F-* /7X, ► PRIMING LIMIT =*,F15.3, 

* ► RTU/HR*) 

C... COMPUTE MIN LIMIT AMI OPERATING POINT CF PIPE 

QMLIM = AL3e*II + A7*Jt46*K+55^L+ JJ + KK*LL+1*I) = AMINKQHLIM, 

* QSL IM,Q£LIM,OBLIM) 

IFTOMLIM .EQ. C) OMLIM = 1 . E- 1C 

PCTLIM = QLAT / CMLIM ► 10 0. 

IF(0ML1M .GE. QLfT) WRITE(6,2MU I, °CTLIM 

IFLQMLIM .LT. QLAT> WRITE(6,290) I, PCTLJM 

260 FORHATt»C PIPE*, 12,* OPERATING AT*,F12.3,* PERCENT OF LIMIT*) 

293 FORMAT(*0 PIPE*. 12,* LCAO ARCWE CAPACITY 9Y*F15.3* PERCENT*) 

C... COMPUTE MASS FLOWRATE 

MOOT = A (32*ll*47*J + 46*Kt55*L*JJ + KK + LL*14-r) = CLAT / HFG 
WRITE (6,300) I, MOOT 

300 FORMAT(*j PIPE*, 12,* MASS FLCW®ATE =*,£12.4-,* L6M/HR*) 


1000 


N1 = 

N2 

N3 = 

NN = 

11 = 

12 
13 

CONTINUE 


N1 + M(1*I) 

N2 + M (2*II,6*J*1+I ) 

N3 * M<h*II*6*J+6*K*1*I) 

NM <■ Mti + I) f MT2*IH-6* ) + l*I» + m(4*IT + 6*J + 6*K*1*I) 

11 + M (25*lIt3«t*Jf33*K*42*L + l+I) 

12 ♦ M (26*IH-34*J + 33*K*42*L+JJ-H+I) 

13 + M (27*IH-34*J+33*K*42*L*JJtKK + l + I) 


RETURN 


ENC 

FSTOP 


SUP.FOUTINE FRONT 


THIS SU9RCUTINE PROVIDES THE LOGIC REQUIRET TC CCNTROL THE 
ITERATIVE FRCNT LCCATICN SOLLTTCN, OUE TO THE VARIATIONS 
IN SENSITIVITY CF WORKING FLUID, NOMCONDENSI RLE GAS, ANO 
GEOMETRY COMRINATICNS, THIS ROUTINE WILL PROBABLY REQUIRE 
MOCIFICATION TC ACHIEVE PROPER OAMPTNG ANO CONTROL. FRONT 
CALLS FLPRNT INTERNALLY, ANO REQUIRES PLACEMENT IN EVECTN 
FOP STEADY-STATE CASES, ANO IN VARBL2 FCR TRANSIENT CASES. 


FCALL COMMON 

I = I TEST 

N = M, 520+1 

Ml = lOUM 

N2 = ODUM 

N3 = KOUM 

L = MDUM 

LOOKUP{M5GO + I ,T131*N1-1, A1012,A758+N1) 

L OOKUP(M510*I,T1C2+N2-1, AlOl 2 , A759 + N2) 

L 00KUP(M52C + I ,T1C3 + N3-1, AlOl 2 , A763 +N3 ) 

010EGKT119+L-1, A1012,A761+I) 

A632+I = A820-tI 

VC = 0. 

DO 5 11=1, N 

5 VC = VC + A526+N3+II-1 * R1 * A689+N3+II-1**2 / 4. 

A836+I = (A527+T - A820+I) / A526«-N3 

IF(A836+I .GT. FLOAT(N)) A836+I = FLOAT(N) 

IF(A836+I .it. 2.) A836 + I = 3. 

FLFRNT( A759 + N2, A761 + I,R7,Tlig*L-l, A809+T,VC, ABO 8* I, 

* A336+I,A763+N3.T103+N3-1,«520+I,A810+I) 

WRITE (6,98) I, A81J+I 

F 98 FOPmatC*; PIPE *,I2,* LA/LC (PPELIK) = *,Fi0.4) 

A82£fl = A527+T * A81J+I 

DELTA = A820+I - A832+I 
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IF (ABStOcLTA) •IT, ,05 ♦ A525<*N3) GO TO 100 
IF(M0D(M761+I|2) .EO, Z) A841<1= DELTA 

IF(MOO(M78i+I,2) .£Q. 1) A842^? = DELTA 

C— IF(KOO(M7814T,3) .£Q. 21 AflA3^X = CELTA 

IF « ( A841 + I.GT . 0, . AND. A84 2+I.LT.0 .) •OR, 

* (A841+I.LT .C, .AND. Afl^rE^I. GT.O. H 

* A8444-I “ A 844+1 ♦ 2. 

IF(0£LTA .GT. A526+N3) A62u+I = A832+I + A526+N3/A844+I 

IF(D£LTA .LT. -A52E+N3) A820+I = A832+I - A5 26+N3/A844+I 

IFCDELTA.GT.O . .ANI. DEL TA .1 T . A526 + N3) 

* A82ti + I = AB32 + I + A526+N3 / A844 + T / 4. 

IF(OELTA.LT.O . .AND. DELTA .GT .-A52 6+N3 ) 

A820 + I = A832 + I - A526 + N3 / A844+I / h. 

A810+I = A82Q+I / A527+I 
WRITE (6, E9) If A31Q+I 

F 99 FORMAT(» PIPE LA/LC (FINAL) = ^.FIC.A) 

IF(A8S( A82C+I-A832+I) .IT. .01 ^ A526+N3) GO TO 60 

IF(M78i+I .GT. 4U) GO TO 53 

M781+I = (M781+I) + 1 

N783+I = 0 

RETURN 



50 

WRITE (6,51) 



F 

51 

FORMAT (♦ J ♦+4++C AOTTON+ ♦♦++ 

FRONT 

LOCATION NOT SOLVED IN 40 ♦ 

F 


^ITE^ATIGhS^) 





GO TO 10 0 




6Q 

WRITE (6,61) 



F 

61 

FORMAT(*u«*»++NCT£+»++* FRO^T LOCATION CHANGING LESS THAN* 

F 


• i PERCENT CF NOOE 

LENGTH 

EACH ITEPATICN, ACCEPTING*) 


10 0 

H7f2*I =1 




RETURN 
F END 

END 

RCO 3EM0 OF DATA 




GJCJICJ //// END OF LIST //// 
CJCJICJ //// END OF LIST //'/ 
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2 . HEAT PIPE MODELING 


A, GOVERNING EQUATIONS FOR FIXED-CONDUCTANCE N(W-ARTERIAL PIPES 

The heat pipe Is a device which transfers heat nearly Isothermally 
through the vaporization and condensation of a working fluid within 
a closed container. Referring to Figure 2-1, the processes occur- 
ring within the heat pipe are: 1) vaporization of the working 

fluid In the evaporator due to heat Input, 2) flow of the result- 
ing vapor toward the condenser due to the pressure difference 
created by the evaporator/condenser temperature difference, 

3) condensation of the vapor in the condenser due to heat re- 
moval, and 4) flow of the condensate toward the evaporator due 
to the capillary pumping action of the wick material. 


EVAPORATOR CONDENSER 



Each of these processes, as well as the more conventional heat 
transfer by axial and radial conduction, will be treated In this 
section. The limitations on heat pipe performance will also be 
briefly discussed. 
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1 . 


Capacitance 


The capacitances (mass x specific heat) of the basic part of 
the heat pipe are computed in subroutines CAPWCK and CAPFIN. 

The capacitance of the vapor is assumed to be negligible due 
to its' low density and heat capacity. The pipe wall capacitance, 

(5 , is 

P 







( 2 - 1 ) 


The capacitance of the liquid/wick combination matrix dn is 
equal to the sum of the liquid capacitance and the wick material 
capacitance: 


^liquid ^wick 




+ (1 - 0 ) p C 
' w p 

I 


w w 


(2-2) 


where the porosity 0^ or void fraction of the wick is used to 
separate the liquid and wick material components of the liquid/ 
wick matrix. 

The inputs required for subroutine CAPWCK are Pp, Cp , tt , Dq^, 

^ip» > ^p^? ^ 

CAPFIN and ARTERY capacitances, from subroutines CAPFIN and ARTERY 
respectively if applicable. The capacitance of each pipe section is 
output. There are no restrictions to the use of CAPWCK^ The sub- 
routine has sufficient flexibility to allow various wick structures 
and dissimilar wick and wall materials. 

The radiator fin capacitance is obtained by calling subroutine 
CAPFIN. This routine provides an automated means of updating the 
fin node capacitances in a loop by the equation 

^ tf f • 

The calling sequence is 

CAPFIN (N, X, Y, Z, DEN, CP, CAP) 


a 
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[ 

h 

t 


\^ere N 
X 
Y 
Z 

DEN 

CP 

CAP 


number of fin nodes 

fin length (% distance between tubes) 
fin thickness 

fin width (condenser length) 

fin density 

fin specific heat 

fin capacitance of node 


2 . Axial Conductance 

The conductances (thermal conductivity x cross-sectional area/ 
heat flow path length) of the heat pipe in the axial direction 
are computed in subroutines AXC0ND and AGC0MB. The conductance 
of the vapor is assimied very small due to its low thermal con- 
ductivity, The pipe wall axial conductance G ^ iri any node 

. • p j 

IS given by 


r L /2 n 

_ . p: 


r L /2 n 

p 

2 2 

k 7 t(d - d: )/4 
L p ' 0 1 J 

P P 

T- 

k 7 t(D^ - )/4 

L-P ^ o JL -J 

^ p p 

1 


G 

p,ax 


1 


r 2 l ”1 


r 2 l “1 


P 


P 


2 9 

k 7 t(D^ - D7 ) ' 

*r 

2 2 

k 7 t(D^ - d: ) 


1— p 0 1 - 

1 

L-P 0^ 1 J 

0 


(2-3) 
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The axial conductance of the liquid/wick matrix is equal 

to the sum of the liquid conductance and the wick material con- 
ductance* However, the liquid conductance is small compared with 
the wick when the axial distance is large compared with the wick 
thickness. Therefore 



1 

/ 

L_ 

t 

1 — ; 

! 

j 

! 

1 

2 2 

'T 

2 2^ 

k 7t(D - D7 ) (1-0 ) 

W 0« 1 ^ ^ 

2 W 

1 

k 7t(D - DT ) (1-0 ) 

W ^ 0 1 ^ ^ w^ 

L- w W 


(2-4) 


Again the porosity 0w is used to isolate the wick material. For 
symmetrical (screen-type) wicks, contact resistance between axial, 
and circumferential wick components effectively removes the cir- 
cumferential components* contribution to axial heat transfer. 
Therefore, the wick cross-sectional area 



is halved to isolate the axial wick components* contribution. 

For unsymmetrical wicks (grooves, axial wires, etc.), the 
appropriate effective porosity is used. Subroutine AXC0M) is 

used to compute the axial conduction of each node; subroutine AGCOMB 
is used to combine the axial conductances of each node in a node center 
to node center overall conductance. 


The inputs to subroutine AXC0ND 


are 



D. , L and 0 . The wall and v^ick conductances are output. The 
1 w w 

w 

only restriction to the use of AXG0ND is its assumption: small liquid 

conductance compared with the wick material. This subroutine is flexible 
enough to accommodate various wick structures and dissimilar wall and 


wick materials. 




OF IHE 
FOO^ 


REPRODXJCIBILITy 
ORIGINAL RAGE IS 


The total axial conductance of the pipe section is then evaluated 
by subroutine AGCOMB. This routine will treat all pipes including 
those containing arteries. The routine must be called after calls 
to ARTERY, AXCOND, and RADFIN in Variables 1. 

The results of the above calls yield the following conductances 
which must be combined reducing to a single axial conductance 
between node centers . 



The calling sequence is 

AGCOMB ( NE , NA , NC , GAE , GAA , GAC , GWAE , GWAA, GWAC , GWKE , GWKA, GWKC , NCTRAP , 
GFINjGTOT) 

where 


N 

GA 

GWA 

GWK 

NCTRAP 

GFIN 

GTOT 

E,A,C, 


= HUMBER OF NODES IN EACH PIPE SECTION 
= ARTERY CONDUCTANCE IN EACH SECTION 
= WALL CONDUCTANCE IN EACH SECTION 
= WICK/ LIQUID CONDUCTANCE IN EACH SECTION 
= COLD TRAP FLAG 
= AXIAL FIN CONDUCTANCE 
= TOTAL AXIAL CONDUCTANCE 
= PIPE SECTION; EVAP, ADIA, OR COND 


« 
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Radial Conductance 


The conductances of the heat pipe in the radial direction are 
computed in subroutine RAC0ND. The radial conductance of the 
pipe wall, computed as 

27rk L 

P = P P 


p) 


(2-5) 


The radial conductance of the liquid/wick matrix, j., is 

computed by 



(2-6) 


where the effective thermal conductivity of the liquid/wick 

matrix k is computed in subroutine EFEBCWK, using the methods 

of Gorring and Churchill for screen wicks (Ref, 1) as discussed 

in Section" 2-A-6 . G ,r and G n* are combined into 

P i:/w,r 


G = 
r 


p,r f/w,r 

the total radial conductance. 


The inputs to subroutine RAC0ND are tt, kp, Lp, Dopj ^ip> ^i/w’ 

Lw> I^o^j ^iw* conductance is output. No restrictions 

exist for this subroutine, and it can handle dissimilar wall and 
wick materials. 

Flow Effects 


The enthalpy change in the flowing liquid in the heat pipe 
Genth computed by 


G . . = mC , 
enth p’ 


(2-7) 
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where the mass flow rate m = -r^ 

hfg 


(2-8) 


The processes of evaporation and condensation occurring within the 
pipe are associated with small pressure drops. Since the vapor 
and liquid are in equilibrium with each other within the closed 
system, the Claus ius-Clapeyron equation relating temperature and 
pressure at saturation conditions is applicable: 


dp 


T 


V 


Pv hfg H 


(2-9) 


Consequently, once a vapor pressure drop in the pipe is determined, 
a corresponding temperature drop can be computed. This temper- 
ature drop can, in turn, be translated into a "conductance" if 
the heat transfer rate is known. Derivations of the vapor pres- 
sure drops in a heat pipe, and the corresponding "conductances," 
follow. 


For evaporation, conservation of axial momentum requires that 


evap _ 
dz 


2c- p 

f V 

D^(3600)^ g^ 


8 d 

D^^3600)2 g^dz 



P ^ 

V V 


rdr, 


(2-10) 


dp 

or the axial pressure gradient evap is equal to the sum of the 

dz 

friction loss and the momentum gradient (Ref o 2) . Since u^ = 0 
at z = 0 and z = L (there is no velocity at the ends of the pipe) , 
the momentum gradient term is zero. Therefore, 


dP "2c p U 

evap _ f V 

D (3600)^ g 
V c 


or 


AP 


evap 


=/ 


^ 9n 

2c p U 
f ^v 


0 0^(3600) 


dz " 


2c^ P U L 
f V 

D (3600)^ g 
V c 


( 2 - 11 ) 


( 2 - 12 ) 


The skin friction coefficient cj can be computed from the follow- 
ing relation derived from Navier-Stokes with wall suction or 
injection, assuming laminar fully-developed flow (Refo 3); 
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Re / ,0481 + 


(4.7 + Re ) 
w' 


(2-13) 


which reduces to the Hagen-Poiseuille solution c^ = — at Re =0 

f Re w 

(no suction or injection) . The Reynolds numbers Re and Re„ can 
be expressed as 


UD p 

T, V V 

Re = > 


(2-14) 


Re = — 

W TTfl^ wix^L hfg 


(2-15) 


Substituting equations 2-13, 2-14, and 2-15 into 2-12 and 
using the result in equation 2-9, the temperature drop ATevap 
associated with the evaporation process can be determined: 


8 T M LQ 

V V ^ 

tt( 3600)^ g^ H hfg^ pj^ I". 0481 + 


(4.7 + iRe^l ) 


Converting this AT into a "conductance" using G = Q/ AT, 


8 T fx L 

V V 

^(3600)^ g^ H hfg^ pj^ r.0481 + 


.0494 

(4.7 + Re I )' 


This conductance value is valid in the evaporator for all positive 
values of Re^. It is also valid in the condenser when Re^^ > -4,5975. 
At larger negative values of Re^ (higher suction rates) , cf drops 
to zero, and the following derivation is used. 

Again, conservation of momentum requires that equation 2-10 be 
satisfied. In the case of zero cf, however, the friction term is 
neglected. Therefore, 


(3600)^ g^ 


p u rdr, 

^v V 
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or 






AP 


cond 


(3600)^ g^ 



2 

u 


V 


rdr. 


(2-18) 


The velocity profile can be expressed as 

r2 


U = 
V 


1 dp 


u dz 


V 

4 


[i - 

V J 


(3600)" g^, 


(2-19) 


from Navier “Stokes for laminar flow in a circular tube of con- 
stant cross-section (Poiseuille flow) (Ref. 4). Substituting 
equation 2-19 into 2-18 and integrating, 

^■■cond ■ 7^ (3600)2 ^2-20) 

8 /X \ / 

V 

For Poiseuille flow, 


ii 

dz 



8m u 

— — , 
(3600)^ 


(2-21) 


from Appendix A. Substituting equations 2-21 and 2-8 into 
2-20 and the result into 2-9 • 


AT 


32 T Q 

V 


cond 


2 2 3 2 4 

7T (3600) g^ H hfg p^ 


or, again converting to a "conductance,” 


( 2 - 22 ) 


cond 


_1 

32 T Q 

V ^ 

2 2 3 2 4 

TT^ (3600)^ g^ H hfg%^^ 


(2-23) 


This conductance is valid for the condenser when the wall Reynolds 
number Re^ < -4.5975. Since the adiabatic section can act as a 
condenser during transients (especially startup) , the relation- 
ships derived for the condenser are applicable to the adiabatic 
section as well. When the vapor temperature exactly equals the 
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adiabatic section liquid/ wick temperature, condensation will 
cease and the conductances will transfer no heat. 

The remaining vapor pressure drop, and corresponding temperature 
drop, is due to the friction loss of the vapor flowing axially 
in the tube. Thus, 



■2c, 


(3600)' 



g D 
c V 


(2-24) 


For this case, the skin friction coefficient Cf « 16/Re, Sub- 
stituting cf and equation 2-14 and integrating equation 2-24 
over the pipe length gives 


AP 

V 


128 Q 

2 4 

7t( 3600) hfg 


From equation 2-9, 


(2-25) 


AT = 
V 


128 Q p. LT 

^ V V 


tt( 3600)^ g^ H hfg^ 


(2-26) 


or 

G = ^ 1 

V 128 p. LT 

V V 

7r(3600)^ g^ H hfg%^^ 

The resultant network is shown in Figure 2-2.* 


(2-27) 


The typically very large vaporization, condensation and vapor 
flow conductances Ggyap> *^cond> 10^ to 10® Btu/hr-®R)» 

in combination with the typically very small axial conduction 
conductances Gp^a ^,a ( ^ 10"^ Btu/hr-®R), result in 

slow MITAS thermal analyzer solutions. Fortran logic to alleviate 
the problem in steady-state has been developed, but the transient 
case presents more difficulty because of the additional heat 
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storage term and the lack of an energy balance criterion. The 
more logical solution to the problem, is to eliminate the large 
conductances, making the vapor and liquid in the pipe isothermal. 
This results in a loss in the capability to conveniently handle 
internal phenomena, such as liquid freezing or wick dryout. How- 
ever, since these cases can be treated with suitable Fortran logic, 
the loss is minimal in most instances. When this compromise is 
made, the resultant simplified network becomes that shown in 
Figure 2-3. 








Heat Pipe Limitations 

Limitations on the heat transfer capability of heat pipes occur 
due to effects that inhibit vapor or liquid flow. The deriva- 
tions of the relations governing these limits are covered ade- 
quately in the literature (Ref. 2,5) and will only be briefly 
summarized here. The effects to be dealt with include the hydro- 
dynamic, sonic, entrainment and nucleation limits* The limita- 
tions discussed below are computed in subroutine LIMITS. The 
inputs to subroutine LIMITS consist of the first data array number, 
evaporation conductance and temperature required for the boiling 
limit calculation, physical constants, and working fluid pressure vs time 
array number. Several limitations are placed upon the use of MITAS 
constants and arrays, primarily that The arrays are 
input in a specific order in a block. This routine will be 
generalized in the future to remove the current limitations. 

a. Hydrodynamic Limit 

The hydrodynamic or “wicking** limit is a result of the pres- 
sure variations due to moving fluids within the heat pipe. 

A circulation is established in the pipe fluids due to the 
capillary pumping capability of the wick. When the pres- 
sure drops due to viscous and momentiam effects become large 
enough to overcome the capillary pumping pressure rise plus 
the effects of body forces, the hydrodynamic limit has been 
reached. The liquid being vaporized cannot be replaced 
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through capillary action at a high enough rate, resulting 
in evaporator dryout (burnout) and a corresponding evapor* 
ator temperature Increase* Thus 


AP 

c 

+ 


m 

r 

< 

L 

+ 

AP 

V 

Capillary 

Head 

+ 

Body 

Force 

Head 

m 

Liquid 

Pressure 

Drop 


Vapor 

Pressure 

Drop 


Each of these contributions will be discussed In the follow- 
ing paragraphs. 

The capillary pumping pressure rise is due to the difference 
In meniscus radius of curvature In the evaporator and con- 
denser wicks. Vaporization In the evaporator results In a 
concave meniscus tending to recede Into the wick. Conden- 
sation In the condenser results In a nearly flat meniscus 
at the wick surface. The surface tension forces attempt to 
flatten the evaporator meniscus, creating a pressure dif- 
ference which pumps liquid toward the evaporator. The 
maximum AP which can be generated due to capillary action 
Is 


AP 

c 


2 g cos 'I' 

R ’ 

pore 


(2-29) 


where Rpore ^ minimum effective pore radius which cor- 
responds to a minimum evaporator meniscus radius, and Is 
usually experimentally determined. ^ Is the wetting angle 
as defined In Figure 2-4. 


i 


i 







I 








The body force contribution is a result of acceleration 
fields in the pipe environment, such as gravity, thrust, 
or rotation. These fields may either augment or diminish 
the performance of the pipe, depending on its orientation 
with respect to the field. The body force pressure drop 
AP^ can be expressed as 

APb = (+ Lcos /3 + sin/3), (2-30) 

p 

where a is the body force acceleration and fi is the angle 
between the body force vector and the pipe axis, as defined 
in Figure 2-5. 



Figure 2-5 Definition of Pipe Orientation Angle g 


Equation 2-30 indicates that the body force pressure drop_ 
is made up of two components, that of the pipe length paral- 
lel to the body force vector, and that of the internal pipe 
diameter parallel to the body force vector. That is, the 
capillary action must raise the fluid a height equivalent 
to the diameter of the pipe if the pipe is horizontal in a 
vertical acceleration field, the height of the pipe if it 
is vertical (and the evaporator is at the top) , or some 
intermediate distance if the orientation is other than 
horizontal or vertical. The (+) in equation 2-30 indicates 
that the length component can be positive or negative de- 
pending on whether the evaporator or the condenser is raised. 

A positive length component indicates a pressure drop; 
negative indicates a pressure rise aiding the capillary action. 
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The third contribution to the pressure balance is the pres- 
sure drop due to liquid viscous flow losses, AP^ . From 
Navier-Stokes for steady, incompressible, fully-developed 
laminar flow in a rectangular cross-section passage, the 
axial pressure gradient dP^/dz is 


dz 


- 






^ 2 


^ X c/ y 

or, in cylindrical coordinates, 


(2-31) 


DP 


1 


dz 





^ u 


1 




JIL 


^ 9 “^ 


(2-32) 


These equations must be solved individually for each geo- 
metry involved. The solutions for five cases — open grooves, 
closed grooves, closed cylinders, full annuli, and partial 
annuli — have been derived for :the model. The solutions of 
equations 2-31 and 2-32 are all of the following forni: 






ra L 


(3600) K A 


(2-33) 


xs 


which, after substituting equation 2-8, becomes 


AP,= 


(3600)^ g^ hfg K A 


(2-34) 


xs 


where KAxs is the (permeability x cross-sectional area) product 

obtained from the solution of Navier-Stokes 

for the appropriate geometryo These KA^s relations have 

been programmed into subroutine PERMBL, which computes 

KA^q for the input geometry. PERMBL is discussed in 

Section 2-A-6* 

The fourth contribution to the pressure balance in equation 
2-28 is the vapor pressure drop due to evaporation, con- 
densation and flow losses detailed in the conductor deriva- 
tions of Section 2-4. The evaporation AP, and the con- 
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densation AP at low suction rates, is given by equation 
2-12 with 2-13 used to eliminate Cf. The condensation 
AP at high suction rates is given by 2-20 with 2-21 used 
to eliminate dP/dz, The vapor pressure drop due to flow 
losses is given by equation 2-25. 

These four contributions to the heat pipe pressure balance 
are inserted into equation 2-28, which is then solved for 
Q|^, the maximum hydrodynamic latent heat transfer capability. 
Since the suction or injection Reynolds number, Re^s is 
dependent upon Q (equation 2-15), an iterative solution 
is the most convenient method of solving 2-28. , 

Sonic Limit 


The second limitation on heat transfer capability in a heat 
pipe is the sonic limit. This limit arises since constant- 
area flow with suction and injection is analogous to con- 
stant mass flow in a converging-diverging nozzle (Ref, 5), 
The end of the injection (evaporator) section of a heat pipe 
is analogous to the throat of the nozzle, and consequently, 
there is a flow velocity limit at the evaporator exit ana- 
logous to choked flow (Mach 1) in the nozzle throat. As 
the heat input is increased above the sonic limit, an 
evaporator temperature increase results since no more axial 
heat transfer can occur. The definition of the Mach number 
gives 
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where U is the bulk vapor velocity and c is the sonic vapor 
velocity. Using 


U = 


m 


P A 
W xs 


(2-36) 


V 

and equation 2-8, and solving for Q, equation 2-35 becomes 


Q = p hfg A M 
^ xs 


Jyp (3600)^ g (144)/ P^ 
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where the cross-sectional vapor flow area Aj^g is equal to 
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The maximum velocity through the '^throat,** or 


evaporator exit, occurs at M « Therefore the maximum 
sonic latent heat transfer capability Qg is 


n D . 
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Entrainment Limit 

The third heat transfer limit is the entrainment limit and 
results from a high velocity vapor stream entraining the 
liquid in the wick, preventing it from returning to the 
evaporator, resulting in evaporator dryout and a temperature 
increase* The capability of the vapor to entrain liquid is 
a function of the wick geometry, e*gp, open axial grooves 
are far more subject to entrainment than a tight mesh screen. 

The force Fq required to tear a particle of liquid from the 
wick is proportional to the product of the dynamic pressure 
and a characteristic length X (representing wick geometry) 
squared: 

Pv 2 

F cC — i X . (2-39) 

(3600)^ g 


The restoring force required to retain that particle is 
proportional to the surface tension and the characteristic 
length; 

^ o- X . (2-40) 

The ratio of these two forces is the Weber number 

p X 

Wb = ^^—2 • (2-41) 

(3600) g^cr 

When the Weber number exceeds unity, entrainment will occur 
and pipe performance degradation will begin* Although* the 
entrainment limit has not actually been reached at Wb = 1, 
the true limit is difficult to compute. Therefore, the con- 
servative approach is to assume it has and to set Wb == 1 in 
equation 2-41. Using 2-8 and 2-36 to eliminate U and 
solving for Q, equation 2-41 becomes 
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Nucleation Limit 


(2-42) 


The fourth heat transfer limit is the nucleation or boiling 
limit and is caused by superheating of the liquid at the 
pipe wall. Since the liquid and vapor are at equilibrium 
at the wick surface, a temperature drop must exist across 
the wick if any heat is to be transferred, and this AT 
increases with increasing loado At some critical AT, 
nucleation will occur near the point of highest superheat, 
the pipe wall. This bubble inhibits replacement of the 
evaporating liquid, causing dryout and a temperature increase. 


For a bubble of radius at equilibrium, the expansion 
force must equal the restoring force, i.e.. 
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where is the saturation pressure at the wall temperature, 
^wall> is the saturation pressure at the wick surface 
temperature^ l^wick* 

Solving equation 2-43 for - Pwick? 
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wick 



(2-44) 


For small P n - P . i , 
wall wick’ 


P -- - P . , hfg p H 

wall wick ^ ^ _ V 

T - . - T . , dT ” T 
wall wick V 


from equation 2-9. Substituting 2-45 into 2-44, 
2 cr T 

*^wall "^wick hfg p R H 

V n 


(2-45) 
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If this critical AT is reached, the bubble will grow and 
nucleation failure will result. The corresponding maximum 
nucleation latent heat transfer capability is 
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For homogenous wicks, Figure IV-6a, 



- (1 - g„) (V 
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(2-51) 


from Reference lo For open axial grooves. Figure 2-6b, k^/v 
is computed assuming two parallel heat flow paths, one through 
the liquid and one through the lands between grooves o The lands 
are assumed to be composed of wall materials 
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where 
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(2-53) 


or the total mean circumference less the total width of the 
grooves multiplied by the pipe length, and 

A. = N WL. (2-54) 

g 

Substituting equations 2-53 and 2-54 into 2-52, 
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Screen-covered axial grooves are treated as a homogenous wick 
(equation 2-51) and open grooves (equation 2-55) in series, 
see Figure 2-6c: 
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The fourth wick type is screen-covered axial spacer wires as 
shown in Figure 2-6d. This type is treated similarly to screen- 
covered axial grooves, except the path corresponding to the wall 
equation 2-52 must be 'handled differently to account for the 




REPRODUCffilUTY OF TEIB 
ORIGINAL PAGE IB POOR 


poor contacts between wire and pipe wall and between wire and 
screen. The method chosen to simulate this poor contact Is to 
treat the wire as three materials In series: a layer of liquid, 

one of wire and another of liquid as shown In Figure 2-7. 


Liquid. 


Wall 

—Wire 


Screen 




Figure 2-7 Simulation of Spacer Wire Contact Conductance ^Liquid 


The liquid layer thicknesses are analytically determined as a 
first approximation, but test data Is probably necessary for an 
accurate The assumption made In computing the liquid 

layer thickness Is that the volume of liquid under the spacer 
wire to the left of line C-C is equal to that to the right of 
C-C in Figure 2-8. 
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Figure 2-8 Liquid Layer Equivalent Thickneea 
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This is analogous to assuming that the same amount of heat that 
has a less difficult path through the thin portions of the liquid 
layer, has a more difficult path through the thicker portions of 
the liquid layer. The liquid layer thickness is obtained in 
Appendix B and is equal to .1738 Rspacer? where Rspacer 
spacer wire radius. The conductance through the spacer wire is 
then 
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The conductance through the liquid between spacer wires is 



spacer 



D ) - (N D )" 

spacer spacer spacer 


Lo 


(2-58) 


The conductance through the screen covering is 
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where y^screeu ^he homogeneous wick effective thermal conduc- 
tivity from equation 2-51. The three contributions are combined 
as follows: 
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and the corresponding k^^^ 

V, ^ c = r 

//w total A total 


is computed from 
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The fifth wick type is the circumferential groove (or threaded 
pipe) as shown in Figure 2-6e. This case is treated as one con- 
tinuous groove and equation 2-52 is used with the following 
values for A. 
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The inputs to subroutine EFFKWK are the type of wick, k^, kwall> 

^ 7 ? ^spacer? ^spacer? 

use is restricted by the assumptions made above for each wick. 

The second auxiliary subroutine is called PERMBL and is used to 
compute the product of permeability and cross-sectional area for 
liquid flow for the five wick types in Figure 2-6. These KA^g 
values are obtained from Appendix A by comparing the specific 
liquid pressure drop solutions contained there to the general 
form of the pressure drop 


“xs ‘•a 


(2-65) 


The KAj^g values needed to satisfy equation 2-65 for each wick 
type are as follows. For homogeneous wicks (Figure 2-6a) : 
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since the permeability of a homogeneous wick is generally known. 
For open axial grooves (Figure 2-6b) , 
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xs g 
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For screen-covered axial grooves (Figure 2-6c) , 


KA = N 
xs g 
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ilcreen-covered axial spacer wires (Figure 2-6d) , are treated 
as annular spaces: 
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for each annular space, where 
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Utility subroutine ANNSPG is used to compute B. 


The values of the outer annulus radius , the inner annulus 

°a 

radius ^i^^ snd the half-angle between spacer wires are computed 
as follows: 
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(2-71) 


The half angle between spacer wires 8 is derived from % :c length y ^ 

number of space^ as expressed in equation 2-71. The total KA^s 
is then the number of annular spaces multiplied by the KA^s 
in equation 2-69. 


Circumferential grooves are treated as one long open groove and 
equation 2-67 applies with the number of grooves Ng = 1. How- 
ever, the pipe length L in equation 2-65 must now be replaced 
by the actual length of the groove 
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such that KA^g is expressed as 
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(2-73) 


Ihe inputs to subroutine PERMBL are the type of wick, 

Ng, W, d, K, 7T, ^spacer Ngpacer* Its restrictions correspond 
to the assumptions made for each wick type. 
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I Another utility routine was written to automate the evaluation 

;i of specific properties which are mono- variant . Subroutine 

[j LOOKUP provides a *’D0*' loop table lookup using the library 

routine DIDEGl. 

f. • 

i; The calling sequence is 

LOOKUP (N,T,A,P) 

i. i 

i- 

^ where 

y N = number of lookups in this call 

' T = independent variable 

'! A = Array of T, P values 

U P = Dependent variable (property) 






GOVERNING EQUATIONS FOR FIXED -CONDUCTANCE ARTERIAL PIPES 

An artery is inserted in the vapor flow passag;e of a heat pipe 
to provide an atixiliary flow passage for the liquid returning to 
the evaporator. By removing most of the liquid flow from the 
pipe wall vicinity, large radial conductances can be achieved, 
minimizing the temperature drops across the wick. At the same 
time, the artery can be designed for minimum pressure drop, 
maximizing the hydrodynamic limit. The result is a decoupling 
of hydrodynamics and heat transfer in the pipe, with correspond- 
ing performance increases. Typical geometries are the pedestal 
artery and its high performance derivative, the spiral artery, 
shown in Figure 2-9- 
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Figure 2-9 Available Artery Geometries 


Capacitance 


Artery capacitance is computed in subroutine ARTERY and is added 
to the liquid/wick capacitance calculated in equation 2-2. 

Artery supports are ignored. The pedestal artery shown in 




Figure 2-Sa has a capacitance equal to the sum of the screen 
capacitance and the liquid capacitance: 




screen 


- 



L (1 - 0 J p , C 
^ art' ^art p 

am 


-d2 

1 


art 


L (1 - 0 J 
' art 


(2-74) 


(2-75) 


The liquid volume is the total artery volume less the screen 
volume. Summing equations 2-74 and 2-75, the total artery 
capacitance C^^t given by 
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screen 
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The second artery configuration is the spiral type, shown in 
Figure 2-9b. Its capacitance is the sum of the screen, spacer 
wires, and liquid capacitances* The screen length is obtained 
from an integration of the equation of a spiral over the total 
angular displacement of the spiral, see Appendix C. The screen 
length is thus 
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where 
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The screen capacitance Cg^reen then the sum of the spiral screen 
volume and the wrapping screen volume, all multiplied by the 
density and specific heat: 
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The spacer wire capacitance is 
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since there is one more spacer wire than there are annular spaces. 
The liquid capacitance is 
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since the liquid volume is the total artery volume less the screen 
and spacer volumes. The total spiral artery capacitance (^art is 
then 


^art screen spacer ^ * 
Axial Conductance 


(2-82) 


The artery contributes to the total axial heat transfer capability 
of the pipe. Its axial conductance is calculated in subroutine 
ARTERY and is added to the liquid/wick conductance from equation 
2-4. The pedestal artery conductance is equal to the sum of the 
screen conductance and the liquid conductance. The screen con- 
ductance is “ 
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where the screen volume fraction (1 - halved since con- 

tact resistance between axial and circumferential results in 
only the axial members of the screen conducting heato The liquid 
conductance is 
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The total conductance of the pedestal artery is then 
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The spiral artery conductance is equal to the sum of the screen, 
spacer wire, and liquid conductances. The screen conductance is 
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where the screen volume fraction (1 - 0art) is again halved since 
only the axial screen members carry heat« The spacer wire con- 
ductance is 


k /D" ^ 

spacer L a \ 4 y 


(2-87) 


and the liquid conductance is 
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The total conductance is thus 
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for the spiral artery configuration. These artery conductances 
are computed for each pipe section** When they are combined with* 



the liquid/wick axial conductance, they are converted to a node 
center-to-center conductance by 


3 . 


4. 


G 


art 
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and then added to the liquid/wick axial conductance linearly. 
Radial Conductance 


The arterial radial conductance is not considered since no radial 
heat transfer occurs across an artery. 

Arterial Effects on Hydrodynamics 

Since the purpose of an artery is to improve the hydrodynamic 
limit of the heat pipe, the liquid pressure drop computation of 
Section 2-A-5 must be revised. The (permeability x cross-sectional 
area) KA^^g must be computed for the artery and added to that of 
the wick. This calculation is performed in subroutine PERMBL, 
using KA^s relations derived in Appendix A and summarized here. 


The pedestal artery contains simple Poiseuille flow and ti 
meability cross-sectional flow area product is thus 
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The spiral artery is extremely complex and must be treated as 

multiple annular spaces with varying internal and external radii 

(Ri and ) and differing angular widths ( S ) . The basic re- 
a a 

lations for the KAxs of an annular space are equations IV-69 and 

IV-70. The values of Rq , Ri > and S, however, require modi- 

a a 

fication for use in a spiraling configuration. The method used 
is a ”D0" loop with initialized values of Ro^^? 8 cor- 

responding to the innermost annular passage. The radii are in- 
cremented using the equation of a spiral to the values correspond- 
ing to the next outward space, and 8 is adjusted if the wire 
spacing changes. KA^g values are computed for each space and 
summed. 
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The general equation of a 
r ■ m9 + b. 

Referring to Figure 2-10, 
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Figure 2-10 Detailed Spiral Artery Geometry 


Therefore, the equation 
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of the spiral artery is 



+ t. 

and 8 are 


(2-93) 








The subtracted term in the 8 computation of equation 2-94 
accounts for the angular contribution of the spacer wire. 

The loop iterates on j from 2 to where is the number 

of annular spaces. The general relations for Rq s % 5 and 8 
in the ’IDO" loop are ^ ^ 



R = R. 4- D , 

o^ spacer 
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where R. is the R. from the previous iteration and the additive 
a a 

term in the R. computation of equation 2-95 represents the AR. 

t a I 

to convert R. to the new R. , Utility subroutine ANNSPC is 
^a ^a 

called in each iteration of the loop, and the resulting KA^g 
values for each space are summed into a total KA^g, 


There is an additional passage in the center of the spiral which 
is treated as Poiseuille flow: 
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This contribution is added to the summed annular KA^g values for 
a total spiral artery KAxsj which in turn is added to the wick 
KA^^g to compute the total liquid pressure drop. The additional 
inputs to PERMBL to account for arteries are artery type, Dq .5 

^^art’ ^wrappings array of 8 values, N^, Nt, Dgpg^Q 03 -, t, and 

The second major effect of the artery on the hydrodynamic limit 
is a reduction in the vapor flow cross-sectional area. A re- 
vised effective vapor flow passage diameter is computed 

in subroutine ARTERY as: 
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art 


and this value is used in place of Dy in the vapor pressure drop 
equations 2-12, 2-20 and 2-25 when an artery is present. 

Subroutine ARTERY requires the artery type, Dq . , D± , 

SI u o j7u 

^wrapping j ^spacer > ^art? ^t» ^^Part’ ^^Pspacer’ ^art> 

kspacerj and ^^ip^^gtal* 

The third effect on heat pipe hydrodynamics concerns the fact 
that the largest pore in the wick/ artery system determines the 
capillary pumping capability. If the artery contains a pore 
with a radius larger than Rpore equation 2-29, the artery 
Rpore ^^st be substituted. A decrease in the hydrodynamic limit 
results o 

Artery Self-Priming Limit 

A new pipe heat transfer limit arises when an artery is intro- 
duced. Since an artery can be depleted of liquid by forces 
arising from acceleration fields, its design must allow repriming. 
The artery will reprime if the capillary head generated by the 
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largest gap in the artery exceeds the body force head, under zero 
loado The existence of a load, however, introduces the additional 
head losses due to flow which must be accounted for. The maximum 
load under which the artery will reprime can be computed using 
equation 2-28, Where AP^ is calculated using the maximum artery 
gap width (excluding the end of the condenser since no liquid/ 
vapor pressure difference exists there) instead of the wick pore 
size Rpore in equation 2-29. This is accomplished by repeating 
the hydrodynamic limit logic in subroutine LIMITS with the 
artery gap size replacing the wick pore size, 

GOVERNING EQUATIONS FOR VARIABLE -CONDUCTANCE PIPES 


The variable-conductance heat pipe (VCHP) utilizes a charge of 
non-condensible gas added to the working fluid to achieve a 
relatively constant heat pipe temperature regardless of heat load. 
This gas is swept to the condenser end of the heat pipe and forms 
a plug which blocks the flow of vapor. (Note the distinction 
between non-condensible gas and working fluid vapor . ) As the 
pipe temperature drops, the vapor pressure of the working fluid 
drops, causing an expansion of the non-condensible gas plug. 

This expanded plug effectively shuts off a portion of the con- 
denser to vapor flow, decreasing the area available for condenser 
heat rejection. This, in turn, increases the pipe temperature. 

In a properly designed VCHP, the pipe temperature can thus be 
made nearly independent of heat load. A reservoir for the non- 
condensible gas at the condenser end of the h,eat pipe provides 
closer temperature control since a greater vapor/ gas interface 
movement can be achieved for a given" temperature (and pressure) 
change. Temperature control of the resertyoir through a heater 
can be used to actively control pipe temperatures by a feedback 
loop, A schematic of the variable-conductance heat pipe is given 
in Figure 2-11. 

The simulation of non-condensible gas control of the operating 
temperature of a heat pipe regardless of load requires more 
sophisticated modeling techniques. The location of the vapor/ gas 
interface ("front”) in a variable-conductance heat pipe must be 
computed based on ideal gas relations. The presence of a gas 
reservoir to provide improved control adds complexity to the 
problem, as does the popular use of a reservoir heater to 
achieve even more precise control. 

The basic assumption in the following analysis is that the effect 
of the diffusion of vapor into the gas plug on the temperature 
gradient is negligible compared to the effect of axial conduction 
in the pipe wall on the temperature gradient. This is generally 
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a good assumption for .high thermal conductivity pipe walls. The 
number of moles' of tipn'^.pondensible gas is a constant n: 


- P V. 


R T 
I. gas 


(2-98) 


summed over the volume elements of the reservoir and inactive 
condenser* There is no gas in the active section of the condenser, 
the adiabatic section, or the evaporator due to the tendency of 
the moving vapor to sweep it towards the condenser* In the 
reservoir, 

T = T 
gas res * 


and 




(2-99) 


P = P 
gas V 


act 
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res 




since l^v^ct the total pipe pressure, and the sum of the vapor 

and gas partial pressures in the reservoir must equal the total 
pressure* 


Equation 2-99 assumes a negligible pressure gradient in the 
pipe. It also assumes the reservoir vapor temperature is equal 
to the reservoir temperature, and that the vapor partial pressure 
is that of saturated vapor at the reservoir temperature. In the 
jth element of the inactive condenser. 


T =T - ■ , 
gas. wall. 

J J 


and 

P = P, " P 
gas . V v . 

^ j act wall. 

J 




> 


( 2 - 100 ) . 


The assumptions for equation 2-99 are used in 2-100 also, sub- 
stituting the condenser wall temperature for the reservoir temp- 
erature. Inserting equations 2-99 and 2-100 into 2-98, 
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where the volume of the jth element of inactive condenser is the 
total condenser volume less the active condenser volume, divided 
by the number of inactive elements J. Dividing by using 

the fact that V^^t/'^cond equal to Laj.j./Lcond ^ constant- 
diameter condenser, and solving for 


“'act = 1 + 
■“cond 
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wall 


Multiple vapor-wall couplings are needed to simulate a variable- 
conductance heat pipe* The couplings in the inactive region are 
zeroed according to the vapor/ gas interface position determined 
in equation 2-102, to simulate no heat pipe action in that 
region. 


This is accomplished by a call to subroutine ADJG. Condenser nodes ^ 
which are within the active length are left alone while the con- 
ductances are set to zero for those nodes which are totally beyond 
the active length. The conductance of the node in which the front 
is located is adjusted based upon the proportion of active length 
calculated. 


The calling sequence is 


ADJG (NGG,NC,LC,LAGT,GCOKD) 
where 

NCG = number of moles of control gas for a VGHF. 
0.0 for a fixed conduction pipe 

MG = number of condenser nodes 

LG = condenser nodal length 

LACT - active condenser length 

GCOND = first condenser wall-to-vapor coupling 
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Since dependent on pipe temperatures and pressures 

and since the temperatures and pressures are dependent on the front, 
an iterative solution is necessary to achieve a satisfactory result. 

In steady-state problems, the iterative solution is handled by multiple 
calls to STDSTL from EXECTN since VAKBLl is called each iteration of 
the network solution. In transient problems, VAkBLl is entered only 
once for each time step; consequently, the BACKUP option is used in 
VARBL2 to force a recomputation of temperatures if the ^^^^^^cond 

based on the most recent network solution does not match the ^^^ti/^^ond 


from the previous solution 
using Pv^^^, R, T 


L /L - is computed in subroutine FLFRNT, 
act cond ^ 


res' 


V 


res- 


’cond. "■ 


J, Pv and T as 

’ wall wall 


inputs . 


Subroutine FRONT controls and dampens the movement of the vapor/control 
gas interface location resulting from subroutine FLFRNT. The results 
of FLFRNT predict the new front location based upon temperatures re- 
sulting from the present time step and the previous front location. 

The max im um front movement allowed is one condenser node length. 

Since the new front location will result in new temperatures, and 
since FLFRNT will predict a new front location based on these new 
temperatures, an iterative process is required. This iterative 
process io repeated until the front location change is reduced to 
1% of the condenser node length or less. At that point the iteration 
is completed and the time step is advanced. Additional logic is 
required in the execution block for steady state runs and in the 
variables 2 block for transient runs to support subroutine FRONT. 

Refer to the program listing for the supporting logic. 


The calling sequence is 


FRONT 


No arginnent list is included for ths subroutine. 

Computation of the last condenser node axial conductance to the 
reservoir is obtained by calling subroutine RSCOND for a VCHP. 
The same assumptions used in subroutine AXCOND also apply. The 
wicks used in the condenser, feedtube and reservoir are assumed 
to be identical. The following network is reduced to a single 
conductor. 
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REPRODUCBILITy OF THE 
OEIGINAL PAGE IS POOR 


... 



Each of the three parallel combinations are simmed and are then 
add by the series combination rule. The artery is not considered 
in this calculation since it does not extend beyond the condenser. 


The calling sequence is 

RSCOND (NCG, PI, LWA, KWA, ODWA, IDWA, LET, KFT, ODR, IDR, LR, 
KR, ODFT, LWK, POROS, KWK, ODWK, IDWK, IDFTW, IDRW, GFTTOT) 


NCG 

LWA 

KWA 

ODWA, IDWA 

LFT 

KFT 


ODR, IDR 

LR 

KR 

ODFT, ID FT 
LWK 
POROS 
KWK 

ODWK, IDWK 
IDFTW 
IDRW 
GFTTOT 


FLAG FOR VCHP 

CONDENSER WALL NODAL LENGTH 

CONDENSER WALL THERMAL CONDUCTIVITY 

CONDENSER WALL OUTER AND INNER DIAMETERS 

FEED TUBE NODAL LENGTH 

FEED TUBE THERMAL CONDUCTIVITY 

RESERVOIR OUTER AND INNER DIAMETERS 

RESERVOIR LENGTH 

RESERVOIR THERMAL CONDUCTIVITY 

FEED TUBE OUTER AND INNER DIAMETERS 

CONDENSER WICK NODAL LENGTH 

WICK POROSITY 

WICK THERMAL CONDUCTIVITY 

CONDENSER WICK OUTER AND INNER DIAMETERS 

FEED TUBE WICK INNER DIAMETER (ODFTW = IDFT) 

RESERVOIR WICK INNER DIAMETER (ODRW = IDR) 

TOTAL AXIAL FEED TUBE CONDUCTANCE 
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GLOSSARY 


A 

B 

c 

G 

P 

P 

F 

G 

H 

J 

K 

L 

M 

N 

P 

Q 

R 

R 

Re 

Re 


area 




w 


T 

U 


defined by equation IV-70 

capacitance (mass x specific heat) j^Btu/°Fj 
specific heat j^Btu/lbra-°Fj 
diameter j^ftj 
force [ibfj 

conductance j^Btu/hr-°Fj 

mechanical equivalent of heat j^t-lbf/BtuJ 

number of inactive condenser sections j^dimensionlessj 

wick permeability 

pipe length [“] 




Mach number 1 dimensionlessj 


number of grooves, turns, ^spacers, etc, (depending on 
subscript) j^imens ionless J 

absolute pressure j^lbf /ft 

latent heat transfer rate |^tu/hr^ ^ 

radius [f 

universal gas constant [l545 ft-lbf/lbmole~°R] 

Reynolds number along tube j^dimensionless^ 

Reynolds number for suction (-) or injection (4*) at 
wick surface j^imensionlessj 

temperature M 

bulk average vapor ' elocity j^t/hr^ 


C-112 




V 

W 

Wb 

X 

a 

c 

c^ 


volume [ft^J 


hfg 

k 

1 

• 

m 

n 

r, 9 

t 

u 

V, w 

X, y 
z 

r 

A 


groove width [ft^ 

Weber number jdimensionlessj 
entrainment characteristic length ‘ 

body force acceleration j^ft/sec^J 
sonic velocity j^ft/hrj 

skin friction coefficient [^dimensionles^ 
groove depth |^ft^ 

conversion constant relating force and mass |32,2 Ibm- 
ft/lbf-secM 

heat of vaporization of working fluid jjBtu/lbmJ 
thermal conductivity j^Btu/hr-ft-^pj 
distance along fin length [ftj 
mass flow rate J^lbm/hrJ 

number of moles of non-condensible gas j^lbmol^ 
polar coordinates 

spiral artery screen thickness j^ftj 

local velocity [^t/hrj 

defined by equation IV-78 

rectangular coordinates 

axial distance along pipe |^ft^ 

pipe orientation angle j^degreesj 

specific heats ratio (Cp/Cv) [^dimensionless^ 

difference 
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8 


<!> 


half-angle width of annular space j^adiansj 

dynamic viscosity [^Ibm/ft-hrJ 

3.14159 

density Q.bm/ft^J 
surface tension [^bf/ft^ 

wick porosity (void fraction) j^imensionlessj 
wetting angle [^egrees^ 


Subscripts 


a 

annular spaces 

act 

active 

art 

artery 

ax 

axial 

b 

body force 

c 

capillary 

cond 

condenser 

e 

entrainment limit 

eff 

effective 

enth 

enthalpy 

evap 

evaporator 

f 

fin 

g 

grooves 

h 

hydrodynamic limit 


i 


inner 


iteration counter 


liquid 

nucleation limit 

outer 

pipe 

restoring 
radial 
rejection 
reservoir 
fin root 
sonic limit 

turns of spiral artery 

vapor 

wick 

cross-sectional 


pipe sections 
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ATTACHMENT 3 

MARTIN MARIETTA CORPORATION 
SPACE RADIATOR PROGRAM 
FIN GRADIENT TECHNIQUE 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


As illustrated in Figure there are essentially three types of gradients 

which are possible in a space radiator. These are: 

a) The classical condition where the adjacent fin root temperatures 
are equal giving a symmetrical temperature distribution between 
nodes. The point where the temperature gradient, dT/dX, is equal 
to zero exists at the midpoint. This condition could have been 
analyzed by the fin effectiveness technique which is commonly used. 
(Ref. 3~1) 

b) The condition where adjacent fin root temperatures are not equal 
giving a nonsymmetrical temperature distribution and the point where 
the temperature gradient dT/dX is equal to zero exists at some point 
other than the midpoint. 

c) The condition where the difference between the adjacent fin root 
temperatures is so large that no point exists between nodes where 
the temperature gradient dT/dX is equal to zero. 

The fin gradient technique used in the MMC radiator program has the capa- 
bility to evaluate all three types of gradients. 

Referring to Figure D-2, the heat which is radiated from an element of the 
fin surface is as follows: 

d ' = € a W dX (t'^ - T (1) 

^rad ^ 


and the heat which is transferred in the fin by conduction is: 


= V AT 
/^cond dX 

Differentiating equation (2) with respect to X, 

dq. 


( 2 ) 


cond _ ,, d ,dTx 

lx lit <lii> 


(3) 


and since the heat lost by radiation represents the difference in heat 
passing across the two internal faces of the element, then equations (1) 
and (2) can be equated as follows: 


^KAD 


= -dq 


cond 


Performing the required substitutions and differentiating yields the follow- 
ing differential equation: 



R 



( 4 ) 


with the boundary conditions 8 = 1 (3 Z = 0 

. 0 = 0„ @ Z = 2 
K 
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where : 


(sometimes referred to as the fin profile 
number) 




3 2 

ex L 


«R - ^r+l'*r 


9 = T /T 

e e r 


T = 
e 


I f A T^+Q 1 ^ 

n n n ^abs — 


Ratio of adjacent fin root temperature to 
the root temperature 

Ratio of equivalent temperature of surrounding 
space to root temperature 

Environmental temperature 


9 = T/T^ 
Z = X/L 


■ d9/dZ = d(T/Tp/d(X/L) 

a = Fin cross sectional area 

T ~ Fin Temperature 

<r = Stephan Boltzmann constant, 

£ == Radiator surface node emissivity 
= Fin root temperature 
T^+1 = Root temperature of adjacent fin 
= Fin length to right side of fin 
~ Fin length on left side of fin 
K = Fin thermal conductivity 
y = Fin thickness 

L = Half the distance between tubes 
= Area of radiator surface node 

jr = Gray Body view factor from other external surface to the radiator 
skin node 

A == Area of other external surface 
n 

= Temperature of other external surface * 


Q, = Absorbed flux from solar, albedo, and earth IR sources 
^abs 


W = Fin width 


A set of tables was available (Ref. D-2) which expressed the dimensionless 
temperature gradient d9/dZ as a function of the three variables Nr, 0r, and 9^, 
These tables were expanded by MMC to also include the condition where the 
equivalent temperature of surrounding space was hotter than the root temper- 
ature ( 9 q 1). This was required for the AM radiator since the segments of 

the radiator which were exposed to the sun actually reheated the coolant as 
it flowed through the radiator. 


The net heat to skin node was determined in the radiator model 

for each time step by a table look-up for the d9/dZ on the right side and the 
left side of the fin root. Then by multiplying these dimensionless gradients 
by Tj;/L, the temperature gradients dTR_/dXR and dT^/dX^ were then determined for 
the right and left side. The net heat loss from the root of the radiator fin 
taken at the stringer was then determined by performing the following energy 
balance on the fin root. 


q ^ = -Ka 
^net 


dT dT 

b ^ R A 

•dx^ dx^ 

J_t K 


(5) 


This technique was included in the subroutine RADFIN which calculated the 
total heat loss for each radiator skin node. The heat loss was first calculated 
on both sides of the tube by using the node temperature at the stringer attach . 
points as the root temperatures. The heat loss was then calculated for both 
ends of the node by assuming that the temperature gradient at the edge was nearly 
adiabatic with 9 r 1.0 and the root temperature was the node temperature at 
the stringer attach point. The subroutine also computed the environmental 
temperature which was printed later for reference. 

To determine the significance of the circumferential conduction, a com- 
parison was made between the two cases, with and without circumferential con- 
duction. As a result, neglecting circumferential conduction gave approximately 
107o better performance than the more accurate method. A side conclusion would 
be that the radiator performance could be improved if the hot side of the radiator 
were isolated from the cold side. Because of the nonlinear nature of radiant 
interchange, this would result in a higher heat rejection from the surface. 
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